DFPD 01/TH/30 



THE PROBLEM OF NEUTRINO MASSES 
O ' IN EXTENSIONS OF THE STANDARD MODEL f] 

o: 

3 : ISABELLA MASINA 

P^ ' Department of Physics G. Galilei, Univ. of Padova 

Via Marzolo 8, 1-35131 Padova, Italy R 

>: 
o. 

(NT 

o. 

O' 

We review the problem of neutrino masses and mixings in the context of Grand Uni- 
fied Theories. After a brief summary of the present experimental status of neutrino 
physics, we describe how the see-saw mechanism can automatically account for the 
large atmospheric mixing angle. We provide two specific examples where this pos- 
sibility is realized by means of a flavour symmetry. We then review in some detail 

.^ \ the various severe problems which plague minimal GUT models (like the doublet- 

triplet splitting and proton-decay) and which force to investigate the possibility of 
constructing more elaborate but realistic models. We then show an example of a 
quasi-realistic SUSY SU(5) model which, by exploiting the crucial presence of an 
abelian flavour symmetry, does not require any fine-tuning and predicts a satis- 
factory phenomenology with respect to coupling unification, fermion masses and 
mixings and bounds from proton decay. 
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Introduction 



It is quite obvious that the Standard Model must be extended. Various arguments 
bring along the important conceptual implication that the Standard Model should 
be seen as a low energy effective theory valid up to some physical cut-off scale. 

It is thus compelling to speculate about the underlying more fundamental the- 
ory. Various ideas have been proposed and developed to answer this question. For 
instance, recently a lot of interest and efforts have been devoted to the idea of 
extra-dimensions, where the different behaviour of the gauge and gravitational in- 
teractions may provide a solution to the hierarchy problem. But the attractive 
proposal of Grand Unification has not lost its charm. 

The concept that all particle interactions unify at very high energies is indeed not 
very recent. The first paper on this subject was published 26 years ago and from 
then on this idea captured the interest of physicists. Many fundamental papers 
that explore the various aspects and consequences of Grand Unification have been 
written. What is then inside this proposal which renders it still so interesting to 
collect the general interest? The main point is that it received phenomenological 
support by: 

A) the success of gauge coupling unification; 

B) recent developments on neutrinos. 

A) Gauge coupling unification works terribly well. The original proposal has 
however notably evolved because now we know that Grand Unification needs Super- 
symmetry Only then gauge coupling unification successfully happens at the scale 
Mqut ~ 10 16 GeV. In addition, in non supersymmetric Grand Unified theories the 
prediction of proton lifetime for the process p — > e + ir° is in conflict with lower ex- 
perimental limits, but including Supersymmetry it turns out to be close to the safe 
side. 

B) Recent data coming from neutrino physics are much more robust and less 
fragmented then they were just few years ago. There is definitive evidence for a 
neutrino mass associated with atmospheric oscillations in the range 0.05 -r- 1 eV. 
There is no room in the Standard Model for neutrinos to have a mass, so that 
neutrinos represent the first true signal and actually give some hints on possible 



extensions of the Standard Model. Baryon and lepton numbers are expected to 
be violated in Grand Unified Theories, leading to the fundamental consequences of 
proton instability and possibility of Majorana mass terms. The previously quoted 
mass range points in favour of lepton number violation at a scale ~ lO^ 14 ^ 15 ) GeV. 

Thus both gauge coupling unification and neutrino masses independently strongly 
suggest that beyond the electroweak scale there is another crucial scale where some 
new physics really arises and which lies just two orders of magnitude below the 
Planck scale. 

Neutrinos are curious also because, in contrast with the observed smallness of 
left handed quark mixings, it is now experimentally well established by Super- 
Kamiokande that they possess a large left handed mixing associated with the neu- 
trino flavours involved in atmospheric oscillations. To explain from the theoretical 
point of view the possible causes which determine the atmospheric mixing to be large 
and the possible consequences is the thread of the present review. 

Grand Unification represents an excellent framework for the study of neutrino 
masses and mixings, because it predicts fermion mass matrices to be related. Several 
Grand Unified models have been constructed which show that the presence of the 
large atmospheric mixing may naturally arise via the see-saw. These models can be 
distinguished into two classes, according to the basic mechanism exploited. 

Most models are based on the possible presence of a large right handed down 
quark mixing, not directly observable in weak transitions, which is strictly correlated 
to a large left handed mixing for charged leptons. The latter can be transferred to 
the Dirac neutrino mass matrix and the desired large mixing is achieved if the right 
handed Majorana neutrinos are not widely split. 

But there is also a particularly interesting alternative mechanism. This mech- 
anism is based on the observation that for a sufficiently hierarchical right handed 
neutrino spectrum, a large atmospheric mixing automatically arises if the left handed 
v^ and v T couple with comparable strength to the right handed v c . Thus, remark- 
ably, this situation does not require the presence of large mixings in the Dirac and/or 
Majorana mass matrices. As opposite to the first mechanism, this second one can 
then be realized even if all left and right mixings of quarks and leptons are very 
small. It is then a possibility for the large mixing to actually arise from nothing, 
that is only as an effect of the see-saw mechanism. This situation is particularly 
attractive from the theoretical point of view because it is compatible with left right 
symmetric scenarios. 

The first two chapters of this review contain a summary of the present experimen- 
tal status and future perspectives about neutrinos, show the two basic mechanisms 
providing an explanation of the large atmospheric mixing and present explicit Grand 
Unified models which realize these two possibilities. 



The two independent supports in favour of the physical relevance of M GUT , i.e. 
supersymmetric coupling unification and neutrino masses, can be hardly interpreted 
as just a coincidence. Armed with the results obtained from the study of the neutrino 
sector, we face the more ambitious task of constructing a sufficiently complete and 
realistic Grand Unified model which could embed one of the two mechanisms for the 
generation of the large atmospheric mixing. After 26 years, one could wonder if it 
is still possible and eventually what kind of informations could be added to further 
increase our knowledge of Grand Unification. There are mainly two new facts: 

1. much stronger limits on proton decay rates, in particular on the channels which 
are expected to be dominant in supersymmetric Grand Unified theories. Super- 
Kamiokande gives a lower limit for proton lifetime of about 2 • 10 33 yr for the 
channel p —* K + u); 

2. the approach itself. Since Grand Unification seems so successful, it is desirable 
to use it as a framework where to address also the flavour problem, that is to 
explain the observed hierarchical spectra and mixings of all fermions, including 
neutrinos. 

The flavour problem is independent of Grand Unification, in the sense that Grand 
Unification alone is not able to provide any light on this. However it is really an 
optimal framework in which to embed possible solutions of the flavour problem, 
because Grand Unification predicts more or less strict group theoretical relations 
between fermion mass matrices. One interesting way to address the flavour problem 
is by introducing an horizontal flavour symmetry, namely a symmetry under which 
the three copies of each fermion transform differently. 

Thus the modern point of view is to construct not too complicated, sufficiently 
realistic and relatively complete models. These models have to be enough detailed 
to actually carry out the calculation of as much low energy observables as possible. 
We are not satisfied with too rough and incomplete models. It is clear that only 
those models which provide predictions in agreement with experiment survive. 

In this respect, minimal versions of Grand Unified Theories (e.g. SU(5) with 
only 24, 5 and 5 of Higgs fields) are not realistic. They suffer from serious problems: 

• doublet-triplet splitting; 

• predictions of proton decay rates; 

• masses and mixings. 

To address some of these points, an enormous fine-tuning is required, e.g. of about 
14 orders of magnitude for the doublet-triplet splitting. But the serious fact is 



that, even if one tolerates to do this unpleasant tunings, minimal models give wrong 
predictions. In fact, minimal models predict protons decaying at a too fast rate 
for the channel p — > K + u, in conflict with experiment (point 1.) and give wrong 
relations between fermion masses (point 2.). 

It is clear that something is lacking in minimal models. In particular it would 
be of crucial importance to extend the Higgs sector in order to: 

• give an explanation to the doublet-triplet splitting; 

• enhance the prediction of proton lifetime; 

• adjust down quarks and charged lepton spectra and correctly describe neutrino 
masses. 

However, the addition of new fields in the Higgs sector must be done in such a 
way not to destroy the only firm starting point, that is coupling unification. Here 
troubles come because a realistic model requires so many improvements over the 
minimal versions that one could think that the beautiful idea of Grand Unification 
is actually not realizable. Note that the effort does not go in the direction of making 
precise predictions on poorly measured observables, but rather of doing predictions 
which are in not rough conflict with experiment. 

The problems and the technical difficulties which are typical of minimal models 
are reviewed in chapter 3, together with a review of the main ideas proposed to 
overcome them. 

Many attempts of constructing realistic models have been done in the context 
of 579(10) but in general they lack of predictivity on unification, which should then 
occur in an accidental way. This is due to the fact that despite SO (10) is interesting 
because it contains right handed neutrinos, it requires a lot of assumptions on the 
Higgs content in order to reduce its rank, thus reaching SU(5) or 577(4) <8> SU(2) ® 
SU(2). 

We instead focus the attention directly on Supersymmetric S77(5), which already 
possesses the same rank as the Standard Model, trying here to address as many 
problems as possible. We present a model which successfully overcomes fine-tunings 
problems and which gives predictions in agreement with experiments. This model 
exploits the crucial presence of an additional U(l) flavour symmetry. 

The fine-tunings linked to the doublet-triplet splitting are avoided by adopting 
the Missing Doublet Mechanism and by using the flavour symmetry to stabilize 
the Higgs doublet mass against large mass corrections due to non renormalizable 
operators. Fermion masses and mixings are reproduced, at the level of the correct 
order of magnitude, by the same flavour symmetry. In the neutrino sector the 
preferred solution is one with nearly maximal mixing both for atmospheric and 



solar neutrinos, which turn out to be compatible with the LOW solution. The 
large atmospheric mixing follows from the presence of large left mixings both in 
the Dirac neutrino and charged lepton mass matrices. The latter is correlated to 
an unobservable large right mixing of the down quark mass matrix. A third effect 
of the flavour symmetry is to adequately suppress non renormalizable operators 
inducing proton decay at a too fast rate. Remarkably the presence of the large 
representations demanded by the Missing Doublet Mechanism for the solution of the 
doublet-triplet splitting problem, makes possible a decrease of the predicted value 
of a^(mz) assuming coupling unification thanks to an increase of the effective mass 
that mediates proton decay induced by dimension 5 operators. As a consequence 
the value of the strong coupling is in better agreement with the experimental value 
and the proton decay rate is still compatible with present bounds, even if it is 
difficult to avoid the conclusion that proton decay must occur with a rate at reach of 
current and next generation experiments. Chapter 4 contains a detailed description 
of this model. We think that this model is interesting because it proves that a 
Supersymmetric SU(5) Grand Unification is not excluded and offers a benchmark 
for comparison with experiment. 



Chapter 1 

Neutrino Masses and Mixings 



Recent data from Super-Kamiokande |l] have provided a more solid experimental 
basis for neutrino oscillations as an explanation of the atmospheric neutrino anomaly. 
Also the solar neutrino deficit, observed by several experiments, is probably an 
indication of a different sort of neutrino oscillations. Neutrino oscillations imply 
neutrino masses, so that it is compelling to look for viable extensions of the Standard 
Model. In order to fulfill such a program, it is important to recognize that the 
extreme smallness of neutrino masses in comparison with quark and charged lepton 
masses indicates that the formers have different nature, linked to lepton number 
violation. 

Thus, neutrino masses provide some insight on the very large energy scale where 
lepton number is violated and on Grand Unified Theories (GUTs). Indeed, GUTs 
constitute a very natural framework for the generation of neutrino masses, because 
they generically predict lepton and baryon number violation. Neutrino masses could 
also give an important feedback on the problem of quark and charged lepton masses, 
as all these masses are related in GUTs. In particular the observation of a nearly 
maximal mixing angle for v^ — > v T is particularly interesting. At present, a large 
solar neutrino mixings seems to be slightly more preferable than a small one. Large 
mixings in the neutrino sector are very interesting because a first guess was in 
favour of small mixings, in analogy to what is observed for left mixings in the quark 
sector. If confirmed, single or double maximal mixings can provide an important 
hint on the mechanism that generate neutrino masses. With neutrino masses settled, 
observation of proton decay will be the next decisive challenge remained to support 
or eventually put in crisis GUTs. 

After a concise summary of the present experimental status of neutrino masses 
and mixings and its future perspective, in this chapter we will review the most 
promising mechanisms which reproduce a large atmospheric mixing. 



8 CHAPTER 1. NEUTRINO MASSES AND MIXINGS 

1.1 Neutrino Oscillations: Present Experimental 
Status and Future Perspective 

Despite the remarkable progress reached during the last few years, the experimental 
status of neutrino oscillations is still very preliminary. While the evidence for the 
existence of neutrino oscillations from solar and atmospheric neutrino data is rather 
convincing by now, the values of the mass squared differences |Am 2 | and mixing 
angles are not firmly established. However, there is a general confidence that future 
experiments will not only succeed in clarifying such a picture but hopefully will also 
find one of the CP violating phases of the neutrino sector and establish the sign of 
the Am 2 . 

1.1.1 Atmospheric 

Present 

Data with impressive high statistics are coming from Super-Kamiokande |1[] which 
are in perfect agreement with the hypothesis of flavour oscillations || Q. Such 
hypothesis is consistent with all Super-Kamiokande data and is also corroborated 
by independent atmospheric neutrino results from the Soudan-2 Q and MACRO || 
experiments, as well as by the pioneering Kamiokande experiment ||. 

Zenith-angle distributions, up-down asymmetries and so on provide plenty of 
smoking guns for the fact that something actually happens to atmospheric v^ as a 
function of LE n with n « — 1, where L is the baseline length and E the energy. 

In the simplest 2z/ case of v^ — > u T , a Super-Kamiokande combined fit indicates 
that: 

(i) the preferred value of Am 2 tm lies in the range (1.5 -j- 5) • 10~ 3 eV 2 ; 
(ii) a large mixing angle appears to be well established: sin 2 20 atm > 0.88 at 90% 
confidence level. 

Interpreted in terms of 2v oscillations, Super-Kamiokande favours v^ — ► v T over 
Vfj, —* v s at 99% confidence level, based on separate analyses of the zenith-angle 
distributions. The same conclusions also come from the rate of tt° production, but 
in this case systematic uncertainties are still large. Pure v^ —>■ v e transitions do not 
provide a good fit to the data, and are independently excluded by the negative v e 
disappearance searches in the CHOOZ ]7| and Palo Verde [§] reactors. 

However, these results have to be taken cum grano salts. Additional oscillation 
channels may be open, as naturally expected in 3u and 4z/ schemes accommodating 
the current phenomenology. In fact, it has been realized that the dominant v^ — > 
v T oscillations plus subdominant v^ -^ v e oscillations are also consistent with 
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combined Super-Kamiokande and CHOOZ data, leading to a much richer three- 
flavour phenomenology, as will be discussed in the following. On the other hand, 



when the data are analyzed [H], |TT| for v^ oscillating into a linear combination of 
v T and v s , one finds such a scenario to be still viable (and also compatible with 
solar neutrino data) for a sizeable region of parameter space. These analysis are 
very important to determine to what extent the 4z/ scenario motivated by the LSND 
experiment can be accommodated. 

It is clear that the definite evidence would be direct r production, but this is 
not an easy task for Super-Kamiokande. Actually, there is not yet direct evidence 
for an oscillation pattern. 

Thus, there is still room to speculate on other interpretations, like neutrino decay 
T2|| , decoherence scenarios |13| and v^ mixing with neutrino states propagating in 



large extra dimensions ||14| , which appear indistinguishable with present data with 
respect to the L/E distribution. Detailed analysis together with high statistics data 
could test these non standard explanations. For instance, the exotic scenario of 
neutrino decay seems to fail ||, [L5| in reproducing the zenith-angle distributions of 
sub-GeV, multi-GeV and up-going muons observed by Super-Kamiokande. 

Future 

So, future experiments are highly needed, among which long-baseline laboratory 

experiments, which require a muon decay source for the neutrino beam, seem to be 

very promising in order to: 

(i) confirm the existence of atmospheric neutrino oscillations by using a controllable 

beam and measuring an oscillation pattern; 

(ii) precisely measure |Am^ 4m | and sin 2 29 atm ; 

(iii) make new measurements possible. 

The K2K experiment |§ (L=250 km, < E v >= 1.4 GeV, from KEK to Super- 
Kamiokande, in operation), finds an encouraging deficit of events with respect to the 
non oscillation hypothesis. The rate is perfectly compatible with Am 2 tm = 3 • 10 -3 
eV 2 and sin 2 26 atm = 1, but the statistics is still low and the significance only about 
2a. 

Long-baseline investigations will be pursued by MINOS |I7| (L=732 km, from 
Fermilab to Soudan, begin in 2002). In addition to (i) and (ii), it should be able 
to discriminate between v^ — > v T and v^ — > v s and to measure | t/ e 3 1 from v e — > v^ 
appearance. 

OPERA and ICANOE @ (L=743 km, from CERN to Gran Sasso, approved), 
have the main purpose of looking for r detection as the smoking gun for v^ — > v T 
oscillations. They also should be sensitive to \Ug\. 



10 
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With respect to long-term future, neutrino factories together with very long 
baselines (L ~ 3000 km) appear [E| to offer the possibility of determining the sign 



of Am 2 atm from matter effects on the Earth crust, measuring with an accuracy of 
tenths of degrees the angle #13, down to values of #13 = 1°. The CP violating phase 
5 could also be measured if the solar neutrino parameters lie in the range of the 
MSW-LMA solution. 

1.1.2 Solar 

Present 

For solar neutrinos the experimental situation and the interpretation in terms of 
frequency and mixing angles is less clear. 

The solar neutrino experiments |2(J sample different v e energy ranges and find 
different flux deficits compared to the Standard Solar Model (SSM) [[H|], as shown 
in Table |P| , so that the v e survival probability is inferred to be energy dependent. 





flux( cx P) / flux( SSM ) 


GALLEX 
SAGE 

Homestake 

Super-Kamiokande 


0.60 ±0.06 
0.52 ±0.06 

0.33 ±0.03 

0.47 ±0.02 



Table 1.1: For the energy range specific of each type of experiment, the ratio between 
the experimentally measured neutrino flux and the flux predicted by the SSM is 
shown. 



Global oscillation fits have been made which include: total rates (the different 
suppression ratio for Homestake plays a vital role), day-night asymmetry, seasonal 



dependence beyond 1/r 2 . Typical candidate solar solutions are given in Table [L2 



Interesting recent developments by Super-Kamiokande (with reduced background 
and threshold of 5 MeV) show that the flux deficit is confirmed but no other direct 
signal of oscillations emerges: the spectrum is compatible with flat, the day-night 
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solution 


AmUeV 2 ) 


sin 2 28 sol 


SMA-MSW 


-5-10" 6 


~5- 1(T 3 


LMA-MSW 


~2- 1(T 5 


~ 1 


LOW 


~ 1(T 7 


~ 1 


VO 


~ iq-( 9 - 10 ) 


~ 1 



Table 1.2: Typical solar solutions from global fits. 



asymmetry is reduced at the 1.3 a level and the seasonal variation is compatible 
with the flux variation from the orbit eccentricity. These last two facts disfavour 
respectively the SMA-MSW and the VO at the 95% confidence level. We are left 
with the LMA-MSW and the LOW, both with large mixing. The pure transition 
u e — ► u s is disfavoured at 95%. It is important to keep in mind that the Am^ 
values of the above solutions are determined by the experimental result that the 
flux suppression is energy dependent. If the Homestake indication were disregarded, 
new energy independent solutions would emerge, with much more extended ranges 
for Am 2 sol . 

All these conclusions need to be confirmed by updated global analysis, without 
forgetting the region with 9 so i > it/A, as well as matter and vacuum effects. It is 
also important to analyze the data in the 3u and 4u schemes, to understand to what 
extent a component u e 



v. 



is tolerated [10|. 



Future 



The experimental situation will be hopefully clarified with the results from future 



experiments like SNO [||, BOREXINO [fH and KamLAND |4 



In the present Phase I, SNO clearly sees the signal from charged current (CC) 
reactions u e +d — > p+p+e~ , with an energy spectrum which seems to be proportional 
to that of the SSM. It also sees a clear signal from elastic scattering (ES) J2i v % + e ~^> 
J2i Vi + e. The ratio CC/ES is expected to be measured in a sufficiently accurate way 
to reveal at 3a if high energy neutrinos with i ^ e are coming from the Sun at a rate 
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consistent with the observed deficit of u e with respect to the SSM expectation. In the 
following Phase II, the neutral current (NC) cross section £\ Vi + d^>p + p + ^\ z/j 
will be measured. With also the ratio of NC/CC, SNO should be able to discriminate 
between different solar neutrino scenarios. 

BOREXINO can measure in real time the "' Be flux and is sensitive to the VO 
solution. 

KamLAND, starting in 2001, should be able to test definitively the LMA-MSW 
and LOW solutions, using neutrinos from nuclear power reactors. It should also 
measure the day-night and seasonal effects for 7 Be neutrinos. 

Thus it seems possible that the solar neutrino problem will be better clarified in 
the near future. 



1.1.3 LSND 



The issue of the additional claim for v e — > v^ neutrino oscillations from the LSND 



experiment [[25j is still open, because it is neither confirmed nor completely excluded 
by KARMEN-2 [2B|. The range still surviving is at small mixing angle and with 
Am 2 = (0.2-7-2) eV 2 . This range is much above the solar and atmospheric neutrino 
frequencies, so that the need to accommodate a third frequency would imply the 
existence of a fourth light sterile neutrino, in addition to the three established weak 
interacting flavours, which should mix with solar and/or atmospheric neutrinos. 
MiniBOONE, starting in 2002, is an experiment designed to solve this problem. It 
will be possibly followed by BOONE, if the signal is confirmed. 

1.1.4 Three Flavour Analysis 

Given the present experimental uncertainties the theorist has to make some assump- 
tions on how the data will finally look like in the future. Here and in all the following 
sections, we tentatively assume that the LSND evidence will disappear. If so then 
we only have two oscillations frequencies, which can be given in terms of the three 
known species of light neutrinos without additional sterile kinds. 

As pointed out, Super-Kamiokande data are compatible with nearly maximal 
v^ — »■ v T two flavour oscillations, but, among the possible deviations from this stan- 
dard picture, three flavour oscillations appear to be perfectly viable. In fact, Super- 
Kamiokande data do not exclude a non- vanishing U e ^ element. 

We then take for granted that the frequency of atmospheric neutrino oscillations 
will remain well separated from the solar neutrino frequency, even for the MSW 
solutions: /S.m 2 atm > Am^. Assuming that two out of three active neutrinos are 
almost degenerate, say m\ ~ m 2 , one can show || that atmospheric neutrinos probe 
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only Am„ tm = m 3 — m{ 2 an d the mixing matrix elements U a %, which satisfy the 
unitarity constraint U 2 3 +Ufa + U 2 3 = 1. The unitarity constraint can be conveniently 
embedded in a triangle plot. The best fit for Super-Kamiokande data (70.5 kTy) is 
found at (Am 2 atm , U 2 3 , Ufa, U 2 3 ) = (3.5 ■ W~ 3 eV 2 , 0.07, 0.57, 0.36). The tolerated 



mixing is constrained by U 2 3 < 0.31 at 90% confidence level. 



V^ — > V e iiiiAiiig, io Luiinnaincu uy u e3 

The weak preference for U 2 3 ^ is suppressed by including CHOOZ data. 
The best fit for combined Super-Kamiokande and CHOOZ data gives (Am 2 vtm , U 2 3 , 
Ufa, U 2 3 ) =(3.0 ■ 10~ 3 ey 2 , 0., 0.5, 0.5), with the further constraint U 2 3 < 0.04 at 90% 
confidence level. 

It appears very difficult to probe, with present atmospheric data, values of U 2 3 
as small as a few %. Constraining U 2 3 is one of the major tasks for future reactor 
and accelerator neutrino experiments. 



1.2 Beyond the Standard Model 

Assuming just three flavours, neutrino oscillations are due to a misalignment between 
the flavour basis v' = [y e , v^, v r ), where v e is the partner of the mass and flavour 
eigenstate e~ in a left-handed weak isospin SU(2) doublet (similarly for v^ and v T ), 
and the neutrino mass eigenstates v = (i>\, v 2 , v^): 

\u') = U\u), (1.1) 

where U is the MNS 3 by 3 mixing matrix j3|. The presence of mixing implies that 
neutrinos cannot be all massless and actually the presence of two different oscillation 
frequencies implies at least two different non zero masses. Neutrino oscillations are 
practically only sensitive to differences Am 2 so that the absolute scale of squared 
masses is not fixed by the observed frequencies. 



1.2.1 Neutrino Masses and Lepton Number Violation 

Neutrino oscillations imply neutrino masses which in turn demand either the ex- 
istence of right handed neutrinos (Dirac masses) or lepton number (L) violation 
(Majorana masses) or both. Anyway, one is forced to go beyond the Standard 
Model (SM). In fact, referring to three neutrino generations, in the SM Lagrangian 
density only the term H7^ T pu L appears, where D^ is the SU{2) <8> U(l) gauge co- 
variant derivative and vl = [y e v^ z/ r ) T is a column of three Weyl spinors in 
a four- component notation (a Weyl spinor has only two independent components) 
satisfying vj_, = 1/2 (1 + 75)^. The standard convention is to assign L(vl) = 1. 



14 CHAPTER 1. NEUTRINO MASSES AND MIXINGS 



If no additional right handed neutrino field Nr = (N e N^ N T ) T (satisfying 
Nr = 1/2 (1 — ^Nr and with L(Nr) = 1) is allowed, there isn't the possibility of 
including in the Lagrangian the Dirac mass term 

C-Dir 3 -VZ T m D N R + h.c. (1.2) 

which, by defining v = Vl + Nr, can be written as Coir 3 —V T m£,v, mo being 
a matrix in flavour space. The above term do not violate L, which is an additive 
quantum number. 

For a massive charged particle of spin 1/2 one needs four states, while only two 
are enough for an intrinsically neutral particle. If L is violated, there is no conserved 
quantum number that really makes neutrinos and antineutrinos different and a new 
type of mass term is possible. We can have a Majorana mass term 

1 T l T 

£-Maj 3 —^iy c )R m L LVL + h.c. = --i/ L Cm LL v L + h.c. 

= --vlm LL v L + h.c. (1.3) 

where {v c )r = CUZ 1 and C is the 4 by 4 matrix in Dirac space that implements charge 
conjugation. We denote by t the transposition in Dirac space, while maintaining 
T for the transposition in flavour space. The Majorana mass term violates L by 
two units. Also, since in the SM vl is a weak isospin doublet, it transforms as a 
component of an isospin triplet. In the following, as we are only interested in flavour 
indices and not in Dirac indices, we will simply denote the Majorana term for left 
handed fields by vJ j 'mLL 1/ L J rh.c., omitting the Dirac matrix C and the transposition 
t. Note that if L is violated and Nr also exists, then a second type of Majorana 
mass is also possible which is 

C M aj 3 -^Nlm RR N R + h.C. , (1.4) 

where we again omitted C and t. Clearly it also violates L by two units, but, since 
N R is a gauge singlet, this term is invariant under the SM gauge group. 

In conclusion, if Nr does not exist, we can only have a Majorana mass ttill if L 
is violated. If iv"^ exists and L is violated, we can have both Dirac tud and Majorana 
masses m^ and m RR . 

If one wants to give masses to neutrinos but to avoid the conclusion that L is 
violated, then one must assume that iv"^ exists and that neutrinos acquire Dirac 
masses through the usual Higgs mechanism as quark and leptons do. Technically 
this is possible. But there are two arguments against this possibility. 
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1) The first argument is that neutrino masses are extremely small so that the cor- 
responding Yukawa couplings would be enormously smaller than those of any other 
fermion. Note that within each generation the spread of masses is by no more than 
a factor 10 3 ^ 5 . But the spread between the t quark and the heaviest neutrino would 
exceed a factor of 10 11 . 

2) A second argument arises from the fact that once we introduce Nr in the theory, 
then the L violating term N^turrNr is allowed in the Lagrangian density by the 
gauge symmetry. In the minimal SM, i.e. without Nr, L and B conservation are ac- 
cidental global symmetries that hold because there is no operator term of dimension 
< 4 that violates B and L but respects the gauge symmetry. On the other hand, in 
the presence of N R , the dimension 3 operator N^turrNr is gauge symmetric but vi- 
olates L. We expect a mass factor in front of this operator in the Lagrangian density, 
and in the absence of a protective symmetry, we expect it to be of the order of the 
cut off, i.e. of order Mqut or larger. Thus, L number violation is naturally induced 
by the presence of Nr, unless we enforce it by hand. Once we accept L violation we 
gain an elegant explanation for the smallness of neutrino masses as they turn out 
to be inversely proportional to the large scale where lepton number is violated. 

1.2.2 Smallness of Neutrino Masses 

Given that neutrino masses are extremely small, it is really difficult from the theory 
point of view to avoid the conclusion that L must be violated. In fact, it is only 
in terms of lepton number violation that the smallness of neutrino masses can be 
explained as inversely proportional to the very large scale where L is violated, of 
order Mqut- 

Assuming three neutrino flavours i/ L = [y e v^ u T ) T , consider the following two sit- 
uations (in the neutrino interaction basis and after having diagonalized the charged 
lepton mass matrix). 

• If L is not conserved, even in the absence of Nr, Majorana masses can be gen- 
erated for neutrinos by allowing in the Lagrangian the presence of dimension 
five operators of the form 

5 2 M 2 M W ,W ' K J 

with = (0 + 0°) T being the ordinary Higgs doublet, L = [y' L li) T the weak 
lepton doublet, A a matrix in flavour space and M a large scale (• is put to 
remember the weak isospin). Neutrino masses generated by O5 are then of 
the order m u ~ v 2 /M for \j ~ 0(1), where (0°) = v ~ 0(100 GeV) is the 
vacuum expectation value of the ordinary Higgs. 
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The existence of a right handed neutrinos is quite plausible because all GUT 
groups larger than SU(5) require them. In particular, the fact that a field 
with the same quantum numbers of a right handed neutrino completes the 
representation 16 of 5*0(10): 16 = 5 + 10 + 1, so that all fermions of each 
family are contained in a single representation of the unifying group, is too 
impressive not to be significant. At least as a classification group SO(10) must 
be of some relevance. Assuming that there are both N R = (N e N^ N T ) T and 



L violation, the see-saw mechanism [28] is possible. Consider the mass terms 



in the Lagrangian corresponding to Dirac and RR Majorana mass matrices 
(we consider LL Majorana mass terms as comparatively negligible) in the 
interaction basis in which charged leptons are diagonal: 

Cseemw = -V~Jm D N R + -N^M RR N R + h.c. (1.6) 

The 3 by 3 matrices mo and M RR are the Dirac and Majorana mass matrices 
in flavour space (M RR is symmetric, M RR = M RR , while mo is, in general, 
non hermitian and non symmetric). One expects mo to come from the Higgs 
mechanism: Cy u k 3 — v' l yDN R {4>°). We expect the eigenvalues of M RR to be 
of order M GUT or more because RR Majorana masses are SU(3) x SU{2) x U(l) 
invariant, hence unprotected and naturally of the order of the cut off of the 
low-energy theory. Since all N R are very heavy we can integrate them away 
and the resulting effective neutrino mass matrix is: 

L ~ = -\^K f ^L + h.c. (1.7) 

with 

m e J f = m D M RR m T D . (1.8) 

This is the well known see-saw mechanism ]^8| result: the light neutrino masses 
are quadratic in the Dirac masses and inversely proportional to the large Ma- 
jorana mass. If some N R are massless or light they would not be integrated 
away but simply added to the light neutrinos. 

Here we assumed that the additional non renormalizable terms from 0% are com- 
paratively negligible, otherwise they should simply be added. After elimination of 
the heavy right-handed fields, at the level of the effective low energy theory, the 
two types of terms are equivalent. In particular they have identical transformation 
properties under a chiral change of basis in flavour space. The difference is, how- 
ever, that in the see-saw mechanism, the Dirac matrix mo is presumably related to 
ordinary fermion masses because they are both generated by the Higgs mechanism 
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(rriD = Udv) and both must obey GUT-induced constraints. Thus, if we assume 
the see-saw mechanism more constraints are implied. In particular we are led to 
the natural hypothesis that m£> has a largely dominant third family eigenvalue in 
analogy to mt, mb and m T which are by far the largest masses among u quarks, d 
quarks and charged leptons. Once we accept that rri£, is hierarchical it is very dif- 
ficult to imagine that the effective light neutrino matrix, generated by the see-saw 
mechanism, could have eigenvalues very close in absolute value. 

Then, assuming that one neutrino mass, presumably the third generation neu- 
trino mass m 3 , is much larger than any other, we have Am 2 atm « m§, which means 
\m^\ ~ 0.05 eV. If we apply the see-saw formula assuming that the Majorana RR 
masses lie around M, then, with \m^\ ~ 0.05 eV, 7Bd 33 ~ 200 GeV (like the top 
quark mass or the Higgs vev), we find M ~ 0.8- 10 15 GeV, a value amazingly close to 
Mqut- This supports the initial assumption that indeed neutrino masses are related 
to the scale of L violation. 

1.2.3 Approximate Mixing Matrix in the Neutrino Sector 

Given the definition of the mixing matrix U in eq. (|1.1|) and the transformation 
properties of the effective light neutrino mass matrix ml 



i^JJ fTi J J 



*/?mf/V r 



p T L OJ T m e JfU)u L 



U T m e J f U = Dmg[e l ^m 1 ,e i42 m 2 ,m 3 ] = m diag (1.9) 

where mi : 2,3 are real and positive, we obtain the general form of rrfj*: 

m e Jf = Um dmg U T . (1.10) 

The MNS matrix U can be parameterized in terms of three mixing angles and one 
phase, exactly as for the quark mixing matrix Vqkm- 



U 



C12C13 S12C13 S13 

-snC23e lS - C12S13S23 ci2C 23 e 4<5 - S12S13S23 C13S23 

Sl2S23e lS - C12S13C23 -Ci 2 S23e tS - S12S13C23 C13C23 



:i.n} 



In addition we have the two phases <pi and 02 that are present because of the 
Majorana nature of neutrinos. Thus, in general, 9 parameters are added to the SM 
when non vanishing neutrino masses are included: 3 eigenvalues, 3 mixing angles 



and 3 CP violating phases p9 



Maximal atmospheric neutrino mixing and the requirement that the electron 
neutrino does not participate in the atmospheric oscillations, as indicated by the 
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Super-Kamiokande and CHOOZ data, lead directly to consider as good approxima- 
tion: \Ufja\ = \U t3 \ = 1/\/2 andU e3 = 0. Hence, c 2 3 = s 2 3 = l/v^2andsi 3 = 0, Ci 3 = 1 
(but different sign conventions are possible plf ) and the following structure of the 
Ufi (/ = e, /i, r, i = 1, 2, 3) mixing matrix is obtained: 



U 



fi 



Cl2 _ Si 2 

-S12/V2 C12/V2 1/V2 
s 12 /V2 -C12/V2 +1/^2 



;i.i2) 



This result refers to the case of arbitrary solar mixing angle S12 = sin# so /, so that 
C12 = S12 — I/a/2 holds for maximal solar mixing. In the limit t/ e 3 = 0, all CP 
violating effects vanish and we can neglect the additional phase parameter 5 generally 
present in [Iff. the matrix U is real and orthogonal and equal to the product of a 
rotation by 7r/4 in the 23 plane times a rotation in the 12 plane. 

While the simple parameterization of the matrix U in eq. ( 1.12 ) is quite useful 
to guide the search for a realistic pattern of neutrino mass matrices, it should not 
be taken too literally. In particular, as stressed in section |1.1.4j , the data do not 
exclude a non- vanishing U e3 element. 



1.3 How to Reproduce the Large Atmospheric 
Mixing 

1.3.1 Mass Hierarchies and Approximate Zeroth Order Tex- 
tures 



Note that since we are assuming only two frequencies, given by 



A so ; = m 2 -m x , A atm 

there are three possible hierarchies of mass eigenvalues 



m l ~ m l,2 » 



HI: 



;i.i3) 



A 
B 
C 



\m 3 \ » |ra 2 ,i| 

\mi\ ~ \m 2 \ » \m 3 \ 

\mi\ ~ |m 2 | ~ \m 3 \ 



;i.i4) 



(in case A there is no prejudice on the mi, m 2 relation). For B and C different 
subcases are then generated according to the relative sign assignments for m lj2 ,3. 
Using eqs. Q1.9Q and ( |1.12j ), for each case, by setting to zero the small masses, it 



is possible to find |32] , |33[| the effective light neutrino matrices which automatically 
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lead to maximal atmospheric mixing []. Such zeroth order texture can be taken as 
guides in the construction of specific models for neutrino masses and mixings. 
For example, in case A, m diag = Diag[0, 0,m 3 ] and we obtain 



A : m di , 



nig 



Diag[0, 0, 1] m 3 



nr 



eff 









0] 


n 


1 


1 




■' 


■> 





1 

2 


1 

2 - 



m 3 



(1.15) 



In this particular case the results are the same for double and single maximal mixing. 
Note that the signs correspond to the phase convention adopted in eq. ( |1.12|) . 

One example of case B follows by taking m diag = Diag[m 1; —mi, 0] with double 
maximal mixing, so that: 



B : m diag = Diagfl, -1, 0] m x 



77) 



eff 



1 1 





o -4* -4 
i 

1 



-4 o 







mi 



;i.i6) 



A strong constraint arises in case C from the non observation of neutrino-less 
double beta decay which requires that the ee entry of nfj' must obey |(m^) ee | < 
As observed in ref. 



0.2 eV B34 



35 



if the average neutrino mass is around 1 eV, it 
can only be satisfied if bimixing is realized (that is double maximal mixing, with solar 
neutrinos explained by the VO solution or maybe by the large angle MSW solution). 



The reason is that, as results from eqs. ( P-.9| ) and ( |1.12j ) 



(m 



\eff 

v)ee 



mic\ 2 + m 2 sl 2 . 



For the required cancellation one needs opposite signs for mi and m 2 and comparable 
values of c? and s? , that is nearly maximal solar mixing. In this case, one example 



'12 

of case C is: 



'12) 



C : m diag = Diag[-1, 1, 1] m 1 



m 



eff 
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1 




V2 


V2 
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/?, 
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V2 
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2 



mi 



;i.i7) 



It is possible to proceed similarly in all the other cases |33|. Once a solution to 
the solar neutrino problem is chosen and a set of suitable small perturbation terms 
is introduced, oscillation phenomena are unable to distinguish between the cases A, 
B and C . However, from the model building point of view, each texture represents 
an independent possibility in a zeroth order approximation. 



^^Note that here we are working in the basis where the charged lepton masses are diagonal, 
and approximately given by mi — Diag[0, 0, m T \. For model building one has to arrange both the 
charged lepton and the neutrino mass matrices so that the neutrino results coincide with those 
given here after diagonalization of charged leptons. 
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1.3.2 Degenerate Neutrinos 

Configurations B and C imply a very precise near degeneracy of squared masses. 
For instance, case C with rrii ~ 1 eV would require a relative splitting \Am/m\ ~ 
Am 2 , tm /2m 2 ~ lCT^ 4 ) and a much smaller one for solar neutrinos, especially if 
explained by the VO solution: \Am/m\ ~ 10~( 10 ^ n \ Case C is the only one in 
which neutrinos could represent a significant component of hot dark matter. An 
intermediate possibility between the case of degenerate neutrino masses and the 
case of hierarchical ones is the situation B in eq. (|1.14j) . With two almost degen- 



erate heavier states and a nearly massless neutrino we cannot reach a mass range 
interesting for cosmological purposes. However, since the experimentally accessible 
quantities are the squared mass differences, case B remains an open possibility. 

As mentioned above, it appears to be unplausible that starting from hierarchical 
Dirac matrices we end up via the see-saw mechanism into a nearly perfect degen- 
eracy of squared masses. Thus models with degenerate neutrinos could only be 
natural if the dominant contributions directly arise from non renormalizable opera- 
tors like O5 in eq. ( |1.5|) because they are a priori unrelated to other fermion mass 
terms. The degeneracy of neutrinos should be guaranteed by some slightly broken 
symmetry, which could also give the small terms required to go beyond the zeroth 
order approximation. Models based on discrete or continuous symmetries have been 



proposed |36 



Notice however that, for the VO scenario, even if arranged at the GUT scale, it 
is doubtful that such a precise degeneracy could be stable against renormalization 
group corrections when running down at low energy unless it is protected by a 
suitable symmetry |37|, |38| . For instance, working only with the light neutrinos, the 
U(1)q charge Q = [L e — L M — L T ) allows at leading order [[$9|] for the texture of 
the example in (|1.16| ). The non- vanishing entries (m u )i2 and (m u )i 3 are expected to 
be of the same order, although not necessarily equal. The vanishing entries can be 
filled by the ratio of the vev (ip) of a scalar field carrying two units of Q and a mass 
scale M providing the cut off to the low-energy theory: 



m„ = m 



e 1 1 
1 e e 
1 e e 



:i.i8) 



where e = (<p)/M and only the order-of- magnitudes are indicated. For values of e 
smaller than 1 one obtains a small perturbation of the zeroth order texture. The 
mixing matrix has a large, not necessarily maximal, mixing angle in the 23 sector, 
a nearly maximal mixing angle in the 12 sector and a mixing angle of order e in 
the 13 sector |40|. The mass parameter m should be close to 10~ ( - lH ~ 2 ) eV to provide 
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the frequency required by the atmospheric oscillations. Finally, m\ — m\ ~ m 2 e. In 
this model the VO solution to the solar neutrino deficit would require a very tiny 
breaking term, e ~ O(10 -7 ). The symmetry protects the degeneracy of m\ and m 2 



against renormalization group corrections |38], [4T |. 

This mechanism is a typical example of such models which are very attractive at 
first glance, but which, after a deeper analysis, turns out to be incomplete. In fact, 
the symmetry introduced is not able to account for the spectra of charged leptons 
- which, in this model, would be roughly degenerate - and of quarks. Thus, in 
order to construct a complete model of fermion masses, one has to do some further 
assumption. For instance one can introduce an ad hoc flavour symmetry which 
fulfills this task without damaging the effective neutrino mass matrix. 



1.3.3 Hierarchical Neutrinos 

We now discuss models of type A with large effective light neutrino mass split- 
tings and large mixings. Reconciling large splittings with large mixing would seem 
difficult. Indeed, one could guess that, in analogy to what is observed for quarks, 
large splittings correspond to small mixings because only close-by states are strongly 
mixed. In particular it is of basilar importance to find a mechanism which leads to 
the large atmospheric mixing. Actually, there are mainly two ways in which this 
situation can be achieved. 



Large Mixing from Stretching 

The first one arises by stretching some parameters |jEj| . For instance, one could 



start with not too small 23 mixings both in the charged lepton and in the effective 
neutrino mass matrices. For instance, take 9i ~ 6 V ~ 9c ~ 13° where 9c is the 
Cabibbo angle. If the relative phases are such that the two angles just add, then 
from two not too small angles one can end up with a large one. Clearly the situation 
is not very satisfactory because the largeness of the mixing is completely accidental 
rather than descending as a consequence of some mechanism. In addition, one can 
expect at most the mixing to be large, but not at all to be exactly maximal. 



Large Mixing from the See-Saw 



Via the see-saw mechanism |£8| , there are other two particularly simple subcases in 
which a large atmospheric mixing can be realized. They are based on the observation 
that, in a 2 by 2 matrix context, the requirement of large splitting and large mixings 
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leads to a condition of vanishing determinant. For example the matrix 



m oc 



2 
X X 



X 



;i.i9) 



has eigenvalues and 1 + x 2 and for x of 0(1) the mixing is large. Thus, in the 
limit of neglecting small mass terms of order 7711,2, the demands of large atmospheric 
neutrino mixing and dominance of 7773 translate into the condition that the 2 by 2 
sub determinant 23 of the 3 by 3 mixing matrix vanishes. The problem is to show that 
this approximate vanishing can be arranged in a natural way without unacceptable 
fine tuning. 

Without loss of generality, leaving aside for the moment the possible presence 
of flavour symmetries, one can go to the basis where both the charged lepton Dirac 
mass matrix mi and the Majorana matrix Mrr of the right handed neutrinos are 
diagonal. For simplicity, we start assuming that the role of the first generation is 
not crucial in the mechanism for the generation of neutrino masses, so than one can, 
with good approximation, work in the 2 by 2 case (in the next section this condition 
will be relaxed). Writing m^ and Mrr, defined as in eq. (|1.6|) , in the most general 
way 



m D 



M, 



RR 



M 2 





M 3 



1.20) 



where v is a vacuum expectation value, 
the see-saw, one obtains: 



a, b and c are Yukawa couplings, then, via 



m 



eff 



V 

Ms 



a" 

~M 2 
a b 
VM 2 



+ 



+ 



c 

Ms 
c 

Ms 



+ 



a b 

M~ 2 

M 2 + M 3 . 
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Ms 

1 



1.21^ 



The request of large splittings among the light neutrino eigenvalues is equivalent to 
demanding that the determinant of vrfj* is much smaller than its trace. It is then 
possible to see at glance that two cases arise exhibiting large mixing, that is when 
the terms with respectively M 3 or M 2 at the denominator are dominant. 



Case I: Large mixing from mp 

One simple example of the first case is realized if M 2 



M 3 and a, b < 1. 



In order to have a large splitting, one must have c ~ 1, that is the right 



handed neutrino of the third generation couples with the same strength |43 
to left handed v^ and v T . The heaviest mass for light neutrinos is then 777,3 ~ 
v 2 /Ms. As already pointed out, in the hierarchical case the data from Super- 
Kamiokande give 7773 ~ 0.05 eV so that, if one assumes that v is a typical weak 
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scale, namely 200 GeV, then M 3 
scale. 



10 15 GeV, just a bit lower than the GUT 



It is worth to stress that this first mechanism makes use of a strongly asym- 
metric Dirac matrix, that is a matrix with a large left handed mixing already 
present. It has been observed 133, M K5L f55L g7| that in SU(5) left handed 



charged lepton mixings correspond to right handed mixings for down quarks 
(more on this later on). Large right handed quark mixings are not in contrast 
with experiment and viable GUT models that correctly reproduce the data on 
fermion masses and mixings can be constructed following this mechanism. 

Case II: Large mixing from nothing 

An alternative possibility j|8| is to have dominance of the terms with M2 at the 
denominator. This is achieved for any c < 1 if a 2 , b 2 > M 2 /M 3 . The request 
for large splitting is then equivalent to the requirement that also a ~ b. Now 
it is the second generation right handed neutrino which is particularly light 
and which couples with the same strength [43j to left handed z/ M and u T . In 
order to be more specific, consider the following example with symmetric ma- 
trices. These matrices are interesting because, for instance, one could want to 
preserve left right symmetry at the GUT scale. Then, the observed smallness 
of left handed quark mixings would also demand small right handed mixings. 
Starting from 



m D = v 



e xe 
xe 1 



Mur 



M* 



r 2 
1 
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where e is a small number, x is of 0(1) and r 2 = M 3 /M 2 , then, via the see-saw, 
it is sufficient that e 2 r 2 ^> 1 in order to have approximately: 
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The determinant is naturally vanishing so that the mass eigenvalues are widely 
split and for x ~ 1 the mixing is nearly maximal. It is exactly maximal if 
x — 1. The see-saw mechanism has created large mixing from almost nothing: 
all relevant matrices entering the see-saw are nearly diagonal ||8|, J49|, that 
is they are diagonalized by transformations that go into the identity in the 
limit of vanishing e. Clearly, the crucial factorization of the small parameter 
e 2 only arises if the light Majorana eigenvalue is coupled to v^ and v T with 
comparable strength, that is x ~ 1. An interesting feature of this second 
case, in connection with a possible realization within a GUT scheme, is that it 
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requires M 3 > v 2 /m 3 , so that one can push M 3 , the scale of L violation, beyond 
the GUT scale. This is a desirable feature because, for instance, it is expected 
in SU(5) if right handed neutrinos are present and also in the breaking of 
SO (10) to SU(5). Summarizing, the second case require a peculiar hierarchy 
in the Majorana eigenvalues in order to work, but it has however the good 
characteristic of being realized with nearly diagonal matrices. 

In the next chapter we present concrete examples which realize Case I and Case 
II in a 3 by 3 context, the desired textures being motivated by the introduction of 
a suitable flavour symmetry. 



Chapter 2 

Models For Hierarchical Neutrinos 
with Large Atmospheric Mixing 



2.1 A Model with Large Mixing from mp 



By left handed mixing we mean non diagonal matrix elements that can only be 
eliminated by a large rotation of the left handed fields. Thus the question is how 
to reconcile large left handed mixings in the leptonic sector with the observed near 
diagonal form of Vckm, the quark mixing matrix. Strictly speaking, since Vckm — 
UlUd, the individual matrices U u and U& need not be near diagonal, but Vckm 
does, while the analogue for leptons, U of eq. ( |1.12| ), cannot be near diagonal. 



However nothing forbids for quarks that, in the basis where m u is diagonal, the d 
quark matrix has large non diagonal terms that can be rotated away by a pure right 
handed rotation. 

Right handed mixings for quarks can correspond to left handed mixings for 
leptons. In fact, in the context of SUSY 577(5) |5D[] there is a very attractive hint 
of how this can be realized. In the 5 of SU(5) the d c singlet appears together with 
the lepton doublet (u, e). The (w, d) doublet and e c belong to the 10 and u c to the 
1 and similarly for the other families. As a consequence, in the simplest model with 
mass terms arising from only Higgs pentaplets, the Dirac matrix of down quarks is 
the transpose of the charged lepton matrix: m d D = (m l D ) T . Thus, indeed, a large 
mixing for right handed down quarks corresponds to a large left handed mixing for 
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charged leptons because, at leading order we may have [) 
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(2-1) 



In the same simplest approximation with 5 or 5 Higgs, the up quark mass matrix is 
symmetric, so that left and right mixing matrices are equal in this case. Then small 
mixings for up quarks and small left handed mixings for down quarks are sufficient 
to guarantee small Vckm mixing angles even for large d quark right handed mixings. 
If these small mixings are neglected, we expect: 
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(2.2) 



The crucial fact is that when the charged lepton matrix is diagonalized the large left 
handed mixing of the charged leptons is transferred to the neutrinos. To see this, 
note that in SU(5) we can always diagonalize the u mass matrix by a rotation of 
the fields in the 10, the Majorana matrix M by a rotation of the 1 and the effective 
neutrino matrix m v by a rotation of the 5. In this basis the d quark mass matrix 
fixes Vqkm and the charged lepton mass matrix fixes U, the neutrino mixing matrix. 

2.1.1 The Choice of Horizontal Abelian Charges 



We give here an explicit example |32| of the mechanism under discussion in the 
framework of a unified SUSY SU(5) theory with an additional U(1)f flavour sym- 
metry |[1]. This model is to be taken as merely indicative, in that some important 
problems, like, for example, the cancellation of chiral anomalies, are not tackled 
here. Most importantly, in the following minimal version, it has the unpleasant 
feature of predicting the same spectrum for down quarks and charged leptons. The 
problem of the differentiation of these spectra will be addressed in the following 
chapters. Despite this, we find it impressive that the general pattern of all what we 
know on fermion masses and mixings is correctly reproduced at the level of orders 
of magnitude. We regard the present model as a low-energy effective theory valid 
at energies close to Mqut "C Mp\. We can think to obtain it by integrating out the 
heavy modes from an unknown underlying fundamental theory defined at an energy 
scale close to Mpi. From this point of view the gauge anomalies generated by the 



^^Here and in the following we write Dirac mass matrices with the convention LmoR- 
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light supermultiplets listed below can be compensated by another set of supermul- 
tiplets with masses above Mqut, already eliminated from the low-energy theory. 
In particular, we assume that these additional supermultiplets are vector-like with 
respect to SU(5) and chiral with respect to U(1)f- Their masses are then naturally 
expected to be of the order of the U(l)p breaking scale, which, in the following 
discussion, turns out to be near M Pi . It is possible to check explicitly the possibility 
of canceling the gauge anomalies in this way but, due to our ignorance about the 
fundamental theory, it is not particularly instructive to illustrate the details here. 

In this model the known generations of quarks and leptons are contained in 
triplets \l/" and \l/g, (a = 1,2,3) transforming as 10 and 5 of SU(5), respectively. 
Three more SU(5) singlets ^/" describe the right handed neutrinos. We start by 
discussing the Yukawa coupling allowed by U(l) ^-neutral Higgs multiplets (p 5 and 
</>§ in the 5 and 5 SU(5) representations and by the flavon X, singlet of SU(5) with 
F = -1. 

We assign to these fields the following F-charges: 



field 


#10 


^5 


*i 


X 


F 


(3,2,0) 


(3,0,0) 


(i,-i,o) 


-1 



Table 2.1: F charge assignments. 



2.1.2 Charged Fermion Mass Matrices 

In the quark sector we obtain Q: 





"A 6 A 5 A 3 " 




"A 6 


A 3 


ml = (m u D ) T = 


A 5 A 4 A 2 


v u , m d D = 


A 5 


A 2 




A 3 A 2 1 




_A 3 


1 


from which we get the order-of-magnitude relations: 






m u : m c : m t = A 6 : A 4 : 


1 






m d : m s 


:m b = A 6 : A 2 : 


1 





A 3 

A 2 
1 



v,i 



(2.3) 



(2.4) 



2 In eq. (2.3) the entries denoted by 1 in m^ and mf, are not necessarily equal. As usual, such 
a notation allows for O(l) deviations. Also remember that we write Dirac mass matrices with the 
convention LmoR- 
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and 



V,, 



A, 



V, 



ab 



x> 



V, 



cb 



y 



(2.5) 



Here v u = {(ps), Vd = (<fs) and A denotes the ratio between the vev of X and 
an ultraviolet cut off identified with the Planck mass Mp\. A = (X)/Mpi. To 
correctly reproduce the observed quark mixing angles, we take A of the order of the 
Cabibbo angle. The fermion spectrum depends on (X) and Mpi only through the 
ratio A. In chapter |] we will see how A can be dynamically generated by minimizing 
the X-dependent effective potential of the theory. For non-negative F-charges, the 
elements of the quark mixing matrix Vckm depend only on the charge differences 
of the left handed quark doublet [{51}] . Up to a constant shift, this defines the choice 

-■'"' " ' " ' -' ' We will 



in Table |2.1| . Equal F-charges for \l/g' are then required to fit mj and m^ 
comment on the lightest quark masses later on. 

At this level, the mass matrix for the charged leptons is the transpose of m d D : 



m 



D 



[m 



d \T 

d) 



(2.6) 



and we find: 



rrif, : m, 



rri-j 



A 



6 . 



A 2 :l 



(2.7) 



The 0(1) off-diagonal entry of m D gives rise to a large left handed mixing in the 
23 block which corresponds to a large right handed mixing in the d mass matrix. 



Thus, this situation corresponds to Case I of section 1.3.3 



2.1.3 Neutrino Sector 

Two different situations arise allowing or not for the presence of the additional anti- 
flavon field X, with F = 1. Let us first consider what happens allowing for the pair 
X and X. 



Two Flavons 

In the neutrino sector, the Dirac and Majorana mass matrices are given by: 



m D 



A 1 


A 2 


A 3 ] 


A 


A' 


1 


A 


A' 


1 



M 



A 2 


1 


A] 


1 


A' 2 


A' 


A 


A' 


1 



M 



(2- 



where A' = (X)/M P i and M denotes the large mass scale associated to the right 
handed neutrinos: M 3> v u d- 
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After diagonalization of the charged lepton sector and after integrating out the 
heavy right handed neutrinos we obtain the following neutrino mass matrix in the 
low-energy effective theory: 



m. 



A 6 


A 3 


A 3 " 


A 3 


1 


1 


A 3 


1 


1 



v 2 

^ (2.9) 

M V ; 



where we have taken A ~ A'. The 0(1) elements in the 23 block are produced by 
combining the large left handed mixing induced by the charged lepton sector and the 
large left handed mixing in mp. A crucial property of m v is that, as a result of the 
see-saw mechanism and of the specific U{\)f charge assignment, the determinant of 
the 23 block is automatically of 0(A 2 ). 

It is easy to verify that the eigenvalues of m u satisfy the relations: 

m 1 : m 2 : m 3 = A 4 : A 2 : 1 . (2.10) 

The atmospheric neutrino oscillations require mj ~ 10~ 3 eV 2 . From eq. fl2.9|) , 
taking v u ~ 250 GeV, the mass scale M of the heavy Majorana neutrinos turns out 
to be close to the unification scale, M ~ 10 15 GeV. The squared mass difference 
between the lightest states is of 0(A 4 ) m 2 , appropriate to the MSW solution to 
the solar neutrino problem. Finally, beyond the large mixing in the 23 sector, m v 
provides a mixing angle s\2 ~ (A/2) in the 12 sector, close to the range preferred by 
the SMA-MSW solution. In general U e ^ is non-vanishing, of 0(A 3 ). 

In general, the charge assignment under U(1)f allows for non-canonical kinetic 
terms that represent an additional source of mixing. Such terms are allowed by 
the underlying flavour symmetry and it would be unnatural to tune them to the 
canonical form. The results quoted up to now remain unchanged after including 
the effects related to the most general kinetic terms, via appropriate rotations and 
rescaling in the flavour space (see also ref. j5"2"|). 

Only One Flavon Field 

A general problem common to all models dealing with flavour is that of recovering 
the correct vacuum structure by minimizing the effective potential of the theory. 
It may be noticed that the presence of two multiplets X and X with opposite F 
charges could hardly be reconciled, without adding extra structure to the model, 
with a large common vev for these fields, due to possible analytic terms of the kind 
(XX) n in the superpotential |53[. It is then interesting to explore the consequences 



of allowing only the negatively charged X field in the theory. 
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It can be immediately recognized that, while the quark mass matrices of eq. 
fl2.3| ) are unchanged, in the neutrino sector the Dirac and Majorana matrices get 
modified into: 



m D 



A 1 


A 2 


A 3 " 


A 





1 


A 





1 



M 



A 2 


1 A" 


1 





A 


1 



M 



(2.11) 



The zeros are due to the analytic property of the superpotential that makes impos- 
sible to form the corresponding F invariant by using X alone. These zeros should 
not be taken literally, as they will be eventually filled by small terms coming, for 
instance, from the diagonalization of the charged lepton mass matrix and from the 
transformation that put the kinetic terms into canonical form. It is however inter- 
esting to work out, in first approximation, the case of exactly zero entries in mo 
and M, when forbidden by F. 

The neutrino mass matrix obtained via see-saw from mo and M has the same 
pattern as the one displayed in eq. ( |2.9| ). A closer inspection reveals that the 
determinant of the 23 block is identically zero, independently from A. This leads to 
the following pattern of masses: 



mi : ^2 : ^3 = A 3 : A 3 : 1 



0(A £ 



777 q 



Moreover the mixing in the 12 sector is almost maximal: 



s 
c 



7T 



+ 0(A 3 



f2.12) 



(2.13) 



For A ~ 0.2, both the squared mass difference (m 2 — m 2 .) /ml and sin 2 29 so [ are 
remarkably close to the values required by the VO solution to the solar neutrino 
problem. This property remains reasonably stable against the perturbations induced 
by small terms (of order A 5 ) replacing the zeros, coming from the diagonalization of 
the charged lepton sector and by the transformations that render the kinetic terms 
canonical. 

We find quite interesting that also the VO solution, requiring an intriguingly 
small mass difference and a bimaximal mixing, can be reproduced, at least at the 
level of order of magnitudes, in the context of a "minimal" model of flavour compat- 
ible with SUSY SU(5). In this case the role played by supersymmetry is essential, a 
non- super symmetric model with X alone not being distinguishable from the version 
with both X and X, as far as low-energy flavour properties are concerned. 

In chapter § we will again exploit this mechanism based on the presence of only 
one flavon X, but with slightly different charge assignments for the matter field ^/g, 
so that the LOW solution to the solar neutrino problem will naturally emerge. 
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2.1.4 Outlook 



We have constructed two examples which realize Case I of section |1.3.3| . Obviously 



the order of magnitude description offered by this model is not intended to account 
for all the details of fermion masses. Even neglecting the parameters associated 
with the CP violating observables, some of the relevant observables are somewhat 
marginally reproduced. For instance we obtain m u /m t ~ A 6 which is perhaps too 
large. However we find it remarkable that in such a simple scheme most of the 12 
independent fermion masses and the 6 mixing angles turn out to have the correct 
order of magnitude. Notice also that our model prefers large values of tan/3 = 
v u /vd- This is a consequence of the equality F(\&f ) = F(\l/g). In this case the 
Yukawa couplings of top and bottom quarks are expected to be of the same order 
of magnitude, while the large mt/rrib ratio is attributed to v u 3> v,i (there may 
be factors 0(1) modifying these considerations, of course). To keep tan/3 small, 
one could suppress m b /m t by adopting different F-charges for the \l/g and ^l Q . 
Alternatively, as will be explicitly shown in chapter |], one can assign F-charges also 
to the Higgs pentaplets. 

Additional contributions to flavour changing processes (both in the quark and 
in the lepton J5|j] sectors) and to CP violating observables are generally expected 
in a SUSY GUT. However, a reliable estimate of the corresponding effects would 
require a much more detailed definition of the theory than attempted here. Crucial 
ingredients such as the mechanism of supersymmetry breaking and its transmission 
to the observable sector have been ignored. We are implicitly assuming that the 
omission of this aspect of the flavour problem does not substantially alter our dis- 



cussion. For instance, it has been pointed out |5lJ that in SU(5) ® U(1)f based 
models with large tan/? and radiatively induced flavour mixing, one generally ex- 
pects the branching ratio of the flavour changing process ft —>■ cy to be comparable 
to the present experimental upper bound (the exact values depending on the partic- 
ular model considered). On the other hand, for tan/? = 0(1), the branching ratio 
turns out to be smaller by about four orders of magnitude. 

A common problem of all SU(5) unified theories based on a minimal Higgs 
structure is represented by the relation m l D = (rn d D ) T that, while leading to the 
successful nib = rn r boundary condition at the GUT scale, provides the wrong 
prediction m^jm s = rrie/m^ (which, however, is an acceptable order of magnitude 
equality). In section |3.5|, we will describe in some detail this problematic feature 



of minimal GUT models and briefly review the proposed solutions. An interesting 
possibility to overcome this problem and improve the picture of fermion masses will 
be adopted in the model of chapter |], where we will exploit the presence of an 
additional supermultiplet transforming in the 75 representation of SU(5). 



32 



CHAPTER 2. MODELS WITH LARGE MIXING 



2.2 A Model with Large Mixing from Nothing 



Without loss of generality we can go to a basis where both the charged lepton Dirac 
mass matrix m D and the Majorana matrix M for the right handed neutrinos are 
diagonal. In fact, after diagonalization of the charged lepton Dirac mass matrix, 
we still have the freedom of a change of basis for the right handed neutrino fields, 
in that the right handed charged lepton and neutrino fields, as opposed to left 
handed fields, are uncorrelated by the SU(2) x U(l) gauge symmetry. We can 
use this freedom to make the Majorana matrix diagonal: M~ l = V T d^iV with 
d M = Dmg[l/M 1 ,l/M 2 ,l/M 3 }. 

In order to understand how to generalize to the 3 by 3 context the mechanism 
already described in Case II of section |1.3.3j (see eq. ( |1.22| )), it is useful to start by 
considering symmetric matrices. These matrices are very interesting because, for 
example, one could want to preserve left-right symmetry at the GUT scale. Then, 
the observed smallness of left handed mixings for quarks would also demand small 
right handed mixings. 

One simple class of examples is the following one |48| 



We start from 



m D = v 



V e 



e 

e 
x e 



x e 



M" 1 = - 
A 



n 





0" 





r-2 











r 3_ 



(2.14) 



where, unless otherwise stated, x and y are 0(1); e, e' and e" are independent small 
numbers and r, = M 3 /Mj. Note that the 23 blocks of these matrices closely match 
eq. ( |1.22|) . We expect e" C e' < e <C 1 and, perhaps, also r\ ^> r 2 ^> r 3 = 1, 
if the hierarchy for right handed neutrinos follows the same pattern as for known 
fermions. Depending on the relative size of the ratios ri/rj, e/e' and e'/e", we can 
have models with dominance of any of the r\p- 

For example, we set x — 1 (keeping y of 0(1)) and assume r 2 e 2 ^> rie /2 ,r 3 , 
together with r 2 e /2 ^> rie" 2 and r 2 e ^> r\e" . 
obtain: 



In this limit, with good accuracy we 



-V 2 



m,. 



A 



r 2 e 



- 1 + 



r\e 



a 



r 2 e 



1 + 



L e 



r 3 
r 2 e 2 



(2.15) 



Note that, for ee" < e' 2 : 



DetK] = (j) 



rir 2 r 3 e 



ia 



(2.16) 
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Approximate eigenvalues, in units of v 2 /A and with m 3 > m 2 > nil, are given by 



m 3 ~ 2r 2 e , m 2 ~ —rie 



12 



rrii 



r 3 e 



a 



for rie > r 3 or, alternatively, 



m 3 ~ 2r 2 e , m 2 ~ -r 3 , 



rrii 



Tie 



/i 



(2.17) 



f2.18) 



for rit 12 < r 3 . Having set x = 1 the atmospheric neutrino mixing is nearly maximal. 
The solar neutrino mixing is instead generically small in these models, being pro- 
portional to e'/e. Thus the SMA-MSW solution is obtained. It is easy to find set of 
parameter values that lead to an acceptable phenomenology within these solutions. 
As an illustrative example we take: 



X 4 



e ~ 

-12 



A fi 



e ~ 



12 



(2.19) 

n ~ A~ 1Z , r 2 ~ A~ y , (2.20) 

where A ~ sin^c, 6c being the Cabibbo angle. The neutrino mass matrices become 



m D 



A 12 
A 6 

A 6 



A 6 


A 6 ] 


A 1 


A 4 


A 1 


1 



M = A 



A 12 











A 9 











1 



and, in units of v 2 /A, we obtain 

m 3 ~ 1/A 



m 2 



nil ~ X 4 



(2.21) 



(2.22) 



The solar mixing angle 812 is of order A 2 , suitable to the SMA-MSW solution. Also 







13 



A 2 . 



It is also possible to arrange the parameters in eq. (|2.14| ) in such a way that 
ri is dominant. Phenomenologically viable models can also be constructed in this 
case, always of the SMA-MSW type. For instance, we can take y = 1 (keeping x 
of 0(1)) in order to obtain a large atmospheric mixing angle, and further assume 
the dominance of the r^ terms, namely: e' 2 ri ^> e 2 r 2 , r 3 , with e" 2 ri ^> e /2 r 2 and 
e'Vi ^> er 2 , r 3 . In this case we can approximate the light neutrino mass matrix as 
follows: 



jii 



c'2 



e 

7 



m v 



A 



rie 



rl 



— 1 



r 2 e 



rie 
1 



rl 



(2.23) 



1 + 



rie 



a 
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As an example we consider the following choice of parameters: 



In this case we have 



e~A 4 



e ~ A , e ~ A , 



7"i ~ A 



-11 



T2 ~ A 



-6 



(2.24) 
(2.25) 



rriD = v 



A 8 


A 6 


A 6 " 


A 6 


A 4 


A 4 


A 6 


A 4 


1 



In units of v 2 /A, we find 



WI3 ~ 1/A 2 



M = A 



A 14 





0" 





A 6 











1 



rri2 ~ 1 ^> mi 



^2.26) 



^2.27) 



and 6*12 ~ A 2 . 

As anticipated in section |1.3.3| , an interesting feature of these textures, in con- 
nection with a possible realization within a GUT scheme, is that m^ is much larger 
than v 2 / A by construction. This means that the lepton number breaking scale A 
can be pushed up to the order of the GUT scale or even beyond. A large value of A, 
A > M GUT , is naturally expected in SU(5) if right handed neutrinos are present and 
in SO(10) if the grand unified group is broken down to SU(5) close to or slightly 
below the Planck mass. We recall that if 777,3 is approximately given by v 2 /A and 
v ~ 250 GeV, then A ~ 10 15 GeV. 

Models based on symmetric matrices are directly compatible with left-right sym- 
metry and therefore are naturally linked with SO( 10). This is to be confronted with 
models - like those described in the previous section - that have large right handed 
mixings for quarks, which, in SU(5), can be naturally translated into large left 
handed mixings for leptons. In this connection it is interesting to observe that the 
proposed textures for the neutrino Dirac matrix can also work for up and down 
quarks. For example, the matrices 



nip oc 






A 6 


A 6 " 


A 6 


A 1 


A 4 


A 6 


A 1 


1 



m D oc 






A 3 


A 3 " 


A 3 


A 2 


A 2 


A 3 


A 2 


1 



(2.28) 



where for each entry the order of magnitude is specified in terms of A ~ sin#c*, lead 
to acceptable mass matrices and mixings. In fact m u : m c : m t = A 8 : A 4 : 1 and 
A 4 : A 2 : 1. The Vckm matrix receives a dominant contribution from 



m d :m s : m b 



2.2. A MODEL WITH LARGE MIXING FROM NOTHING 



35 



the down sector in that the up sector angles are much smaller than the down sector 
ones. The same kind of texture can also be adopted in the charged lepton sector Q 

It is simple to realize that models of this sort cannot be derived from a broken 
U(1)f horizontal symmetry, with a single field that breaks spontaneously U(1)f- For 
symmetric matrices, left and right U(l) charges must be equal: then, for example, 
if m 33 ~ l 2q3 ~ 1 or g 3 = 0, then m 22 ~ l 2q2 is different from m 23 ~ l q2 etc. 

However, the symmetry requirement is not really necessary to produce a large 
atmospheric mixing angle from neutrino mass matrices characterized by small mix- 
ings. Indeed, by considering for simplicity the 2 by 2 case, instead of eq. Q1.22J ), we 
could equally well have started from a Dirac mass matrix of the kind: 



mo oc 



e 
xe 



x 
1 



M" 1 



x 



r 2 
1 



(2.29) 



with x' any number smaller than or equal to one. Once the condition e 2 r 2 ^> 1 is 
satisfied, the neutrino mass matrix of eq. ( |1.23|) is obtained again, independently 
from x 1 . The generalization of this mechanism to the full 3 by 3 case is straight- 
forward, and analogous to the one discussed above for symmetric matrices. This 
leaves much more freedom for model building and does not necessarily bound the 
realization of these schemes to left-right-symmetric scenarios. 

The solar mixing angle is generically small in the class of models explicitly dis- 
cussed above. However, small mixing angles in the Dirac and Majorana neutrino 
mass matrices do not exclude a large solar mixing angle. For instance, this is gen- 
erated from the asymmetric, but nearly diagonal mass matrices: 



m D = v 



A 6 





0" 


A 6 


A 1 








A 1 


1 



M = A 



A 12 





0" 





A 10 











1 



They give rise to a light neutrino mass matrix of the kind: 



(2.30) 



m v 



A 2 


A 2 


A 2 


1 





1 



A 2 A 



which is diagonalized by large 6* 12 and # 23 and small #13. 
suitable to the LMA-MSW solution. 



(2.31) 



The mass hierarchy is 



3 Notice that the possibility of treating on equal footing all fermion mass matrices is also present 
in the framework of an extended flavour democracy. See for instance Ref. |5q] and references 
therein. 
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2.2.1 Renormalization Effects 

These results are stable under renormalization from the high-energy scale A where 
the mass matrices are produced down to the electroweak scale. Indeed, suppose 
that we start with a supersymmetric theory whose superpotential w at the scale 
A > Mqut is given by: 



w = Qy u U c H u + Qy d D c H d + Ly v N c H u + Ly e E c H d + -N C MN C 



where 



o o 

y T 



y"(A) 



e x e 
x e 1 



M(A) 



M 2 
M 3 



(2.32) 



(2.33) 



Here we focus on the two heaviest generations. The matrices y u ^ e ' v describe the 
Yukawa couplings of the theory, related to the mass matrices through m UyU = y u,u v u 
and ra<z, e = y d ' e Vd, where v Ujd are the vacuum expectation values of the neutral scalar 
components in the superdoublets H ud . The mass scales are ordered according to 
A > M 3 ^> M 2 3> Myy, My/ denoting the electroweak scale. At the scale M 2 , after 
integrating out the superfields A^ c associated to the heavy right handed neutrinos, 
the mass matrix for the light neutrinos, including the renormalization effects, is 
given by: 



mJM 2 ) = — 

vy 2J M 2 



\2 



°w 



^2.34) 



j/22 2/32 {y h 32j 

The couplings y u 22 an d 2/^32 are evaluated at the scale M 2 . If the running of M is 
computed by by taking e = 0, then M 2 stays approximately constant when going 
from A down to M 2 . Notice that, if we neglect the terms of 0(1/M 3 ), the determi- 
nant of m u (M 2 ) vanishes, even including the renormalization effects. These effects 
modify the mixing angle at the scale M 2 and they depend on how the Yukawa cou- 
pling y" 2 2 an d y U 3 2 change from A to M 2 . The running of the neutrino Yukawa 
couplings is governed by the equation: 



dt 



167T 2 



tr(3y u y^ + y^) 



-47r(3a 2 + -a 1 ) + 3yV t +2/V t 
5 



V" , (2-35) 



from /i = A to ft = M 3 , with t = log//. At M 3 , the superfield N£ is integrated out. 
This produces terms of order 1/M 3 , which are neglected in the present discussion. 
From ft = M 3 down to fj, — M 2 the running of y u 22 and y u 32 is described by an 
equation formally identical to the previous one, with 



y = 



2/22 2/23 
2/32 2/33 



y 



2/22 
yk 



(2.36) 
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To assess the size of the renormalization effects on y u 22 and y u 32 , we have integrated 
numerically the above equations in the approximation of constant y u,e and «i j2 . This 
should provide a rough estimate of the effects. As an example we take A = 10 18 
GeV, M 3 = 10 16 GeV and M 2 = 10 8 GeV, a 2 = (5/3)a x = 1/24; we choose y T = 0.6, 
corresponding to the large tan j3 = v u /vd regime, which enhances the renormalization 
effects of y u . Assuming as initial conditions y y 22 (A) = 2/^23 0V) = y v '32(A) = 0.05 
and y u 33 (A) = 1, we find that y l ' 22 and y u 32 are modified by about 10% at the scale 
M 2 . The mixing angle, nearly maximal at the scale A, remains close to maximal: 
sin 2 2^23 ~ 0.98 in the example at hand. 

The running from M 2 down to the electroweak scale does not appreciably affect 
neither the vanishing determinant condition, nor a nearly maximal neutrino mixing 
angle. This can be seen by analyzing the renormalization group equation for m u , 
which, in the one-loop approximation, is given by: 



dm. 



dt 8vr 2 
The solution is: 



tr(3y"V) - 47r(3a 2 + -ai) 



m, 



+ - [m u y e] y e + (y e] y e ) T m u ] \ 



m v {M ] 



w j 



C 



m v22 (M 2 ) m u32 (M 2 ) B 
m v32 (M 2 ) B m v33 (M 2 ) B 2 



(2.37) 



(2.38) 



where 
C = exp 

and 



c{t')dt' 



to 



B = exp 



b(t')dt' 



to 



c(t) 



■iTl" 



t = log(M w ) 



tr(3j/ u V)- 47r(3a 2 + - ai ; 

5 



(*) , 



b(t) 



y 2 T (t) 

16vr 2 



t = log(M 2 ) 
(2.39) 

(2.40) 
(2.41) 



We see from eq. ( |2.38|) that if the determinant of m u vanishes at the scale M 2 , then 
it vanishes also at the weak scale. The mixing angle at the weak scale is given by: 



sin 2 (26) (Mi 



w , 



4m,| 2 (M 2 )S 5 



4m u 2 32 (M 2 )B 2 + (m u33 (M 2 )B 2 - m u22 (M 2 )y 



(2.42) 



The condition of maximal mixing at the scale M 2 is m l/33 (M 2 ) = -m^ 22 (M 2 ). It is 
easy to see that, if this condition is met, then the first correction to sin 2 (26 l )(Mvi/) is 
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of second order in the parameter (B — 1). Numerically, we find that sm 2 (29)(Mw) 
is reduced only by a of few percent. 

Summarizing, the renormalization effects do not appreciably modify the results 
obtained at the tree level for eigenvalues and mixing angle of the neutrino mass ma- 
trix. They manifest themselves as 0(1) coefficients that, for all practical purposes, 
can be absorbed in the definition of the classical parameters. An analogous analysis 
can be carried out in the 3 by 3 case, with similar results. This is in agreement 
with the analysis of ref. |57|], as far as the case of hierarchical neutrino masses is 
concerned. 



2.2.2 An Explicit GUT Example with Broken Flavour Sym- 
metry 

From our general discussion it is clear that some amount of correlation between the 
Dirac and the Majorana sectors is required in order to reproduce, at the order-of- 
magnitude level, the observed pattern of oscillations (assuming from now on the 
SMA-MSW solution for solar neutrinos). In order to show that this correlation 
could be induced by an underlying flavour symmetry and, at the same time, to 
present an explicit grand unified example of the class of textures considered here, 
we [Q sketch the features of an SU(5) model with a spontaneously broken flavour 
symmetry based on the gauge group U{1)a x U{1)b- 

The three generations are described by \l/io a and \l/g a , (a = 1, 2, 3) transforming 
as 10 and 5 of SU(5), respectively. Three more SU(5) singlets \l/i a describe the 
right handed neutrinos. Here we focus on the Yukawa coupling allowed by Higgs 
multiplets <p 5 and </>§ in the 5 and 5 SU(5) representations and by three flavons, ip^, 
ip B and (/?£, taken, to begin with, as SU(5) singlets. The charge assignment of the 
various fields is summarized in Table |2.2| . We assume that the flavon fields develop 
the following set of vevs: 



(<p A ) ~ A A 



(<p B ) ~ A A 



(<Pb) ~ A 



(2.43) 



where A denotes the cutoff of the theory and provides the mass scale suppressing 
the higher-dimensional operators invariant under the gauge and the flavour groups. 



In the quark sector we obtain ^ 

■a 8 a 6 a 4 



m 



D 



(m 



i)i 



A 6 
A 1 



A b 

A 1 
A 2 



A 2 
1 



v„ 



m 



D 



A 6 


A 1 


A 3 


A 5 


A 3 


A 2 


A 3 


A 


1 



Vd 



(2.44) 



4 In eq. fl2.44 ) the entries denoted by 1 in m% and m d D are not necessarily equal. As usual, such 
a notation allows for 0(1) deviations. 
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*10 


*5 


*i 


V?5 


¥5 


Va 


<Pb 


¥ + B 


5*7(5) 


10 


5 


1 


5 


5 


1 


1 


1 


u(i) A 


(3,2,0) 


(3,1,0) 


(6,3,0) 








-1 








U(1) B 


(1,0,0) 


(-2,-2,-1) 


(0,2,0) 











-1 


+ 1 



Table 2.2: Quantum numbers of matter and flavon multiplets. 



from which we get the order-of-magnitude relations: 

m u : m c : m t =- A 8 : A 4 : 1 
m d : m s : m h 

and 



A 6 : X s : 1 



V u , 



(2.45) 



A, V uh ~\ 6 , V cb ~\ z . (2.46) 

Here v u = (</?s), t>d = (v^s)- Notice that V u b and V^ s are dominated by the down- 
quark contribution. The mass matrix for the charged leptons is the transpose of 



m: 



D- 



m l D 



[m 



d \T 

d) 



and we find 



m e : m u '■ va T = A 6 : A 3 : 1 



(2.47) 



f2.48) 



At this level we obtain the well-known prediction m^ = m T , together with the 
unsatisfactory relation m^/m s = m e /m^ (which, however, is an acceptable order of 
magnitude equality). 

In the neutrino sector, the Dirac and Majorana mass matrices are given by: 



m D 



A 9 


A 6 


A 3 ] 


A 7 


A 1 


A 


A 6 


A 1 


1 



M 



A 12 
A 11 


A 11 
A 10 


A 6 " 
A 5 


A 6 


A 5 


1 



A 



(2.49) 



Notice that while in the model under consideration m^ is asymmetric, it is diagonal- 
ized by small unitarity transformations, i.e. transformations that go to the identity 
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in the limit of vanishing A, in both the left and the right sectors. After diagonaliza- 
tion of the charged-lepton sector and after integrating out the heavy right handed 
neutrinos, we obtain the following neutrino mass matrix in the low-energy effective 
theory: 

^A 4 A 2 A 2 ^ 
in,. = A 2 1 1 
A 2 1 1 



v 2 
A 2 ! 



(2.50) 



The 0(1) elements in the 23 block are produced only by the interplay of the left 
handed mixing in mjj and the hierarchy in the Majorana sector. The contribution 
from the charged lepton sector is subleading. An important property of m u is that, 
as a result of the see-saw mechanism and of the specific U(l)p charge assignment, 
the determinant of the 23 block is of 0(A 2 ) (in units of (t> 2 /A 2 A) 2 ). This gives rise 
to the desired hierarchy between atmospheric and solar frequencies. It is easy to 
verify that the eigenvalues of m v satisfy the relations: 

mi : m 2 : m 3 = A 8 : A 2 : 1 . (2.51) 

The atmospheric neutrino oscillations require m 2 , ~ 3.5 x 10 -3 eV 2 . From eq. ( f2.50[ ), 
taking v u ~ 250 GeV, the mass scale A of the heavy Majorana neutrinos turns out to 
be close to 2- 10 16 GeV. The squared mass difference between the lightest states is of 
0(A 4 ) m 2 , appropriate to the MSW solution of the solar neutrino problem. Finally, 
beyond the large mixing in the 23 sector, m v provides a mixing angle 9 12 ~ A 2 in 
the 12 sector, close to the range preferred by the SMA-MSW solution. In general 
U e s is non-vanishing, of 0(A 2 ). We find it encouraging that the general pattern of 
all that we know on fermion masses and mixings is correctly reproduced at the level 
of orders of magnitude. 

This model is to be taken merely as an illustration, among many other possi- 
bilities, of the scenario outlined: small mixing angles for the fundamental fermions 
giving rise to a large mixing angle for the light neutrinos. For this reason we do 
not attempt to address some important problems, such as the cancellation of chiral 
anomalies, which we implicitly postpone to an energy scale higher than the unifi- 
cation scale. We have not dealt here with the problem of recovering the correct 
vacuum structure by minimizing the effective potential of the theory. Also, the or- 
der of magnitude description offered by this model is not intended to account for all 
the details of fermion masses. Even neglecting the parameters associated with the 
CP violating observables, some of the relevant observables are somewhat marginally 
reproduced. For instance, as already emphasized, a common problem of all SU(5) 
unified theories, based on a minimal Higgs structure, is represented by the relation 
m l D = (m d D ) T , which, while leading to the successful m^ = m T boundary condi- 
tion at the GUT scale, provides the wrong prediction m^/m,, = m e /m^. We might 
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overcome this problem and improve the picture [57j by assuming that the flavon 
fields (p^, ip% or a suitable subset of them transform in the adjoint representation 
of 577(5). The product (24) q 5 (where q is a positive integer) contains both the 5 
and the 45 S77(5) representations and the couplings would allow for a differentiation 
between the down quarks and the charged leptons. Again, since the purpose of the 
present model was only illustrative, we do not insist here on recovering fully real- 
istic mass matrices. Finally, additional contributions to flavour- changing processes 
and to CP-violating observables are generally expected in a supersymmetric GUT. 
However, a reliable estimate of the corresponding effects would require a much more 
detailed definition of the theory than attempted here. Crucial ingredients such as 
the mechanism of supersymmetry breaking and its transmission to the observable 
sector have been ignored in the present note. 

In this example we have given up a possible symmetric structure of the textures 
in favour of a relatively simple and economic realization of the underlying flavour 
symmetry. It is possible, within the 577(5) or 5*0(10) gauge theories, to provide 
slightly more elaborate examples where this symmetric structure is manifest. Here 
we have insisted on requiring both the Dirac matrices and the Majorana matrix 
M to be quasi- diagonal (that is with all off-diagonal entries ij suppressed with 
respect to the largest of the diagonal ii and jj ones). In some alternative models, 
all Dirac matrices are indeed quasi-diagonal, but large mixings are present in the 
Majorana matrix M. For instance, the idea of a large atmospheric mixing angle from 
intrinsically small mixings in the Dirac mass matrices has recently been advocated 
in the context of a U(2) flavour theory [BS|. In this model the Dirac neutrino mass 



matrix is the one displayed in eq. (|2.14j) , with y of 0(e) (i.e. {7x10)12, <S (^0)12) and 
e" ~ e /2 . The U(2) symmetry favours an r% dominance, by adequately suppressing 
the M11 element. However, to have both v^ and v T with approximately the same 
coupling to the light Majorana mass eigenstate, this should possess a component of 
order e' along the v\ direction. In turn, this property is only exhibited if also the 
M33 entry is sufficiently small []. As a consequence, there are two heavy Majorana 
states with similar masses and there is a large mixing in the 23 Majorana sector. In 
this case, large mixings, which are absent from the Dirac mass matrices, arise from 
the Majorana sector, a physically similar possibility but technically different from 
that examined here. 



technically, this suppression can be realized in the framework of an 5*0(10) gauge group, by 
assuming that there are no flavour singlets that couple to v^v^ . 
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2.2.3 Outlook 

Summarizing, in most models that describe neutrino oscillations with nearly max- 
imal mixings, there appear large mixings in at least one of the matrices run, Tn l D , 
M (i.e. the neutrino and charged-lepton Dirac matrices and the right-handed Ma- 
jorana matrix). Here we have discussed the peculiar possibility that large neutrino 
mixing is only produced by the see-saw mechanism starting from all nearly diagonal 
matrices. Although this possibility is certainly rather special, we have shown that 
models of this sort can be constructed without an unrealistic amount of fine tuning 
and are well compatible with grand unification ideas and the related phenomenology 
for quark and lepton masses. 



Chapter 3 



Problems of SUSY GUTs 



The idea that all particle interactions merge into a unified theory at very high 
energies p| is so attractive, for both aesthetical and phenomenological reasons, 



that this concept has become widely accepted by now. Indeed, the quantitative 
success of coupling unification in SUSY Grand Unified Theories (GUTs) |K)| can be 



hardly interpreted as just a coincidence. In addition, support has also come from the 
developments on neutrino oscillations |l[ , which point to lepton number violation at 
a scale close to the coupling unification one. 

Despite the fact that they constitute promising ideas beyond the standard model, 
SUSY GUTs are in a not so safe ground as at first sight could appear. This chapter 
is devoted to a critical reanalysis of those problematic aspects found in constructing 
realistic models. 

Assuming gauge unification, minimal models of GUTs based on SU(5), 5*0(10), 
... have been considered in detail. Also, many articles have addressed particular 
aspects of GUTs models like proton decay, fermion masses and, recently, neutrino 
masses and mixings. In the previous chapter, indeed, the attention has focused on 
the phenomenology of neutrino masses and mixings. But on the one hand, minimal 
models are not plausible as they need a large amount of fine tuning and are therefore 
highly unnatural, for example with respect to the doublet-triplet splitting problem 
or proton decay. On the other hand, analyses of particular aspects of GUTs often 
leave aside the problem of embedding the sector under discussion into a consistent 
whole. 

So the problem arises of going beyond minimal toy models by formulating suffi- 
ciently realistic, not unnecessarily complicated, relatively complete models that can 
serve as benchmarks to be compared with experiment. More appropriately, instead 
of "realistic" one should say "not grossly unrealistic" because it is clear that many 
important details cannot be sufficiently controlled and assumptions must be made. 
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Finally, the success or failure of such program can decide whether a stage of gauge 
unification is a likely possibility or not. 

In this chapter we will review SUSY GUT's basic weak-points and analyze crit- 
ically the current status of solutions. As will be seen, most of these solutions are 
partial, have some problematic feature and are in general quite complicated. The 
crucial problem is to find a consistent picture of the Higgs sector which leads to 
light doublets but triplets heavy enough to escape current bounds on proton decay 



6jJ . In addition, the field content must not destroy coupling unification. During 
this survey, we will devote special attention to those solutions that will adopted in 
the following chapter, where an example of semi-realistic model will be proposed. 



3.1 The Doublet-Triplet Splitting Problem 

The phenomenological motivation for weak scale SUSY is the solution of the hier- 
archy problem. While the mass of the Higgs doublet should be around 100 GeV, 
both coupling unification and non observation of higgsino mediated proton decay 
require the Higgs triplet mass to lie around 10 16 GeV, the precise value depending 
on the details of the model under discussion. Achieving this doublet-triplet splitting 
B2j is the most compelling problem one has to face in GUTs, because it requires to 



fine-tune a parameter at the level of 14 orders of magnitude. It is worth to remind 
that, with respect to this problem, SUSY GUTs stay in a much better situation 
than non SUSY GUTs. Thanks to non-renormalization theorems, in the former the 
fine-tuning has to be done only once, that is at tree level, while in the latter it has 
to be done at all orders of the perturbative expansion. 

After a short presentation of the transparent way in which the doublet-triplet 
splitting problem appears in the case of minimal S77(5), this section has the pur- 
pose of summarizing the status of the most interesting mechanisms which avoids 
this fine-tuning also insisting on their weak-points. In many models, in fact, even 
if one succeeds in achieving the desired splitting at tree level, this could be upset 
by a lifting of the doublet mass, due to either radiative corrections or non renor- 
malizable operators. In addition, one has to be careful because some of the possible 
mechanisms can destroy coupling unification. 



3.1.1 The Doublet-Triplet Splitting Problem in Minimal 577(5) 

In minimal SUSY 577(5) the Higgs sector consists of one chiral superfields in the 
adjoint representation, E, and one pair of fundamental and antifundamental, H, H. 
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The breaking of 577(5) is provided by the following superpotential terms: 

w m99a = M 24 E^E* + A E£E*E$ + c H \Y, A B H B + M 5 H A H A , (3.1) 

where A,B,... = 1, ..., 5 and summation over repeated indices has to be understood. 
In the limit of unbroken SUSY, minimization of this part of superpotential gives the 
following SU(3) (8> SU(2) ® U(l) symmetric vevs Q: 

(E) = |^ ^(2,2,2,-3,-3) , (3.2) 

(77) = , (3.3) 

(77) = . 

This superpotential also ensure heavy masses to those £ multiplets which are not 
eaten by the Higgs mechanism (see Table |3.7|) . Due to the structure of the £ vacuum, 
doublets and triplets acquire different masses: 

M D = M 5 --^M 24 , M T = M 5 + ^-M 2i . (3.4) 

A 3A 

To have light doublets one has to impose that the parameters M 5 and M 24 , which 
should assume values around 10 16 GeV, satisfy 

M 5 = ^M 24 (3.5) 

with an accuracy of at least the 14 th decimal, so that Mr, = O(lOOGeV). This 
hierarchy problem is inherent to all GUTs. SUSY doesn't explain the origin of the 
two scales, but renders this constraint stable against radiative corrections, so that, 
at least, the fine-tuning has to be done only once. Non SUSY GUTs, in contrast, 
need the tuning to be repeated at all perturbative orders. 

Similar considerations also apply to the case of minimal 5*0(10), where 77, H lie 
in the 10 representations of 50(10). 

3.2 Solutions to the Doublet-Triplet Splitting Prob- 
lem 

The actual problem is to achieve the doublet-triplet splitting in a natural way, 
without fine-tuning. Some solutions that overcome this problem have been proposed, 

1 Actually there are three supersymmetric degenerate minima: in one of them ST/ (5) is unbroken, 
while in the other two it is broken respectively to SU(4) <g> 17(1) and SU(3) ® 517(2) <g> 17(1). 
The degeneracy could eventually be lifted and the MSSM vacuum chosen due to supergravity 
corrections. 
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based on the introduction of Higgs representations with useful group-theoretical 
properties, which do the job of allowing a mass only for triplets but not for doublets. 
Several attempts have been done and someone of the most important will be 
summarized in the following. Despite the fact that basic ideas are indeed very ele- 
gant, emphasis will be given to some problems which remain unsolved. In addition, 
it turns out that embedding these solutions in a realistic and not too cumbersome 
framework is not an easy task. 

3.2.1 The Missing Doublet Mechanism 

One of the most economical solution the doublet-triplet splitting problem is repre- 
sented by the Missing Doublet Mechanism (MDM) |63|| . It works in SU(5), which is 
the lower rank group containing the standard model, and is realized just replacing 
the 24 of Higgses by a 75, Y, and by further adding a pair of 50, 50 Higgses. Unfor- 
tunately, this proposal suffers of some weakness which, however, can be overcome. 

Description of the Mechanism 

The part of the superpotential responsible for the breaking of SU(5) reads 

w„ iggs = Cl Y& YEF YfE + M Y Y£E YlE , (3-6) 

where, as usual, A,B,... = 1,2, ...,5. These terms provide Y with a vev (Y) = 
My/(2ci) ~ Mqut and give a mass to all physical components of Y, i.e. those that 



are not absorbed by the Higgs mechanism (see Table |4.2|) . The SUSY vacuum of 



the 75 breaks SU(5) uniquely to SU(3) <g> SU(2) <g> U(l): 



{Y£S) = ~(Y) 



6«8§ + 28 a c 8 b d -h$8E) (A- B) 



(3.7) 



where Greek indices assume values from 1 to 3 and Latin indices from 4 to 5. 
The doublet-triplet splitting is realized by the following superpotential terms, 

_„ uA v DEtt BCFG c i _ T7 vBCtt ADEFG 

WmDM — C2 tl r BC U 50 £ADEFG + c 3-n A* DE n WBCFG e 

+ M 50 #50 ' H mABCD . (3.8) 

Note that the direct term M 5 HaH a must not to appear. The splitting occurs be- 
cause the 50 of chiral superfields contains a (3, 1, 1/3), i.e. a coloured antitriplet, 
SU{2) singlet (of electric charge 1/3) but no colourless doublet (1,2,-1/2). Denot- 
ing the colour triplets contained in H, H, H 50 , H-^ by H Su , H 3 d, H' 3d and H' 3u 
respectively, after SU(5) breaking: 

wmdm 3 -4V3c 2 (Y) H' 3d H 3u - 4v / 3c 3 (Y) H 3d H' 3u + M 50 H' 3d H' 3u + ... (3.9) 
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where dots stand for mass terms for the other multiplets contained in 50, 50. Thus, 
wmdm doesn't produce any mass term for the doublets but, on the other hand, the 
Higgs colour triplets mix with the analogous states in the 50. The resulting triplet 
mass matrix is of the see-saw form: 



rriT 



-AVSc 2 (Y) 

-AVSc 3 (Y) m 50 



(3.10) 



,t 



The eigenvalues of the matrix m^viVj, are the squares of: 
1 



Mf + 48 (c 2 + c 3 ) 2 (Y) 2 ± \ M 2 Q + 48 (c 2 - c 3 ) 2 (Y) 



(3.11) 



Note that mnin^ = 48 c 2 c 3 (Y) 2 , so that triplets are naturally heavy. The 
effective mass that enters in the dimension 5 operators with \AB = 1| is the following 
combination of physical masses: 

m Tl m T2 48 c 2 c 3 (Y) 2 

mr = — 7T = 7T • 3.12 

M 50 M 50 v ; 



In the literature there is confusion |64], ^, ^6], ^7| on the numerical coefficients 
of the previous relations and also on the values for the heavy masses of the 75. In 
particular, many authors seem to forget the coefficient 4v^3 in eqs. Q3.9Q , ( |3.10| ) 



and following. But this factor is important, because it gives in eq. ( |3.12|) the large 
coefficient 48 that cannot be neglected when evaluating the effective triplet mass 
which suppress proton decay amplitudes. The formulae above have been rederived, 
by carefully paying attention to the group-theoretical factors. Appendix |C| contains 
a brief summary of this procedure. 

Possible Criticisms 

Criticism against the MDM is often moved by the common fear that its large rank 
representations are dangerous, mainly for two arguments. The first one is that, from 
the point of view of string theory, it is really hard to find examples containing such 
huge representations. From a more conservative point of view, however, embedding 
GUTs into a string theory is for sure interesting but not at all compelling. The sec- 
ond argument is that the large content of matter fields in the MDM causes the gauge 
coupling to blow up before Mp\. In principle, a theory which is non perturbative 
before Mpi is acceptable. 

Apart from the two previous considerations, the MDM in its simplest version 
has mainly a weak point. No matter what symmetry is used to forbid the direct 
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mass term for HH in the superpotential, the doublet Higgs could take mass from 
non renormalizable operators of the form [BIJ : 



HHYY 

^r ■ (313) 

If one believes Planck suppressed operators consistent with the theory are present, 
the doublet has a too big mass. The problem is exacerbated in MDM because the 
coupling blows up at a low scale and it is probably the associated strong scale which 
would suppress such operators. A promising way to overcome this obstacle was 
introduced in ref. |67j and is based on the introduction of an additional anomalous 



U(l) flavour symmetry. 

As will be clarified, the solution of the doublet-triplet splitting problem in terms 
of the 50, 50 and 75 representations is very suggestive because it automatically leads 
to improve the prediction of a^{mz) and at the same time relaxes the constraints 
from proton decay. 

3.2.2 The Sliding Singlet Mechanism 

New possibilities emerge when introducing an additional singlet field S. The simpler 
one |)9| is the Sliding Singlet Mechanism. It is realized by just adding to eq. ( |3.1|) 



the term 

fSH A H A . (3.14) 

where / is an adimensional coefficient. At tree level and considering also soft SUSY 
breaking terms, the vev of S can be arranged to cancel the large contributions to 
the doublet mass coming from the other terms. However, it is well known that this 
mechanism is unstable against radiative corrections [7(J, because tadpole diagrams 



induce the doublet to take a too big mass ~ ^/rn^pMcur- 
3.2.3 The Missing VEV Mechanism 



The Missing VEV Mechanism |7T| (MVM), which apply to the case of 50(10), 
is based on the presence of a 45 Higgs representation. With respect to SU(5), 
45 = 1 + 24 + 10 + 10, so that, in practice, having a singlet and a 24, it is based on 
the same idea of the Sliding Singlet. This can occur because the 45 representation 
is not constrained to be traceless. In fact, the 45 can have a vev of the form 
(45) ~ diag(l, 1, 1,0,0) (g> ia. The term 10i 45 IO2, gives mass to triplets but not 
to doublets. Note that two 10 are required due to the antisymmetry of 45, so that 
there are four doublets and four triplets. 
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This mechanism seems very nice at first glance, but worrisome at second. In fact 
the problem arises of giving mass to just two doublets (or two linear combinations 
of them). This could be done, for instance, by introducing a term M2IO2IO2 and 
forbidding any direct coupling with IO1IO1. Heavy doublets have then a mass ~ M2 
and are in no way connected with the light ones. But there is also the need for 
triplets heavy enough to agree with experimental bounds on proton decay, and this 
is not very easy. For instance, by adding M2IO2IO2, the triplets mass matrix % 
turns out to be of a see-saw form (exactly as in the case of the MDM but with the 45 
instead of the 75 of SU(5), see eq. (|3.12|) ). The effective mass entering proton decay 
amplitudes is then proportional to (45) 2 /M 2 . But in order to achieve the needed 
suppression of proton decay amplitudes, we need M 2 < 1CT( 2 ^ 3 )(45) ~ M GUT . Apart 
from the fact that the not-so-heavy doublets would give dangerously large threshold 
corrections, this requirement is clearly a fine-tuning. 

The simplest versions of the MVM are then not realistic, but at the price of 
introducing a good deal of new fields, various models have been realized in the past 
in which the parameter playing the role of M2 is suppressed because of symmetry 
reasons. Also the new fields produce large threshold corrections, so that one cannot 
exclude the accidental possibility of a general compensation of all their contributions. 

However, dealing with a large number of GUT fields is not very attractive, be- 
cause the predictivity of the model is highly reduced and, most importantly, the 
success of coupling unification within the MSSM would be merely a lucky coin- 
cidence. In addition, due to the increased number of Higgs representations that 
generally appears in many models, there is the problem of actually implementing 
a potential to get the desired minimum, paying attention to possible flat directions 
and without destabilizing the minimum with higher order terms. In fact, in order 
to lower the rank of 5*0(10), one needs also one (or more) 16. But the couplings 
between 16 and 45 can destabilize the the desired vev of the 45. The literature is 



rich of interesting models ]72, [T^, |73] where this sort of compromise between the 



internal elegance and the necessity of agreeing with phenomenology appears clearly. 



3.2.4 Higgses as pseudo-Goldstone Bosons 

An interesting approach to the doublet-triplet splitting problem is obtained if one 
consider the Higgs doublets to be pseudo-Goldstone bosons related to a sponta- 
neously broken accidental global symmetry of the Higgs sector. This distinguishes 
the doublets from the triplets in a very non trivial way. In this scenario the Higgs 
sector of theory does possess a spontaneously broken accidental global symmetry 
and the Higgs doublets are the associated Goldstone bosons. The accidental global 
symmetry is explicitly broken in the Yukawa sector. The Higgs doublets become 
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pseudo-Goldstone bosons with a mass determined by the scale of soft SUSY break- 
ing terms. Because of non-renormalization theorems, this mass is protected against 
radiative corrections. 

Despite that the idea is very elegant, the difficulty one finds in its actual realiza- 
tion is to motivate the accidental symmetry in a natural way. For instance, models 
involving a gauged SU(5) symmetry and a global S77(6) |75| unfortunately require 
an unpleasant fine-tuning of some parameters of the superpotential. In order to 
justify the accidental symmetry, an alternative possibility J7B| is to assume that it 



arises because two sectors of the chiral superfields responsible for gauge symmetry 
breaking do not mix. But such a separation has to be imposed ad hoc. 



3.3 Proton Decay 

Proton decay is maybe the most interesting consequence of GUTs. Since a long 
time it has been representing a very hard challenge from the experimental point of 
view. The Kamiokande detector was construct for this purpose and, after Super- 
Kamiokande, maybe a third-generation detector will be constructed. From the the- 
oretical perspective this subject is not less engaging: the actual value of proton 
lifetime (if any) is determined by such an enormous amount of physics, from the 
GUT scale down to 1 GeV, that it offers the possibility of exploring various aspects 
of particle physics at the same time. The theoretical prediction is thus not at all 
accurate and strongly model dependent, but the common opinion is that it should 
lies just in the range already at hand of experimental searches or not very far. 
A lot of weighty informations have already come from these studies: 

• non SUSY SU(5i) is excluded because, due to the lightness of heavy gauge 
bosons, it predicts a lifetime induced by dimension 6 operators which is much 
lower than present bounds. Same conclusions hold for greater groups which 
break in one step to SU(5), like 50(10). 

• minimal SUSY SU(5) is in very serious troubles: with SUSY particles lighter 
than 1 TeV, dimension 5 operators would induce proton to decay mainly into 
K + i? at an exceedingly fast rate, in conflict with experiment. The case of 
SO (10) is less definite, because the prediction of proton lifetime is linked to 
the particular Higgs content chosen to realize the breaking of 50(10). How- 
ever, it turns out that also in those models where as few as possible Higgs 
representations are introduced flJ2| , the prediction of proton lifetime should 
not be far from the range explored by present searches. 
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Channel 


Lower Limits (90 % CL) 
r/B(10 33 yr) 


K+U 

7T + V 

e + vr° 

/i+TT 


1.9 

0.025 

4.4 

3.3 



Table 3.1: Lower Limits on Proton Decay Channels [77]. 



Prompted by the recent Super-Kamiokande improvements on proton partial life- 
times bounds J77[, a renewed interest on this subject has arisen. Since GUTs without 
SUSY turned out to be disappealing, the attention has focused mainly on the ques- 
tion whether it is possible to construct examples of realistic SUSY GUT models. 
In such models, achieving a sufficient suppression of proton decay amplitudes from 
dimension 5 operators is of crucial importance. In fact, since the actual prediction 
of proton lifetime depends significantly of the details of the model, one is interested 
in arranging things in such a way to escape conflict with bounds. 

To understand how this can be done reasonably, that is without introducing fine- 
tuning, it is useful to have a close look at the calculation of proton decay lifetime 
in the simplest possible case, namely minimal SUSY SU(5). In this framework 
we will summarize the principal steps of the calculation of proton lifetime induced 
by dimension 5 operators [|79| . We will treat independently Yukawa couplings of 
doublets and triplets. Also down quarks and charged leptons mass matrices will 
not be correlated, so that the calculation can be easily adapted also to non minimal 
models. We will show in some detail that, generically, the amplitude is dominated 
by only two terms of the Lagrangian and we will discuss, also quantitatively, how 
large can be the effect of the two CP violating phases (in addition to the CKM 
one) which enter this amplitude. We will comment on the source of the various 
uncertainties and give numerical examples. These have to be compared with the 
bounds collected in Table |D] from ref. |77| . 
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3.3.1 Dimension Five Operators 

The starting point is the most general renormalizable superpotential which could 
come from SU(5). Yukawa couplings of doublets and triplets will be considered un- 
correlated as well as down quarks and charged lepton mass matrices. The discussion 
can thus be adapted for any SU(5) based model. The possible presence of any GUT 
superfield with vanishing vev will not alter the following discussion. 

It is useful to write down the superpotential in terms of the MSSM fields. The 
triplets H 3u ~ (3, 1, — 1/3), H 3d ~ (3,1,1/3), together with the two light Higgs 
doublets H 2u , H 2d stay in a pair of SU(5) pentaplets H, H. Omitting colour and 
weak isospin indices, in the interaction basis the terms of the superpotential which 
are relevant to proton decay read: 

w 3 Q T y u U c H 2u + Q T y d D c H 2d + L T y e E c H 2d 

+ Q T AQH 3u + U cT BE c H 3u + Q T CLH 3d + U cT DD c H 3d (3.15) 

+ "mTH 3u H 3d , 

where Q, L, U c , D c and E c denote as usual the chiral multiplets associated to 
the three fermion generations, y u ,d,e, A,B,... are the GUT scale Yukawa matrices of 
doublets and triplets and rriT is the triplets mass. 



By integrating out the colour triplets |78], ^, |8T| one obtains the following effec- 
tive superpotential 



1 

w eff = 



QAQ QCL + U C BE C U C DD C 



+ ... , (3.16) 



dots standing for terms that do not violate baryon or lepton number. The former 
(latter) term, which involves four fields with only left (right) chirality, is referred to 
as LLLL (RRRR) operator. w e ff generates dimension 5 operators with two fermions 
and two sfermions in the Lagrangian. 

At lower scales these dimension five operators give rise to the four-fermion op- 
erators relevant to proton decay, via a "dressing" process, that is through 1-loop 
diagrams with exchange of charginos and neutralinos. It turns out that the relevant 
contributions arise from chargino exchange J73L I7DI, B3I. When considering the op- 



erator QQQL, the most important dressing comes from wino exchange. There are 
four different terms in the Lagrangian, here called Of, (i = 1, ...,4), coming from 
the exchange of a wino 0. Charged higgsino exchange provides instead the most 



important dressing of U C E C U C D C [85]. Also in this case four Lagrangian terms arise, 



2 Gluino dressing contributions cancel among each other in case of degeneracy between first two 
generations of squarks [g3|, Q . 
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here called Of, (i = 1,...,4). The conventions adopted, the definition of the loop 
function / and all these operators are explicitly given in appendix [A] and in Table 
POL 



3.3.2 Hunting Dominant Operators 

The situation would then appear quite complicated. However, in the case of minimal 
577(5), where 

-2A = B = y u = y T u , -C = D = y d = yj , (3.17) 

it turns out that, in each category, only one term gives the relevant contribution to 
proton decay amplitudes and that the contributions from the other three terms are 
suppressed by at least two powers of the Cabibbo angle, A 2 . This is expected to be 
approximately true also in those models where the matrices A,B,... do not differ 
significantly from those of the minimal case. 

From eqs. ( |3.17| ), also the following relations follow 



Rl = L u K ,L e = R\ ,R e = L\ , (3.18) 

where R u ,d, e , L u ^, e are defined in Appendix [A] and K is a diagonal matrix of phases 
|86]. |8l) with determinant one. The two degree of freedom it contains are CP violating 
phases, in addition to the CKM one. Since proton decay amplitudes are insensitive 
to the overall phase, it is not restrictive to choose K to be of the following form 

K = diag(e i4 '\e i<t '\l) . (3.19) 

In the mass eigenstates basis and making use of the previous definitions, the 
four-fermion operators of Table |A.1| get greatly simplified. In fact, it turns out that 



the combinations in square brackets can be expressed just in terms of the diagonal 
Yukawa couplings of doublets and in terms of the CKM matrix, as shown in Table 



3T2| . This means that no trace is left of right handed mixings of fermions. 

The relative importance of the various operators can be approximately under- 
stood expressing them in terms of powers of the Cabibbo angle 9c = A. If the masses 
of the sparticles in the loop are such that they allow / (defined in appendix |A]) to 
vary in a small range, so that it can be quite safely factored out, and remembering 
the approximate GUT scale relations 

m u : m c : m t = A 8 : A 4 : 1 

rrid : m s : m& = A 4 : A 2 : 1 (3.20) 
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[L^A + A^L^LlCL^ - 


[V^K+V^V^fX, 


= r L 
— c l ijkl 


ok 


[LKA + A^L^ILJCL^ - 


[yi d) K^V]^ yd X 


= r L 

— ^2 ijkl 


ok 


[L T u (A + A T )L d ] l3 [L T u CL e ] kl - 


[V^K+V^yf^ 


— c 3ijkl 


ok 


[Ll{A + A T )L u } l0 [L T d CL e } kl - 


[y^K^^yiX 


= c L 


of 


[L + d ylBytLl] kl [RlDR + d ] i3 - 


+ [V + {y ( ffVyfUKV*y ( d \ 


— °1 klij 


o? 


[KBytLlUL+ylDRJ-]^ - 


- [yi d WX[v + yi d) KV*y d % 


— c 2 klij 


of 


[RlBR+] kl [L + d ylbyiL* u } i3 - 


- [y^v\ k i[v + y^KV\ y ffv% 


— n R 

— c 3 klij 


of 


[L+ylBRtUKDy+Ll^ - 


- [V + (yi d) ) 2 V] kl [KV*(y d d) rV% 


— n R 

— c 4 klij 



Table 3.2: In the case of minimal SU(5), the combinations in square brackets of 
Table |A.1| can be expressed in terms of the diagonal Yukawa couplings of doublets 
and in terms of the CKM matrix V = Vckm- See appendix |A] for the conventions 
adopted here. 



Vckm 



1 A A 3 
A 1 A 2 
A 3 A 2 1 



(3.21J 



one immediately obtains Table 3.3 



One can see at glance that Of, Of are, within their groups, the dominant op- 
erators. The contributions to proton decay amplitudes from other operators are 
suppressed by at least A 2 with respect to them, and can be safely neglected in the 
present estimate. The calculation then proceeds as summarized in Appendix |A[ 

Generally, the dominant channel is p — > K + v. However, also p — > ir + u is im- 
portant. In fact, the difference of a factor A ~ 0.25 doesn't produce a factor 1/16 
in rates. The ir + channel is favoured by phase space for a factor of 2, as appears 
from Table IA.2L Another factor 4 comes from Wick theorem. It results that 7r + is 
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ci 


c 2 


„R 


C 2 


v- 


^K+v 


A 10 tan/3 


A 12 tan/? 


A 10 tan 2 /? 

(Pr) 


A 22 tan 2 /? 


p 


— > TT + U 


A n tan/3 


A 14 tan/3 


A 11 tan 2 /3 

(fir) 


A 25 tan 2 /? 

(*W) 




c 3 


4 


4 


~J2 

c 4 


v- 


^K°l+ 


A 14 tan/3 
00 


A 12 tan/? 


A 17 tan 2 /3 
(e+) 


A 14 tan 2 /? 
(/^ + ) 


p 


^7r°/+ 


A 13 tan/3 
(m + ) 


A 14 tan/? 
(e+) 


A 18 tan 2 p 
(e+) 


A 15 tan 2 /? 
(e+) 



Table 3.3: Relative importance of the combinations in square brackets appearing in 
Table |0|, expressed in terms of powers of the Cabibbo angle 9c = A. Note that we 
summed over neutrino flavours, which are experimentally indistinguishable. 



suppressed by only a factor of 1/2 in rate, with respect to K + . As will be explained 
later, under certain circumstances linked to the two CP violating phases, the channel 
with 7r + could even be the dominant one. 



As pointed out in [35|, Of, Of can be comparable. For large tan/? and/or m~ h ^> 
rriu,, Oi dominates. But large values of tan/5 are discouraging because they lead 
to faster proton decay, exacerbating a possible contrast with experiment. For small 
tan/3 and m~ h < m^, Of gives the dominant contribution. 
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3.3.3 Proton Decay from Non-Renormalizable Operators 



It has been pointed out |87| that also non-renormalizable operators can lead to 
proton decay. If one views GUTs as effective rather than fundamental theories, 
then one has to consider also the role of non-renormalizable operators. The following 
operators, where it is understood that we are taking light generations, are completely 
allowed by gauge symmetry, 

QQQL U C E C U C D C 
Weff 9 c L — - — + c R (3.22) 

and lead, after the dressing process already described, to proton decay. Consider 
for instance the first term. It will be dressed as the LLLL operator arising from 
triplet exchange, so that, evaluating the ratio between its proton decay amplitude 
A nr and the amplitude from heavy triplet exchange A3, loop functions cancel each 
other. The difference between the first terms of ( 3.16Q and ( |3.22|) resides only in 



their multiplicative coefficients and in the mass at denominator. Approximately, for 
the dominant channel p — > K + u, with A ~ Mpi, one finds 

A 3 ~ M„A 10 tan/? l ' ' 

For small tan/3 and my ~ 2 • 10 16 GeV we have something like A nr /A 3 ~ c^A^ 7 , 
which means a reduction of a factor 10 8_; ~ 9 in lifetime for cl = 0(1). With the former 
values of tan/3 and my, as we will see, the expectation for proton lifetime is roughly 
10 31 yr, which is lower then actual bounds. For non renormalizable operators one 
should then expect a lifetime of around 10 22 yr. To restore agreement with the 
bounds, one has to impose c^ to be at least smaller than A 8 ~ lO - ^ 6 ). This is 
indeed a very large fine-tuning. The most natural solution is to impose flavour 
symmetries, which could give a valid reason to justify the required smallness of cl 
thus adequately suppressing these disturbing non renormalizable operators. Flavour 
symmetries are usually largely employed in model building, because they provide a 
clue to understand the observed picture of fermion masses and mixings. So, it would 
be interesting to study if, with the tool represented by flavour symmetries, one can 
solve both problems. 



3.4 How to Prolong Proton Lifetime 

There are several parameters entering the formula for proton lifetime in minimal 
SU(5), so that one can should be careful not to underestimate theoretical errors. 
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For instance, by applying the chiral Lagrangian technique, one obtains the for- 
mulas for the calculation of proton decay rates which are given in Table |A.2|. One 



can see that the hadronic matrix elements a = —f3 are very important ingredients, 
because proton rate scale as their second power. Unfortunately, these quantities are 
not very well known. Commonly, their quoted range |88[ is 0.03 -j- 0.003 GeV 3 and 



the smallest value used in the calculations. However, it is worth to stress that with 
the largest value, the rate would be increased by roughly two orders of magnitude. 
However, one could also construct elaborated models where one or more of the 
important parameters change with respect to the minimal model. So, it is of crucial 
importance to analyze how much the prediction for proton lifetime is affected if one 
moves a certain parameter in its possible range. In what follows we will briefly give a 
survey of the most interesting circumstances in which proton lifetime can be raised. 
Finally some numerical example will be provided. These are to be compared with 
those given in the next chapter, where a concrete model will be examined. 

3.4.1 Altering Triplets Yukawa Couplings 

Clearly, the previous discussion about dominant operators, apply only to those mod- 
els where the matrices A,B,... don't differ significantly from those of the minimal 
case. If, for some reason, such matrices are much more suppressed than in the 
minimal case, then proton decay rates would be reduced. This possibility has been 



explored |89| in the literature and, in general, the suppression of certain matrix 
elements follows from the group-theoretical construction of the model. One can 
also investigate if the same strategy employed to suppress the Yukawa couplings of 
triplets could also be effective in differentiating charged leptons and squarks spec- 
tra. Even if this is certainly an attractive possibility, it seems not very effective in 
increasing proton lifetime. Indeed, when implemented in a realistic framework, it is 
difficult, by only this mean, to adequately suppress all channels at the same time. 

3.4.2 Raising the Triplet Mass 

A really crucial parameter is the mass m^. The dependence is strong, t p ~ m^ 4 , so 
that a large value of m^ would drastically decrease the rates for all channels. As will 
be shown in the following, when considering the effect of heavy thresholds on the 
running of the gauge couplings, it follows that in minimal SU(5) and with SUSY 
particles lighter then 1 TeV, ttlt should not exceed 10 16 GeV. This value causes 
proton decaying at a too fast rate, in conflict with experiment. Minimal SUSY 
SU(5) is then in serious troubles. Interestingly, a large value for my is instead 
preferred in the MDM. Notice that in this case rrij- is not a physical quantity but 
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an effective one, given by eq. ( |3.12|) , and could in principle even exceed Mpi. In the 



next chapter we will discuss a realistic version of the MDM which maintains this 
attractive feature. On the contrary, in those realizations of the MDM and MVM 
where a large number of heavy representations lead to huge threshold corrections 
(often of both signs), the correlation between the value of rriT and the prediction for 
a 3 (mz) is lost. 

3.4.3 Moving the Two CP Violating Phases 



There are phases entering the proton decay amplitudes [^fj] which don't affect any 
of the low energy observables, like for instance fermion masses and mixings [], but 
which play a really interesting role, because they can produce sizeable effects on 
proton rates. In fact, proton lifetime can easily span over one order of magnitude, 
just allowing the phases to vary in their range. 

As usual, it is instructive to see this firstly in minimal SU(5), and then comment 
on its extensions. 

The phases of minimal SU(5) 

In minimal SU(5) the two phases can be clearly identified. Following the definition 
of eqs. ( [3.18| ) and ( p,19|) , after a bit of work, one can see how they enter in the decay 
amplitude. In particular, considering the dominant operators Of, 0±, one realizes 
that the dominant amplitudes are proportional to the following quantities 

A{K + u^) oc (e i<fc af + af) 

A{K + u T ) oc (e^ 2 af + af + e^af) , (3.24) 

where af , a~ and af; are complex numbers which contain the piece of amplitude 
coming respectively from c and t for wino dressing and from t c for higgsino dress- 
ing, af , af are comparable, so that, depending on </> 2 , they can interfere construc- 
tively or destructively in the amplitude. One can have three possibilities. Firstly, if 
the contribution from higgsino exchange which enters in A{K + v T ) is smaller than 
the contribution from wino, a global suppression of the channel p — > K + u can be 
achieved [S^| for certain values of 02- Note that, in this situation, the role of <pi 



is irrelevant. Secondly, if also af c is of the same order of a~,a~, also 0i, one can 



never achieve at the same time a suppression of both A(K + u fJ- ) and A{K + v T ). Sum- 
ming their corresponding rates, there isn't a sizeable variation of the total rate for 



3 In the diagonalization of the mass matrices, these phases can always be absorbed through 
suitable fields redefinitions, so that only one phase is left, the Vckm one. 
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p — > K + v over the whole region of the space spanned by (0i,02)- The third and 
last possibility is that a? dominates over the other contributions, in which case any 
effect of interference is totally lost. However this situation would require very large 
values of tan (3, leading to unacceptably fast proton decay. 

In the first case we obtained a suppression of the main decay channel p — > K + v. 
What about the channel p — * tt + u, for which the rate is in general expected to 
be just a factor of 2 smaller? In the first case, for those <p2 values which produce 
destructive interference in the K + v rate, it becomes the dominant channel. In the 



second case, instead, it is always subdominant |85 

The phases in generic models 

When constructing realistic SU(5) models which break the relation y^ = yj, then, 
in general, A,B, ... are no more linked to fermion masses. However, it is reasonable 
to expect their entries to be of the same order of magnitude, that is: 

-2A~B~y u , -6~D~y d ~yT . (3.25) 

In this situation one again finds the same results of the minimal case. In fact, the 
dominant amplitudes can be still written as 

A(K + ^) oc (b^ + bf) 

A(K + u T ) oc (b± + bf + b? c ) , (3.26) 

where b^,bj: and bf c are complex numbers which contain the piece of the ampli- 
tude coming respectively from c and t for wino dressing and from t c for higgsino 
dressing. Thus, all the comments made when discussing eq. (|3.24|) still apply. For 
instance, when bf c is not important, interference occurs depending on the relative 
phase between b~,b~;. 

3.4.4 Playing with SUSY Particles 

In the proton decay amplitude the loop function /, defined in eq. (|A.3|) , plays an 
important role. It is obtained by dressing the dimension 5 operators in order to 
convert them into four-fermions operators. This function depends on the masses of 
SUSY particles which run in the loop. To have an idea of the range of values in 
which / can lie, it is not a too bad approximation to assume that all squarks and 
sleptons are degenerate, with common mass m. In this case, the loop function for 
the exchange of the chargino c of mass m^ greatly simplifies 

,2" 

(3.27) 



f /; ' 

Jc 



167T 2 (m 2 — m~) 



,2 2 2 



m z — mi — mi log — ^ 
' m\ 
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The behaviour of / is the following: by fixing m, it grows until 2m £ and then 
decreases; by fixing 777.5, it grows when m decreases. In the limit m^ <C m, then 



/cOC 



rrir 



171' 



(3.28) 



This means that if sfermions are much more heavy than charginos, proton decay 
rates get greatly suppressed. However, it is not possible to push m towards values 
greater then 1 TeV, if one wants SUSY to be the solution of the hierarchy problem. 
The hypothesis of degeneracy (or nearly degeneracy) of sfermions is useful be- 
cause the various / can be factorized, that is, summing over v flavours, one can 
symbolically write: 

T(p^K + u)^\Ua L + f- h a R \ 2 (3.29) 

where a^R represent the contribution to amplitude coming respectively from wino 
and higgsino dressing. If ma, = m~ h the cll contribution is the dominant one. Dom- 
inance of dR occurs only if rriw <C m~ h and/or tan/5 is large. To quantify the un- 
certainty on r p induced by the SUSY spectrum, one has to study the behaviour of 
/ for reasonable values of m^ ^ and sfermion masses. Considering typical SUGRA 
derived spectra with values of the soft breaking parameters in the range from 100 
GeV to 1 TeV, one finds a variation of about one order of magnitude in rates. 



3.4.5 Numerical Examples 

The previous observations are clarified by the following example, where we com- 
pare two different set of textures for fermion masses which lead to exactly equal 
phenomenology, that is they both give the same masses and the same Vqkm- 
The first set is the well known Georgi-Jarlskog proposal |KJ 








a\ 6 





Vu J = 


a\ 6 





b\ 2 







bX 2 


1 



yf 






g A 3 0] 


g'X 3 


j\ 2 





1 



(3.30) 



For tan/? = 3 and with A = 0.23, the values a = 0.799, b = 0.813 and g = 0.005 + 
0.372i, g' = 0.004-0.2852, j = 0.387 give: {m U} m C} m t ) = (0.0012,0.28,154) GeV, 
(m d ,m s ,m b ) = (0.0009,0.24,1.17) GeV, \V ud \ = 0.22, \V ub \ = 0.027, \V cb \ = 0.043 
and the Jarlskog invariant Iq = 2.2 • 10~ 5 . These GUT scale values, when run to 
low energies, are compatible with the measured fermion masses and mixings. The 



4 In our convention Dirac mass terms are given by LmR. 
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second set is obtained by replacing y^ J with the following matrix possessing a large 
right-handed mixing in the 23 plane: 



vf° 



o 

A* A 3 




A\ 3 

1 




B\ 2 

1 



(3.31) 



This particular form for the down quark mass matrix has been studied in detail 
by Hagiwara and Okamura |JT 



For tan/? = 3 and A = 0.274 + 0.4752, B = 
—0.217 — 0.711z, one obtains exactly the previously quoted values of masses and 
mixings. These two set, due to the different values of CP violating phases, give 
different predictions for proton lifetime. Keeping my = 2-10 16 GeV, a = —f3 = 0.003 



m„- 



GeV 3 , 



with m~f 



= 250 GeV, m~ h = 125 GeV and assuming sfermions to be degenerate 
837 GeV, then one predicts the values shown in Table 





r p /B(p - K+v) 
10 31 yr 


T p /B(p -> Tl + V) 

1Q 31 yr 


1 st set (GJ) 
2 nd set (HO) 


5.6 
0.7 


2.8 
1.2 



Table 3.4: Comparison of the predictions for proton lifetime between the two set of 
mass matrices. 

The values of the various parameters adopted here are suitable for the case of 
minimal SU(5). As pointed out, they suffer from large theoretical uncertainties. 
However one can recognize that minimal SU(5) is in serious troubles, predicting a 
mean value for proton lifetime which is smaller then the present bound by two or 
three orders of magnitude. 

In addition, one can also conclude that, even in the best constrained case of 
minimal SU(5), results depend on the textures utilized, even if they give completely 
equivalent low energy phenomenology. The predictions differ for about one order 
of magnitude and also the main channel is different. This discrepancy has to be 
attributed entirely to the two CP violating phases. In fact, it is possible to recognize 
that the first case is an example of destructive interference for the amplitude of 
K + v (see eq. ( [3.24J )). The dominant mode is then tt + u. In the second case, where 
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no destructive interference occurs, proton lifetime is smaller by about one order 
of magnitude and the mode tt + v is slightly subdominant. The result is that an 
uncertainty of about one order of magnitude arises just from the ignorance of the 
input textures. 

One can wonder about the effect on rates induced by taking appropriate textures 
for y e , and not just yj, as required by minimal SU(5). As will be discussed in the 
next section, one simple and well known possibility is to assume that certain elements 
of these matrices arise from a 45 of Higgses and not from a 5. In the former case, 
the element of y e will be multiplied by a factor —3 with respect to the corresponding 
one of yj. The matrices C, D, arising from the interaction of the 5 with matter, will 
have zero in the corresponding entry For the first set, Georgi and Jarlskog pointed 
out that realistic Yukawa couplings for charged leptons can be obtained if the (22) 



element comes entirely from a 45, so that m e = 1/3 m^m^ = 3m s ,m T 
the second set, these relations can be achieved if yd 2 3 comes from a 45: 



mb. For 



,GJ 



g'X 3 

g\ 3 -3j'A 2 

1 



,HO 



A*X 3 

AX 3 1 

-3B\ 2 1 



(3.32) 



Correspondingly, C^ J = D^ = and C^f = D^f = 0. It easy to adapt the calcu- 
lation of proton decay lifetime to this particular case. Under the same assumptions 
made to obtain Table 13.41 one now obtains Table |3~5. 





r p /B(p -► K+v) 
10 31 yr 


T p /B[p -> 7T+Z/) 

10 31 yr 


1 st set (GJ) 
2 nd set (HO) 


5.8 
0.8 


2.9 

1.5 



Table 3.5: Comparison of the predictions for proton lifetime between the two set of 
mass matrices in the realistic case. 



The results are very similar to those of the minimal case. Indeed this is a simple 
and very particular example, however, it is true also in general that the use of simple 
realistic textures also for charged leptons doesn't have very sizeable effects on proton 
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lifetime. Achieving a certain degree of suppression through a differentiation between 
triplet and doublet Yukawa couplings is a technically realizable possibility only for 
much more complicated models, where a number of ad hoc assumptions must be 



made |89| . It is clear that one is interested in looking for such mechanisms that lead 
to construct non minimal models where proton lifetime can be pushed comfortably 
above experimental bounds. As pointed out, a resolutive way to do this, would be 
to have triplets much more heavy that in the minimal case. In the next chapter we 
will discuss one of these models. 

3.5 Fermion Masses and Mixings 

The aesthetically attractive fact of putting in the same representation conceptually 
different objects like quarks and leptons has mainly two consequences. The most 
apparent one is the fact that quarks and leptons mass matrices turn out to be strictly 
correlated. Despite the difficulties found, this can give interesting starting points 
to face the flavour problem. The second consequence is the possibility of violating 
baryon and lepton numbers. As shown, this leads to the striking prediction of proton 
lifetime, with values in the region already partially excluded by present bounds. 
Minimal models fail in correctly predicting both mass relations and proton lifetime. 
In the following we will review the usual methods which can be adopted to 
share some light on the resolution of the first problem. It is clear that finally one 
is interested in understanding if these mechanisms are compatible or even help in 
the solution of the other problems like proton decay and the other typical GUT 
problems. 

3.5.1 Mass Relations 

The evident problem of minimal models is the prediction of wrong mass relation for 
down quarks and charged leptons at M GUT . Assuming an exact i?-parity discrete 
symmetry under which matter fields are odd and H, H are even, if the superpotential 
contains the following Yukawa interaction terms 



WYuk . 



-M* G u ^d H E t abode + V2 K B G d ^5A H B (3.33) 



where A,B,... = 1,...,5, H, H are the pentaplets which contain the pair of light 
doublets and G u , Gd are 3x3 matrices in flavour space. Writing Wy u k in terms of 
MSSM multiplets, the interactions with doublets read: 

w Yuk 3Q T G u U c H 2u + Q T G d D c H 2d + L T G e E c H 2d . (3.34) 
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Then we get the following Dirac mass matrices at Mqut'- 

m u = G u -j= , m T e =m d = G d —= , (3.35) 

where v u , v d parameterize the vevs of H, H. The sharp transposition relation 
mf = m d is in conflict with what is expected on phenomenological grounds, by 
evolving the values of down quarks and charged lepton masses from their values 
measured at the electroweak scale to the GU scale, i.e. the Georgi-Jarlskog relations: 

1 
m e = ^rn d , m M = 3m s , m r = m h . (3.36) 

This problem is even exacerbated in minimal 5*0(10), where rh u = m^ = rh V D oc 
m d = mj = rh e . However, at least in SO (10) there is naturally the place for 
right-handed neutrinos. 

On the contrary, in SU(5) the inclusion of neutrino masses represents an ad- 
ditional problem. A very simple way out is the addition of a gauge singlet field 
\&i, which contains right-handed neutrinos. If one views 577(5) as descending from 
50(10), this is of course not only natural but also necessary. Then, in the superpo- 
tential these terms are allowed 

*! G„ V- 5A H A -^MV X G M tf i (3.37) 

which realize the see-saw mechanism. If M is just a bit smaller than Mqut, then 
the atmospheric neutrino mass scale is naturally obtained. However one should em- 
phasize that, because the term which gives mass to right-handed neutrinos doesn't 
break 577(5), then M should in principle be expected of order of the cut-off of 
the theory, that could even be as large as M Pi . The advantage of groups like 
50(10), SU(A) <g> SU(2) <g> SU(2) is that they directly possess a particle with the 
same quantum numbers of right-handed neutrinos, whose mass scale thus results for 
dynamical reasons of order of the GUT scale. 

3.5.2 How to Achieve Realistic Spectra 

The flavour problem consists in explaining the origin of the difference between the 
three fermionic families, difference which causes their masses and mixings to be 
organized in strongly and bizarre hierarchical structures. Probably this problem is 
linked to the more fundamental one: why just three families? It is clear that the 
idea of GU cannot alone answer these deep questions. The only thing it can offer 
is a much more constrained framework with respect to the MSSM, where one can 
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look for the criterion according to which the different mass matrices are linked one 
to each other. One can thus learn a lot about what would be allowed to be present 
at scales higher then the electroweak one and what could not at all be there. 

To understand the situation it is convenient to take SU(5) as reference. This is 
not at all restrictive because 5*0(10) can be always thought in terms of its SU(5) 
content. The problems one has to face are two: firstly a mechanism to produce 
a progressive suppression in the masses of the three generation, secondly a way 
to reproduce the correct values of masses and mixings. The usual road taken to 
overcome the first problem is by the use of flavour symmetries [|51[], for instance 
a continuous U(l). This strategy is very successful, even if one is left with the 
discomfort of having done ad hoc assumptions for the charges. Apart from this, it 
turns out that it is very difficult to realize the correct spectra without tuning some 
parameter, even in the more complicated extensions of minimal GUT models. 

Possible ways for accounting for neutrino masses and mixings have been already 
discussed. Here we will review the main strategies one can use to break the bad trans- 
position relation between down quarks and charged leptons. Since 10x5 = 5 + 45, 
a mass term for quarks and leptons can be realized only with a 5 and a 45. The 
5, which is of course nothing that the pentaplet containing one of the light Higgs 
doublets, predicts the sharp transposition relation. So, to have a realistic spectrum 
there is the absolute need for a 45, which succeeds in breaking the unwanted trans- 
position. The 45 can be a "true" or an effective one. The first scenario, with a true 
45, was proposed a long time ago. Now it seems not very attractive because the 
true 45 contains a true doublet, which needs to be massive. One can then introduce 
another 45 to construct a direct mass term (and also to canceling anomalies). But 
it is well known that the two additional doublets cannot be light. If this were the 
case then the prediction for « 3 (m^) assuming coupling unification would be given 

by: 

56 — 25uh 

a 3 {m z ) = <Xem{m z ) — . r — — r 3.38 

(120 + Qbn H )sin z 6(m z ) — (24 + 3dn H ) 

where Snn is the number of light extra doublets, in addition to the two of the M SSM. 
Taking a em (mz) = 1/128 and sin 2 8(mz) = 0.23, one has: 

a 3 (m z ) = 0.12 5n H = (MSSM) (3.39) 

a 3 (m z ) = 1.13 5n H = 2 . (3.40) 

It is then absolutely unavoidable: to give a GUT scale mass to the four triplets; to 
maintain at the electroweak scale two linear combinations of the doublets in 5, 45 and 
5,45 but raising at the GUT scale their two orthogonal combinations. Here come 
the difficulties, because this asymmetric splitting between doublets and triplets is a 
much more complicated version of the usual doublet triplet splitting problem. 
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To escape from these difficulties one can make complete use of the possibility 
of working with non renormalizable operators and explore the field offered by the 
effective 45. The simplest possible mechanism is the well known proposal made by 
Ellis and Gaillard |57| of considering the 45 of the product 5 x 24 = 5 + 45 + 70. 
This is indeed a very economical solution, because the 5 is already there and the 
24, which must have non vanishing vev, can be naturally the one responsible for 
the breaking of 577(5). Thus, this mechanism uses just the field content of minimal 
SU(5). The following operator 



A 



V2 9tf t G dij *- 5Cj H B { -^1 (3.41) 



where A is the cut-off, gives Gdij — Udij = —^/tyeji- Note that both channels 5 
and 45 of the product 5 x 24 contribute to the mass. Also operators with more 
insertions of S can be useful in realistic examples. The number of the invariants 
obtained by contracting the indices in all possible ways grows with the number of 
insertions made, so that one can use the freedom in choosing their relative coefficient 
to fit the masses. This is however a kind of proof that realistic mass spectra are 
realizable rather than a genuine prediction of the model. 

There are other possibilities, in addition to the use of a 24. In the MDM, 
for instance, the 75 responsible for the breaking of SU(5) can be successfully well 
employed, without the introduction of additional fields. In fact, since 5 x 75 = 
45 + 50 + 280, in this case only the pure 45 channel contributes. For instance 

V2 K B > G dij *scj H d ^LL (3.42) 

gives Gdij — Udij = ~^/^Veji- In the next chapter the strategy of inserting the 75 
to obtain a phenomenologically acceptable spectrum of fermions will be employed 
in the context of a realistic version of the MDM. 



3.6 Coupling Unification 

The most outstanding support to the idea of Grand Unification comes from the 



quantitative success of coupling unification in SUSY GUTs J60|. Assuming just the 
field content of the MSSM (with two light Higgses), in the leading order approxima- 
tion (LO), that is at one loop, the RGE for the running coupling measured at irtz 
read: 

1 +^l»m (3.43) 



aci(Q) oci{m z ) 2ir \m z 
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with 

fci = 

5 



S3 
6i = -=- , &2 = 1 , &3 = -3 . (3.44) 



Nowadays it is possible, thanks to the quite precise measurements of the three 
running gauge coupling at mz, to single out a not too broad range for the scale of 



SUSY unification Mqut- Fig- 3.1 shows the flow of the gauge coupling in the one 



loop approximation starting at mz with the PDG data of last summer |92| , that 
is, in the WS scheme: a em (m z ) = 127.934 ± 0.027, sin 2 6(m z ) = 0.231 17 ± 0.00016, 
m z = 91.1872 ± 0.0021 GeV and the world average a 3 (m z ) = 0.1184 ± 0.0012. As 
can be clearly seen, unification is very successful and the associated scale Mqut 
turns out to be in the small range (1.9 -=-2.1) • 10 16 GeV, with the previously quoted 
errors. 

Now it is also well established that GUT without SUSY is neatly disfavoured, 
as appears clearly from the comparison of figs. |3.1| and |3.2| . For the SM equations 
analogous to flglfl) hold, but with b x = 41/10, b 2 = -19/6, 63 = -7. 



Going beyond LO turns out to be necessary because the corrections from two 
loop running on the prediction for oc^imz) are of the order of 10%. It is reasonable 
however to disregard three loop running, which would give correction of only about 
1% on a${mz)- In the following we will describe the calculations relative to the two 
loop prediction of a^{mz) starting from ax(mz), OL2{raz) and assuming coupling 
unification. We will give numerical results for the case of minimal SU(5). They will 
be useful for comparison with the analogous predictions obtained in the MDM case, 
whose analysis will be present in detail in the next chapter. 

Different GUT models give different predictions for otz{mz) because the contri- 
bution to its running coming from threshold effects is sizeable. If one is convinced 
that the success of coupling unification in the LO approximation is not just a coin- 
cidence, then, in model building, one is constrained to introduce as few as possible 
GUT fields, which in addition have to lie in not too large representations. It is clear 
that working with a lot of field is equivalent to have a lot of degrees of freedom 
at disposal, which can be arranged ad hoc to reproduce the data. The elegance of 
the model is completely lost, but at least one possess a sort of existence proof of 
the realizability of low energy data starting from a GUT scenario. Working with 
few fields, instead, it is much more hard to give a completely satisfactory picture 
of particle phenomenology. However, this latter point of view is the only one from 
which a true deeper understanding of fundamental constituents could come. 
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SUSY Standard Model 
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Figure 3.1: Evolution of the running gauge couplings within the MSSM. The curves 
are obtained in the 1 loop approximation by evolving the gauge couplings from their 
measured values at m z ||: ai(m z ) = 0.016945 ± 0.000007, a 2 (m z ) = 0.033813 ± 
0.00003 (from a em {m z ) = 127.934 ± 0.027, sin 2 fl(m z ) = 0.23117 ± 0.00016, m z = 
91.1872 ± 0.0021 GeV) and a 3 (m z ) = 0.1184 ± 0.0012. 
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Figure 3.2: The same as fig. [Obut within the SM. 



70 CHAPTER 3. PROBLEMS OF SUSY GUTS 

3.6.1 Coupling Unification Beyond the Leading Order Ap- 
proximation 

It is well known that, going beyond the LO approximation, in the minimal version 
of SUSY SU(5) the central value of a^mz) required by the constraint of coupling 
unification is somewhat large: a^{mz) ~ 0.13 |64|, [5(| [37], |100| , |101| . 



From the one-loop renormalization group evolution of the couplings in the MSSM 
with only two light Higgs doublets, eqs. (|3.43|) , by evolving two couplings from their 
values measured at mz and imposing the constraint of unification, one can run back 
to low energies with the third coupling thus giving a prediction for its value at mz- 
In this procedure it is preferable to take as input the two best measured couplings, 
that is ati(mz), Q-iimz)- In fact, defining in the MS scheme a = olqed^z)-, sin 2 9 = 
sin 2 0(mz), 

oci(m z ) = - -tb > «2K) = — -2T (3.45) 

1 — sin 9 sin 9 

and the predictions are: 

\n(M§l T /m z ) - 



% 



(0) 
5 

(0) 



vr(3 


- 8 sin 2 9) 


14a 

28a 


36 sin 2 9 - 
28a 


3 



(3.46) 
60 sin^ 9-12 y ' 



where a 3 = as(mz), a$ = oc^Mgut) and the superscript (0) refers to a quantity 
evaluated in the LO approximation. It is worth to stress that the LO results are the 
same in minimal SUSY SU(5) and in all its extensions where the MSSM (3 functions 
apply until M GUT . 

By using a" 1 = 127.934 and sin 2 9 = 0.231 one obtains af ] = 0.118 and M^ T = 
2.1 10 16 GeV. 

To go beyond the LO one must include two loop effects in the running of gauge 
couplings, threshold effects at the scale msusY, close to the electroweak scale, and 
threshold effects at the large scale Mqut- The unification scale Mqut is not uni- 
vocally defined beyond LO and here we choose to identify it with the mass of the 
SU(5) superheavy gauge bosons my. Following ref. 0] it is possible to include 
these effects by writing: 

1 1 bi , f M GUT \ Si 

— + -i- In —^ + - (3.47) 



aci(mz) a 5 2ir \ m z J vr 

where a$ = o/.§{Mgut) and the next to leading corrections are contained in Si, i — 
1, 2, 3 and can be separated according to their origin. The logarithmic contribution 
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coming from two loop running, <5j(2), and the contributions of SUSY and GUT 
thresholds, respectively 8{(l) and 8i{h): 

S i = S i (2)+S i (l) + S i (h) (3.48) 

The LO eqs. ( |3.46| ) are thus improved by: 



m z J \ m z 

1 1 

a (0) 
(1 + 4 ; 5 3 ) 

where the various <5x,5,3 are linear combinations of Sf. 

5x — — ($2 — Si) 

28 

8 5 = -4(3352-550 
zo 

<5 3 = ^-(5<*i-12* 2 + 7<J 3 ) (3.50) 

In order to obtain the prediction for a 3 it is necessary to specify the various 
contributions which enters 83 = 83(2) + 83(1) + 83(h). 

Two-loop running 

The one from two-loop running is finite and given by: 

^)-E|>»(^|4°'.,(^)) (3.51) 

where 

199/50 27/10 44/5 

/;,,= 9/10 25/2 12 

11/10 9/2 7 



J 'j 



(3.52) 



In particular 

<J 3 (2) = -0.823 . (3.53) 
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SUSY Thresholds 

The one loop thresholds from SUSY particles read: 



*,(,) = -X>?>0)>n(^) 



(3.54) 



where j runs over all light particles of mass rrij and b\ (j) can be read from Table 



(0, 



3.6 



sparticle 


b?U) 


4°(J) 


bl\j) 


gluinos 








1 


winos 





2/3 





higgsinos 


1/5 


1/3 





extra Higgses 


1/20 


1/12 





Ql 


1/60 


1/4 


1/6 


Ur 


2/15 





1/12 


d R 


1/30 





1/12 


T L 


1/20 


1/12 





e R 


1/10 









Table 3.6: SUSY thresholds. 
To give a concrete example, one has to specify a certain SUSY spectrum. We 



will consider as representative the one given in Table [B.l| , where the additional 
freedom related to the parameters mo, mi/2, ft, ttlh can be fixed by choosing 0.8mo = 
O.87711/2 = 2/i = mil and taking as a definition msusY = m>H, so that all particle 
masses can be expressed in term of msusY- 
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The result is then 



*s(Q 



19 / 
-0.503 + ^— In ( 



msusY \ 



287r \ rriz J 



(3.55) 



GUT Thresholds 

Threshold contributions from heavy particles can be expressed as: 



~ \m GUT J 



(3.56) 



with C\ = 0, C2 = 1/6, C3 = 1/4 and the sum is over all heavy particles. 

To specialize to the case of minimal 577(5) one has to consider its heavy spectrum. 
After the Higgs mechanism, 12 complex fields of the 24 representation (the 24 would- 
be Goldstone bosons) are eaten by the 24y of gauge vector bosons, which thus 
become massive. The other 12 complex fields in the 24 acquire instead a heavy 
mass. In addition to them, there are the colour triplets of the 5 and 5 representation. 
Thresholds can be evaluated according to: 



state (SU(3),SU{2),U(1)) 


b[ h \j) 


bf\j) 


bf\j) 


mass 


(8,1,0) 
(1,3,0) 

(1,1,0) 








1 



3/2 




0.2m s 


(3, 1,-1/3), (3, 1,1/3) 


1/5 





1/2 


ttlt 



Table 3.7: GUT thresholds, mj = 5M 24 , with M 2 a defined as in section |3.1.1| . 
For this case one obtains 



5™ n {h) 



-0.0114 



9 



14tt 



In 



rriT 



M, 



GUT 



+ ^ln / "^ 



147T 



M, 



GUT 



(3.57) 
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Results 



The first of eqs. ( p. 49 ) turns out to be very useful. The contributions to its correc- 



tive term 5x coming from two-loop running and SUSY thresholds induced by the 
spectrum of Table |B.1| are: 



S x {2) = 0.469 (3.58) 

Sx(l) = -0.036 -gl»(^f) • (3.59) 



Thresholds from GUT particles within minimal SU(5) give: 



6x(h) = 0.060 - A In (J?*-) + 1 In { ^- ) . (3.60) 



l GUT J J-^ \ 1V1 GUT 



Here in minimal SU(5) it turns out clearly, but it is true in general that, keeping 
fixed m S usY, the first of eqs. Q3.49] ) can be considered as a relation between the 



GUT masses. In this case it relates my = Mqut to m^^mx- 

Substituting in 63(h) the expression thus obtained for Mqut, one can always 
write 63 according to the following expression: 

*3 = k + - In ( ^*) - A In (™) (3.61) 

2tt V rn z J 5tt \M® t ) l ' 

This expression for <5 3 is general because we have subdivided the various contri- 
butions in a numerical coefficient plus logarithmic contributions. The logarithmic 
term with rrisusY comes from particles with masses near msusY- Similarly, the 
logarithmic term with rnj- arises from particles with masses of order m? ~ Mqut- 
The definition of ttit is the mass of Higgs colour triplets in the minimal model or 
the effective mass that, in realistic models, plays the same role for proton decay 



64] , |96| . Notice that that there isn't any dependence on m^. The term indicated 
with k = k(2) + k(SUSY) + k(M GUT ) contains the contribution of two loop diagrams 
to the running couplings, fc(2), the threshold contribution from states near ttisusy, 
k(SUSY), and the threshold contribution from states near M GUT , k(M GUT ). This 
way of subdividing the contribution is useful because the threshold contributions 
would vanish if all states had the mass msusY or M G ut so that only mass splittings 
contribute to k(SUSY) and k{M G ur)- 

The values of k(2) and of k(SUSY) are the same in minimal SU(5) and in 
its realistic extensions, because only the heavy spectrum changes. Typical values 
are k{2) = -0.733 and k(SUSY) = -0.510. The value for k(SUSY) = -0.510 
corresponds to the representative spectrum displayed in Table [B.l| . Note that the 
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Figure 3.3: Contours of a 3 = a s {m z ) in the plane (tususy, ™>t), for minimal SU(5). 
The SUSY spectrum is parametrized as in Table [B.l| . 
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value of k(McuT) is practically zero for the 24 of the minimal model. The reason 
is that in minimal SU(5) all coloured GUT multiplets are degenerate, as appears 



clearly giving a look to Table 3.7. The result is thus: 



k = -0.733 - 0.510 = -1.243 minimal model (3.62) 

In fig. |0] the prediction for a 3 in the plane (rn SUSY ,m T ) is shown for the case 
of minimal SU(5). It can be seen that, in this case, we need to take msusY as 
large as possible and m-r as small as possible. But the smaller is mr, the faster 
is proton decay. The best compromise is something like msusY ~ 1 TeV and 
rriT ~ Mq(j T , which leads to a^rrtz) ~ 0.13 which is still rather large and proton 
decay is dangerously fast, in conflict with experimental limits. In fact, with this 
value for %, we already know from the discussion of section |3.4.5| that proton 



lifetime for the channel p — > K + i> is expected to be around 10 31 yr and that, in the 
minimal model, there is no way to push it above the experimental limit of ~ 2 • 10 33 
yr. Thus, we are lead to conclude that minimal SUSY SU(5) is ruled out. 

As already pointed out, models based on 50(10) are highly model dependent, 
because the choice of Higgs fields which are introduced to realize the breaking of 
5O(10) affects many predictions, among which proton decay and coupling unifica- 
tion. For instance, in the quite economic 5*0(10) model of ref. ||42| , the corrections 
to a 3 (m z ) coming from representations like 10, 10, 16, 16 and 45 are expected to be 
positive and of order 10% with respect to the experimental value of a 3 (mz), so that 
one needs to invoke negative corrections of the same order arising from Planck scale 
effects or other thresholds. 



Chapter 4 

A realistic SUSY SU(5) GUT 
Model 



It has emerged from the survey of the previous chapter that, despite the idea of 
GUT is really very attractive, its actual realization is not at all precisely defined. 
Consider also that at the Planck scale, Mpi, the unification of gravity with gauge 
interactions is a general property of superstring theories. But in simple GUT models 
gauge interactions are unified at a distinctly lower mass scale Mqvt within the 
context of a renormalizable gauge theory. However, this neat separation between 
gauge unification and merging with gravity is not at all granted. The gap between 
Mqut and Mpi could be filled up by a number of threshold effects and several layers 
of additional states. Also, coupling unification could be realized without gauge 
unification, as suggested in some versions of superstring theory or in flipped SU(5). 
The model we aim at should not rely on large fine tunings and must lead to an 
acceptable phenomenology. This includes coupling unification with an acceptable 
value of a s (mz), given a and sin 2 9w at mz, compatibility with the more and more 
stringent bounds on proton decay |61j, |7?fl , agreement with the observed fermion 



mass spectrum, also considering neutrino masses and mixings and so on. 

Prompted by recent neutrino oscillation data some new studies on realistic GUT 



models have appeared in the context of SO (10) or larger groups ||2], [93|]. Here we 
address the question whether the smallest SUSY SU(5) symmetry group can still 
be considered as a basis for a realistic GUT model. We indeed present an explicit 



example of a realistic SU(5) model [94], which uses a (7(1) flavour symmetry as a 
crucial ingredient. In principle the flavour symmetry could be either global or local. 
We tentatively assume here that the flavour symmetry is global. This is more in 
the spirit of GUT's in the sense that all gauge symmetries are unified. The associated 
Goldstone boson receives a mass from the anomaly. Such a pseudo-Goldstone boson 
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can be phenomenologically acceptable in view of the existing limits on axion-like 
particles |95| . 



In this model the doublet-triplet splitting problem is solved by the MDM pB 



stabilized by the flavour symmetry against the occurrence of doublet mass lifting 
due to non renormalizable operators. Relatively large representations (50, 50, 75) 
have to be introduced for this purpose. A good effect of this proliferation of states 
is that the value of a s (mz) obtained from coupling unification in the next to the 
leading order perturbative approximation receives important negative corrections 
from threshold effects near the GUT scale. As a result, the central value changes 
from a s (mz) ~ 0.130 in minimal SUSY SU(5) down to a s (mz) ~ 0.116, in better 



agreement with observation []64], |96|, J97fl . The same U(l) flavour symmetry that 
stabilizes the missing partner mechanism is used to explain the hierarchical structure 
of fermion masses. In the neutrino sector, mass matrices very similar to the ones 



already presented in the first model of section 2.1.3J (and proposed in ref 



are obtained. The large atmospheric neutrino mixing is due to a large left handed 
mixing in the lepton sector that corresponds to a large right handed mixing in 
the down quark sector. In the present particular version maximal mixing also for 
solar neutrinos is preferred. A possibly problematic feature of the model is that, 
beyond the unification point, when all the states participate in the running, the 
asymptotic freedom of SU(5i) is destroyed by the large number of matter fields. 
As a consequence, the coupling increases very fast and the theory becomes non 
perturbative below M Pi . In the past models similar to ours have been considered, but 
were discarded just because they contain many additional states and tend to become 
non perturbative between M GUT and M Pi . We instead argue that these features 
are not necessarily bad. While the predictivity of the theory is reduced because 
of non renormalizable operators that are only suppressed by powers of Mqut/ 'A 
with A < Mpi, still these corrections could explain the distortions of the mass 
spectrum with respect to the minimal model, the suppression of proton decay and 
so on. However, it is certainly true that also in this case, as for any other known 
realistic model, the resulting construction is considerably more complicated than in 
the corresponding minimal model. 



4.1 The Model 

The symmetry of the model is SUSY SU(5) (g> U(l). We define Q the charge as- 
sociated with U(l) and assume that the global U(l) flavour symmetry is a good 
symmetry of the superpotential, that is it is not violated by non-perturbative effects 
through operators appearing in the superpotential w. The superpotential of the 
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model can be thought as composed of three parts: 

w = W1+W2 + W3 , (4.1) 

where W\ accounts for the breaking of 577(5), w 2 realizes the doublet-triplet splitting 
by means of the MDM and W3 contains the Yukawa couplings. 

4.1.1 577(5) Breaking 

The W\ term contains only the field Y of the representation 75 of £77(5) with Q = 0: 

Wl = Cl Y £E yEF YfE + M Y y£E y c a E , (4-2) 

where A,B,... = 1,2,..., 5. W\ has the effect of providing Y with a vev (Y) = 
My/(2ci) ~ Mqut and of giving a mass to all physical components of Y, i.e. those 
that are not absorbed by the Higgs mechanism. The 75 uniquely breaks S77(5) down 
to 377(3) <g> S77(2) ® 17(1) through the following vev: 



{YcS) = ~{Y) 



S^ + 28 a c S b d -ls^SE) ~{A+>B) 



2 
where Greek indices assume values from 1 to 3 and Latin indices from 4 to 5. 

4.2 Doublet-Triplet Splitting 

The W2 term induces the doublet-triplet splitting: 

r A v DETj BCFG r 1 „ T7 V BC tj ADEFG 



(4.3) 



W 2 - C 2 n Y BC -H50 tADEFG -- C 3 ll A Y DE ft '50 BCFG € 

150ABCD 



c 4 H 5Q ABCD H- 0ABCD X . (4.4) 



Here H and H are the usual pentaplets of Higgs fields, except that now they carry 
non opposite Q charges. It is not restrictive to take c 2 , C3 and C4 real and positive. 
The renormalizable couplings that appear in w 2 are the most general allowed by the 
S77(5) and Q assignments of the fields in eq. ( |4.4j ) which are given in Table |4.1| . 
The value of q will be specified later. 

4.2.1 Limit of Unbroken SUSY 

At the minimum of the potential, in the limit of unbroken SUSY, the vevs of the 
fields 77, 77, TTso and TTgQ all vanish, while the X vev remains undetermined. These 
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field 


Y 


H 


H 


Hso 


H m 


X 


SU(5) 


75 


5 


5 


50 


50 


1 


Q 





-Q 


q-l 


Q 


l-q 


-1 



Table 4.1: Q charge assignments. 



can be easily understood by writing symbolically the equations for the minimum of 
the superpotential Wi = 0, the subscript i standing for derivation: 



Wy 


= 2M Y Y + 3 Cl Y 2 + c 2 HH 50 + c 3 HH m 


W 5 


= c 2 HY + c 4 XH m 


™m 


= c 3 HY + Ci XH 50 


W H 


= c 2 YH 50 


W H 


= c,YH m 


W X 


= c^H^H-^ 



(4.5) 



For any (X), a minimum is achieved by taking (H) = (H) = (H50) = (-H50) = 
and (Y) as shown in eq. ( |4.3| ). 

As we shall see, when SUSY is softly broken the light doublets in H and H 
acquire a small vev while the X vev will be fixed near the cut-off A, of the order 
of the scale between Mqvt and Mp\ where the theory becomes strongly interacting 
(we shall see that we estimate this scale at around 10 — 20 Mqut, large enough that 
the approximation of neglecting terms of order Mqvtj 'A is not unreasonable). The 
MDM occurs as discussed in eqs. ( |3.9|) and ( 3.1 Op with the only difference that the 



mass term for H 50 and H-^ originates from the vev of X. The U{1) flavour symmetry 
protects the doublet Higgs to take mass from radiative corrections because no HH 
mass term is allowed. Also no non renormalizable terms of the form HHY m X n 
(m, n > 0) are possible, because X has a negative Q charge. As anticipated in 
section |3.2.1| , this version of the MDM was introduced in ref. |)7]] and overcomes 
the observation in ref. p8[ that, in general, non renormalizable interactions spoil 



the mechanism. The Higgs colour triplets mix with the analogous states in the 50 
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and the resulting mass matrix is of the see-saw form: 



rriT 



-4 v / 3c 2 (y) 

-4V3c 3 (F) c 4 (X) 



(4.6) 



Defining m<p = c±{X) the eigenvalues of the matrix rriTrh T are the squares of: 



m,Ti,2 = ^ 



y/ml + 48 (c 2 + c 3 ) 2 (Y> 2 ± ^/m 2 + 48 (c 2 - c 3 ) 2 (Y) 2 J . (4.7) 



and mT\V(iT2 — 48 c 2 c 3 (V) 2 . The effective mass that enters in the dimension 5 
operators with \AB — 1| is l/^m^ 1 )].!: 



rriT 



m Tl m T 2 48 c 2 c 3 (Y)' 



7770 



c 4 (I) 



(4i 



4.2.2 Broken Supersymmetry 

The breaking of SUSY fixes a large vev for the field X by removing the corresponding 
flat direction, gives masses to s-partners, provides a small mass to the Higgs doublet 
and introduces a \i term. In the (50, 50, X) sector, the terms that break SUSY softly 
are: 



C 



soft 



2 1 |2 

m x \x\ 



+ m 50 \c\ + 



2 |-|2 



+ 777 3 / 2 c 4 Ai (cca; + /i.e.) + 



(4.9) 



where c, c and x denote the scalar components of H^o, H-^ and X respectively, 
mj, 77250, "nigQ their soft breaking masses, m 3 / 2 the gravitino mass and A4 a trilinear 
coupling. Note that C4 has here the same role played in SUSY radiative electroweak 
breaking by the top Yukawa coupling. Dots stand for the remaining soft breaking 
terms, including mass terms for the scalar components of the Higgs doublet fields 1 . 
Consider first the limit m x = m 50 = m-^ = 777,3 11 — m an d C4A4 — 1. In the 
SUSY limit and neglecting the mixing between the (50, 50) and the (5, 5) sectors, 
fermions and scalars in the 50 and 50 have a common squared mass c\x 2 . When 
SUSY is broken by the soft terms for each fermion of mass c^x there are two bosons 
of squared masses c\x 2 + m 2 ± mx. The x terms in the scalar potential at one loop 
accuracy are given by: 



V 



2 2 1 

777 x + 



50 
64tt 2 



2(c 4 x +777 +7773;)' 



log 



C 2 X 2 + 777 2 + 7773; 

K 2 



3/2 



1 To find the minima of the scalar potential it is not restrictive to set to zero the imaginary part 
of x, which will be understood in the remaining part of this section. 
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<->/ 2 2 i 2 \2 

2{c A x + m — mx) 



4afx 4 



log 



c\x 2 
"A 2 " 



log 

3/2 



c 2 aX 2 + m 2 



mx 



A 2 



3/2 



(4.10) 



A numerical study of this potential in the limit m C A and for c 4 of order one, 
shows that the minimum is at x of order A (somewhat smaller than A but close to 
it). The expression for V in eq. ( |4.10| ) provides a good approximation of the scalar 
potential only in a small region around A. Outside this region the perturbative 
approximation will break down. We take our numerical analysis as an indication 
that the minimum occurs near A. 

To analyze in more detail the situation we consider the differential equation 
describing the evolution of the x mass and we identify (X) with the scale at which 
m 2 vanishes. With the soft terms of (|4.9|) , we find that the following set of differential 
equations is appropriate to describe the running in range between M GUT and the 
cut-off: 



dm 2 
dlnQ 

dlnQ 

dc4 
dlnQ 
rf(m 3/2 A 4 ) 
dlnQ 
dm 



c\(m 2 x + m 



25 
i^ 2 
dm%r 

dlnQ " 

42 
— — ac 4 + 

07T 



-)0 



+ m m + m m A A ) 



84 

5tt 

13 
4^° 4 

13 



am 2 + 



57T 



2/ 2 1 2 

cAm x + m 



2^4 



50 



+ m m + m s/2 



Al) 



(4.11) 



84 _ 
■— am+ ., 

07T 2TT Z 



m 3/ 2C 4 A 4 



dlnQ 

da 
dlnQ 



52 
— i 
2tt 

52 

2T 



-am 



-a 



where a is the running gauge coupling and m are gaugino masses, which are supposed 
to be degenerate. The value of (X) depends on the boundary conditions assigned 
to the relevant parameters at the cut-off scale, in particular c 4 and the trilinear 
ccx coupling in C so f t . To give a numerical example, if A = 1.1 • 10 17 GeV and 
a(A) = 0.88, then (X) ~ A/4 is obtained by choosing the following reasonable 
values of the parameters at the cut-off scale: m x = m 50 = m-^ = m 3 / 2 = 800 GeV, 

0.2 and A4 = 0.7. This particular situation is shown in fig. |0 

2 



m = 2 TeV, c 4 



where the running of m 



2 ml 



mz 



, ,.. ,,(, 50 , iii^q and m 2 ^ = (m 3 /2C 4 y4 4 ) 2 is displayed. In the next 
section we assume that the true minimum occurs at x — (X) = 0.25A, that is, by 
defining A = (X)/A, at A = 0.25. 
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Figure 4.1: Running of ml,ml ,m^ and rr? m = (7773/2C4A4 



lines correspond to: A 
GeV, m = 2 TeV, c 4 
A = (X)/A ~ 0.25. 



1. MO 17 GeV, a(A) = 0.88, m x 
0.2 (c 4 = 0.4) and A 4 = 0.7. 



The solid (dashed) 

= m 50 = ^50 = m 3/2 = 800 

The solid red line gives a 
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We have studied the dependence of (X) upon both the boundary values of c 4 
and A4. We find that the dependence on 04(A) is logarithmic and that there is 
a mild linear dependence on the trilinear ccx term evaluated at A. The above 
dependence results mainly from the interplay of the differential equations describing 
the evolution of the x mass and the C4 coupling, the quantities with the faster 
running in the region close to the cut-off A. We obtain: 5(X)/A « 0.4 Sc^/c^. This 
means that, with the same boundaries of the previous numerical example but with 
c 4 doubled, the value of (X) becomes greater by a factor 2.2, as can be seen from 



fig. fO|. The vev of X seems thus to be not dangerously unstable with respect to 
the boundaries. This situation differs from the analogous mechanism of radiative 
electroweak symmetry breaking, where a doubling of y t can cause a variation of 
several orders of magnitude in the vev of the Higgs doublet. We have checked that, 
in our case, the mild dependence on C4 is due to the running of c 4 itself which, in 
the appropriate energy window, is not negligible. 

The complex field x describes two physical scalar particles. The one associated 
to the real part of x has a mass of order m.3/2 and couplings to the ordinary fermions 
suppressed by 1/A. The particle associated to the imaginary part of x is massless 
only in the tree approximation. Due to the anomaly of the related U(1)q current 
the particle acquires a mass of order f 7T m n /(X). Cosmological bounds on (X) have 



been recently reconsidered in ref. |)5| where it has been observed that (X) of the 



order of the grand unification scale is not in conflict with observational data. 

An alternative possibility is to assume that the U(l) symmetry is local |67 
In this case the supersymmetric action contains a Fayet-Iliopoulos term and the 
associated D-term in the scalar potential provides a large vev for x, of the order of 
the cut-off scale A. 

A jjl term for the fields H and H of the appropriate order of magnitude can 
be generated according to the Giudice-Masiero mechanism |B9J. Assume that the 
breaking of SUSY is induced by the 9 2 component of a chiral (effective) superfield 
S, singlet under SU(5) with Q = 0. In the Kahler potential a term of the form 

S^X^HH , 
K = A2 + h.c. (4.12) 

is allowed. The vevs of S and X, (S) ~ 9 2 m 3 /2Mpi, (X) = AA lead to an equivalent 
term in the superpotential of the desired form [iHH with \x ~ Ara 3 / 2 Mpj/A which 
can be considered of the right order. 
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4.3 Coupling Unification 

As was pointed out previously, it is well known that in the minimal version of SUSY 
SU(5) the central value of a 3 (m^) required by the constraint of coupling unification 
is somewhat large: a 3 (mz) ~ 0.13 |64], |96], ^, |100| , |101|| . In the model discussed 



here, where the doublet-triplet splitting problem is solved by introducing the SU(5) 
representations 50, 50 and 75, the central value of a^{mz) is modified by threshold 
corrections near Mqut which bring the central value down by a substantial amount 
so that the final result can be in much better agreement with the observed value. As 
discussed in refs. |64|, |96| , this remarkable result arises because the 24 of the minimal 
model is replaced by the 75. The mass splittings inside these representations are 
dictated by the group embedding of the 75 singlet under SU{2>) ® SU{2) Cg> U(l) that 
breaks SU(5). The difference in the threshold contributions from the 24 and the 75 
has the right sign and amount to bring a^imz) down even below the observed value. 
The right value can then be obtained by moving msusY and m? m a reasonable 
range. The difference in the favored value of my with respect to the minimal model 
in order to reproduce the observed value of a^{mz) goes in the right direction to also 
considerably alleviate the potential problems from the bounds on proton decay. Note 
that there are no additional threshold corrections from the 50 and 50 representations 
because the mass of these states arise from the coupling to the field X which is an 
SU(5) singlet. Thus there are no mass splittings inside the 50 and the 50 and no 
threshold contributions. While there is no effect on the value of a^mz) the presence 
of the states in the 50 and 50 representations affects the value of the unified coupling 
at Mqut and also spoils the asymptotic freedom of S77(5) beyond Mqijt- We find 
it suggestive that the solution of the doublet-triplet splitting problem in terms of 
the 50, 50 and 75 representations automatically leads to improve the prediction of 
otz(mz) and at the same time relaxes the constraints from proton decay. We now 
discuss this issue in more detail. 



4.3.1 Estimate of ^(^zOIms 

The calculation proceeds exactly as in the previous chapter. With the same SUSY 
spectrum, all the differences on the values of 05, Mqut, 03, with respect to minimal 
SUSY 517(5), are caused by the term 6i(h) of ( ggg) . 

Following the steps already encountered, £3(2), 63(1) are those of eqs. ( |3.53| ) and 
( |3.55| ), but 83(h) is now obtained substituting in eq. ( |3.56| ) the b\ reported in Table 
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state (SU(S),SU(2),U(1)) 


hf\j) 


b?(j) 


hf\j) 


mass 


(8,3,0) 





8 


9/2 


m Y 


(3, l,5/3)(3, 1,-5/3) 


5 





1/2 


0.8m Y 


(6,2,5/6)(6,2,-5/6) 


5 


3 


5 


0Am Y 


(1,1,0) 











0.2m Y 


(8,1,0) 








3/2 


0.2m Y 


(3, 1,-1/3), (3, 1,1/3) 


1/5 





1/2 


m T1 


(3, 1,-1/3), (3, 1,1/3) 


1/5 





1/2 


m T2 


50,50 


173/10 


35/2 


17 


m^ 



Table 4.2: GUT thresholds for the MDM. m Y = 5M Y . 



For this case one obtains 



S 3 (h) = 2.047 



14vr 



In 



rriT 



M, 



GUT 



+ 



14vr 



In 



m Y 



M, 



GUT 



(4.13) 



where mx is the effective triplet mass. Note that this equation is very similar to eq. 
( 3.57 ), apart from the non-logarithmic term. We can still use the first of eqs. ( p. 49 ), 
with the difference that thresholds from GUT particles within the MDM now give: 



«*<*> = - a993 -n k te 



| 1 ln ( m T 

14 n \M GUT 



(4.14) 



Thus, keeping fixed msusY, the first of eqs. ( |3.49D can be considered as a relation 
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between the GUT masses m v = M GUT ,m Y ,m T . 
M GUT 



In 



: 



m z J V m z 

+ [0.469] (2) + 



25 , f rriQTjQY\ 
-0.036 -In' ] ] 



+ 



-0.993 ln , 

14 V Mgut 



4 \ mz J 

my \ 1 



14 V M G(/T 



(4.15) 



(h) 



Substituting it in 63(h) for Mgut, as pointed out in [3.6. 1| , we can write the 
prediction for £3 in the following form, 



2% \ m z J 



07T 



ttit 



M, 



(0) 



(4.16) 



GUT, 



where again k = k{2) + k(SU SY) + k(M GUT ) and it contains the contribution of two- 
loop diagrams to the running couplings, fc(2), the threshold contribution of states 
near msusY, k(SUSY), and the threshold contribution from states near M G ut, 
k(M G jjT). Remember that the threshold contributions would vanish if all states 
had the mass m S usY or M GUT so that only mass splittings contribute to k(SUSY) 
and k(M G uT). The values of k(2) and of k(SUSY) are essentially the same in the 
minimal and the realistic model. As in the previous chapter we take: k(2) = —0.733 
and k(SUSY) = —0.510, corresponding to the representative spectrum displayed in 
Table [B.l| . The value of k(MauT) was zero for the 24 of the minimal model, but 
now we have k(M G uT) — 1-857 for the 75 of the realistic model. The 75 contains in 
fact non degenerate coloured multiplets. Thus we obtain: 



k 



k = -0.733 - 0.510 = -1.243 
-0.733-0.510 + 1.857 = 0.614 



minimal model 

realistic model (4.17) 



This difference is very important and makes the comparison with experiment 
of the predicted value of at^mz) much more favorable in the case of the realistic 
model. In fact for k large and negative as in the minimal model we need to take 
itT'Susy as large as possible and mr as small as possible. But the smaller is mx, the 
faster is proton decay. The best compromise is something like msusY ~ 1 TeV and 
rriT ~ M GUT , which leads to a s (m z ) ~ 0.13 which is still rather large and proton 
decay is too fast. This is to be confronted with the case of the realistic model where 
k is instead positive and large enough to drag a s (mz) below the observed value. 
We now prefer my to be larger than M GUT by typically a factor of 20-30, which 
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Figure 4.2: Contours of a 3 = a s {mz) in the plane (msusYyfnT), for the missing 
doublet model. The SUSY spectrum is parametrized as in Table |B.l . 
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means a factor of 400-1000 of suppression for the proton decay rate with respect 
to the minimal model. For example, for tususy ~ 0.25 TeV and mj> rj 6 • 10 17 
GeV we obtain a 3 {m z ) ~ 0.116 which is acceptable. The predictions for a 3 versus 



{msusY, %) in the minimal and the missing doublet model are shown in figs. 3.3 



and 4.2 



Due to the large matter content the model is not asymptotically free and the 
gauge coupling constant a 5 blows up at the scale 

M GUT esxp(^\ . (4.18) 

Taking into account that threshold effects modify a 5 with respect the LO value 
a 5 ~ 1/24, we find that the pole occurs near 10 17 GeV, the precise value depending 
on the details of the heavy spectrum. 

4.3.2 Theoretical Uncertainties 

Clearly, there is an uncertainty on k(SUSY), related to the possibility of varying 

both the parameterization used and the relations assumed between the parameters 

mo, wii/2-i A* , iwjj. An overall uncertainty of ±0.005 on a 3 has been estimated in ref. 

100|| by scanning many models consistent with electroweak symmetry breaking, a 



neutral lightest supersymmetric particle and sparticle masses above the experimental 
bound and below ~ 2 TeV. 

We stress that in the present analysis we are assuming universal soft breaking 
parameters at the cut-off scale. If we relax this assumption, then a larger range for 
a 3 could be obtained. Indeed a 3 is particularly sensitive to wino and gluino masses, 
that, at the electroweak scale, are approximately in the ratio 1:3 when a universal 
boundary condition on gaugino masses is imposed. It has been observed that, by 
inverting this ratio, a negative contribution of about —0.01 to a 3 is obtained ||102| . 



Another source of uncertainty on a 3 is related to the unknown physics above 
the cut-off scale of the grand unified theory. There can be threshold effects due to 
new heavy particles or even non perturbative effects that arise in the underlying 
fundamental theory. These effects can be estimated from the non-renormalizable 
operators, suppressed by (V)/A, that split the gauge couplings at the scale Mq(j T 



103|| . For generic, order one coefficients and barring cancellations among different 
terms, the presence of these operators may affect a 3 by additional contributions 
of about ±0.005Mpz/A [|100]| . As we will see the model under discussion requires 



A rj O.lMpi. Therefore, to maintain the good agreement between the experimental 
and predicted values of a 3 , we need a suppression of this contribution by about a 
factor of 10, which we do not consider too unnatural. 
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In view of the large theoretical uncertainties on 03 we cannot firmly conclude 
that the gauge coupling unification fails in the minimal model while it is completely 
successful in the realistic one. However we find very encouraging that, by solving 
the doublet-triplet splitting problem within the missing partner model, acceptable 
values of 03 can be easily obtained and that they are directly related to a proton 
lifetime potentially larger than in the minimal model. 

4.4 Yukawa Couplings 

The W3 term contains the Yukawa interactions of the quark and lepton fields \&io, 
\&5 and \l/i, transforming as the representations 10, 5 and 1 of SU(5) respectively. 
We assume an exact .R-parity discrete symmetry under which \&io, ^5 and \l/i are 
odd whereas H, H, H 50 and H-^ are even. The W3 term is symbolically given by 

^3 = ^i G u (X,Y)^ 10 H + * 10 G d (X,Y)V z H + ^G^X^^H 

+ Mf 1 G M (X,Y)V 1 + y w G m (X,Y)V 10 H m . (4.19) 

The Yukawa matrices G u , Gd, G u , Gm and G-^ depend on X and Y and the asso- 
ciated mass matrices on their vevs. The last term does not contribute to the mass 
matrices because of the vanishing vev of H-^q, but is important for proton decay. 
The pattern of fermion masses is determined by the U(l) flavour symmetry that 
fixes the powers of A = (X)/A for each entry of the mass matrices. In fact X is 
the only field with non vanishing Q that takes a vev. The powers of A in the mass 
terms are fixed by the Q charges of the matter field \l/ and of the Higgs fields H and 
H. We can then specify the charge q that appears in Table JO] and the Q charges 
of the matter fields \1/ in order to obtain realistic textures for the fermion masses. 
We choose q = 2, so that we have: 

Q(H) = -2 and Q(H) = 1 , (4.20) 

and, for matter fields 

Q(* 10 ) = (4,3,l) , Q(* s ) = (4,2,2) , Q(^) = (1, -1,0) . (4.21) 

The Yukawa mass matrices are of the form: 

G r ((X),(Y)) ij = A n «G r (00) o -, r = u,d,u,M . (4.22) 

We expand G r ((Y))ij in powers of (Y) and consider the lowest order term at first. 
Taking G r (0)ij of order 1 and n^ as dictated by the above charge assignments we 
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vA 4 
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maj 



A 2 1 A" 




1 


M . 


A 1 




(4.23) 


id A ~ Ac 


w 0.22, 



For a correct first approximation of the observed spectrum we need A 

Ac being the Cabibbo angle. These mass matrices are very similar to those already 



introduced in section |2.1.3| , with two important differences. 

We have here that tan/3 = v u /vd ~ mt/rribX 4 , which is small. The factor A 4 is 
obtained as a consequence of the Higgs and matter fields charges Q, while in chapter 
|2] the H and H charges were taken as zero. We recall that a value of tan/3 near 1 
is an advantage for suppressing proton decay. A small range of tan/? around one is 
currently disfavored by the negative results of the SUSY Higgs search at LEP [98|. 
Of course we could easily avoid this range, if necessary. 

Remember that the zero entries in the mass matrices of the neutrino sector occur 
because the negatively Q-charged X field has no counterpart with positive Q-charge. 
Neglected small effects could partially fill up the zeroes. As explained in ref. | 32[ 
these zeroes lead to near maximal mixing also for solar neutrinos. At variance of 
the example of section |2.1.3| , due to the different charge assignments for ^^ 1 , the 
LOW solution to the solar neutrino problem emerges. 

As pointed out, a problematic aspect of this zeroth order approximation to the 



mass matrices is the relation rrid 



mt 



which is good as an order of magnitude 



relation because relates large left handed mixings for leptons to large right handed 
mixings for down quarks |32], [0],|55|, ^(| [I7|] . However the implied equalities mb/m T = 
m s/""V — md/me = 1 are good only for the third generation while need to be 
corrected by factors of 3 for the first two generations. As already mentioned in the 
discussion about the mass matrices presented in section |3.5| , the necessary corrective 
terms can arise from the neglected terms in the expansion in (Y) of G r ((Y))ij 
57| . The higher order terms correspond to non renormalizable operators with the 
insertion of n factors of the 75, which break the transposition relation between m^ 
and m e . For this purpose we would like the expansion parameter (Y)/A to be not 
too small in order to naturally provide the required factors of 3. We will present in 



2 In our convention Dirac mass terms are given by LmR and the light neutrinos effective mass 



matrix is m u m ma jml 
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the following explicit examples of parameter choices that lead to a realistic spectrum 
without unacceptable fine tuning. 

4.4.1 A Numerical Example 

To reproduce the fermion mass spectrum we must further specialize the Yukawa 
couplings introduced in eq. ( |4.19| ): 



W 3 = |*10 G u *io H e ABODE + ""^"^lO G 50 *10 H 50 ABC D 

+ V2 K B G d ^- 5A Hb + \V2 *#* F d m- bc H D Y2E 

+ *5 A G t ,* 1 H A -±MV 1 G M *i + ... (4.24) 



where G r (r = u, d, v, M, 50) is proportional to G r ((X},0) of eq. ( |4.19| ). We have 



explicitly introduced a term linear in Y, whose couplings are described by the (X)- 
dependent matrix Fd. There are other terms linear in Y/A, not explicitly given 
above. In particular we may insert Y/A in the renormalizable term providing masses 
to the up type quarks. We neglect such a term since, on the one hand, the matrix 
G u is already sufficient to correctly reproduce the up quark masses and, on the other 
hand, this operator would not significantly modify the results for proton decay. As 
we will show, Y/A is close to 0.1 in our model and higher order terms in the Y/A 
expansion can be safely neglected. The interaction term involving H-^ ABCD does 
not contribute to the fermion spectrum, but it will be relevant for the proton decay 
amplitudes. 

The term linear in Y in the previous equation is sufficient to differentiate the 
spectra in the charged lepton and down quark sectors. We get the following Dirac 
mass matrices: 

Vu Vd (A OK\ 

m U)U = y U)V — = , m d)6 = y dfi — , (4.25) 

Vu = G U y v = G u , 

(Y) (Y) (426) 

y d = G d + ^F d y T = Gd -3^F d , 

where v u , v d and (Y) parameterize the vevs of H, H and Y respectively. The fermion 
spectrum can be easily fitted by appropriately choosing the numerical values of the 
matrices G u , G d , F d , G v and Gm- The most general fitting procedure would leave 
a large number of free parameters. Here we limit ourselves to the discussion of one 
particular example. In agreement with eq. (|4.23| ), we take: 
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G, 



;-0.51 + 0.61i)A 6 
(0.42-0.70i)A 5 
(0.27 + 0.86i)A 3 



(0.42 - 
;-0.39 
'-0.30 



0.70i)A 5 
F0.52i)A 4 
- 1.14i)A 2 



(0.27 + 0.86i)A 3 

;-0.30-1.14i)A 2 

1.39 



(4.27) 



Ga — A 



(2.39 — l.lli)A 5 

(0.87 + 0.55i)A 4 
;-1.50 + 0.94i)A 2 



(0.33 - 0.59i)A 3 

(2.76 + 0.89f)A 2 

(0.45 + 1.78a) 



(0.13 + 0.45i)A 3 

(0.69 - 0.5l/)A 2 

1.94 



(4.28) 



00 

A 



A J 



(0.38-0.18i)A 5 
(-0.08-0.05i)A 4 
(-0.09 + 0.12i)A 2 



(-0.16-0.06i)A 3 (0.07 + 0.04i)A 3 
(-0.20-0.15i)A 2 (0.15 - 0.11i)A 2 
(0.07 + 0.14i) -0.19 



(4.29) 



G h 



-0.78-0.19i)A 3 
-1.23 - 0.34z)A 

(0.45 + 1.18i)A 



(0.52 - 0.34i)A 






;i.38 + 0.39i)A 2 
\l.04+1.31z) 
0.8 + 1.2i 



(4.30) 



G 



M 



(1.50 + 0.55z)A 2 

(1.41 + 1.19z) 

(0.35 - 1.53i)A 



;i.4i + i.i9i) 







(0.35 


-1.53i)A" 







1.26 + 1.48z 



(4.31) 



While the generic pattern of G u , G d , F d , G v and Gm is dictated by the U(l) 
flavour symmetry, the precise values of the coefficients multiplying the powers of A 
are chosen to reproduce the data. We take A = (X)/A = 0.25, tan/3 = 1.5 and 
M = 0.9 • 10 15 GeV. In SU(5) the matrix G u contains two additional phases [5B|, 
0i and 02 that have been set to zero in eq. G4.27T ). These phases do not affect 
the fermion spectrum but enter the proton decay amplitude. When discussing the 
proton decay we will analyze also the dependence on <\>\ and 02- 

From the above matrices we obtain, at the unification scale: 



(4.32) 



t = 200 GeV m c 


= 0.27 GeV m u -- 


= 0.9 MeV 


b = 1.0 GeV m s 


= 26 MeV m d = 


= 1.1 MeV 


T — 1.1 GeV m„ 


= 71 MeV m e = 


= 0.34 MeV 



|K S | = 0.22 \V ub \= 0.0022 |Vrf,l = 0.052 J = 1.9 • 10~ 5 , (4.33) 

where J is the CP- violating Jarlskog invariant. In the neutrino sector, we find: 



m 1 = 0.81 • 10~ 3 eV , m 2 = 0.88 • 10" 3 eV , m 3 = 0.061 eV 



(4.34) 
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More precisely: 

Am 2 sol = m\-m\ = 1.1- 1(T 7 eV 2 , Am 2 tm = m\-m\ = 3.7- KT 3 eV 2 .(4.35) 
The neutrino mixing angles are 

#i2~^ , #23~^ , #13 = 0.06 . (4.36) 

Until now we have not specified the value of (Y)/A. We know that (Y) should be 
around Mq{j T . The cut-off A cannot be too close to (Y), otherwise most of the 
spectrum of the model would lie beyond the cut-off. At the same time A cannot 
be too large: it is bounded from above by the scale at which the SU(5) gauge 
coupling blows up, which, as we see from eq. ( |4.18|) , occurs more or less one order 



of magnitude below the Planck mass. This is welcome. If, for instance, we take 
(Y)/A = 0.05 — 0.1, it is reasonable to neglect multiple insertions of Y in the Yukawa 
operators. At the same time, in the example given in eq. ( |4.29| ), the coefficients of 
the powers of A in Fd remain of order one, even if (Y) / A is as small as 0.05. A too 
large cut-off would have been ineffective in separating down quarks from charged 
leptons unless we had chosen unnaturally large coefficients in the allowed operators. 
What is usually considered a bad feature of the missing partner mechanism - the 
lack of perturbativity before the Planck mass - turns out here to be an advantage 
to provide a correct description of the fermion spectrum. 

In the neutrino sector, we obtain a bimaximal neutrino mixing with the so- 
called LOW solution to the solar neutrino problem Q. Within the same U(l) flavour 
symmetry considered here we could as well reproduce the quasi-vacuum oscillation 
solution, by appropriately tuning the order one coefficients in G u and Gm- We recall 
that the value of M required to fit the observed atmospheric oscillations is probably 
somewhat small in the context of SU(5), where a larger scale, closer to the cut-off 
A, is expected. This feature might be improved by embedding the model in 50(10) 
where M is directly related to the B — L breaking scale. 

In conclusion, the known fermion spectrum can be reproduced starting from a 
superpotential with order one dimensionless coefficients. Mass matrix elements for 
charged leptons and down quarks match only within 10 — 20%, due to (Y)/A m 0.1, 
and this produces the required difference between the two sectors. The neutrino 
mixing is necessarily bimaximal in our model, with either the LOW or the vacuum 
oscillation solution to the solar neutrino problem. 



3 As pointed out, in chapter 1, latest preliminary results from Super-Kamiokande (77J], including 
constraints from day-night spectra, seem in fact to prefer bimaximal neutrino mixing and to revamp 
interest in the LOW solution. 
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4.5 Proton Decay 

Similarly to the case of minimal SU(5), we expect that the main contribution to 



proton decay comes from the dimension five operators []79| originating at the grand 
unified scale when the colour triplet superfields are integrated out |78], |80], |5I|. We 
denote the colour triplets contained in H, H, H 5Q 



#50 by H 3u , H 3d , 



H' 3d and H' 3u , 



respectively. The part of the superpotential depending on these superfields reads: 



</.' 



-QG U Q + U C G U E C 



= H$u 

+ H 3d 

+ H! iu 

- AV3c 2 (Y) H' 3d H ?ya - 4v / 3c 3 (F) H 3d H' 3u + c 4 (X) H' 3d H' 3u 



QCL + U C DD C 



\QG m Q + U c G m E< 



where 



C 



~G n 



(Y) 



D = G n 



00 



... (4.37) 



(4.38) 



A _u ' " -" A " 

and Q, L, U c , D c and E c denote as usual the chiral multiplets associated to the 
three fermion generations. Notice that, at variance with the minimal SU(5) model, 
an additional interaction term depending on H' 3u is present, which originates from 
the interaction of 50 with the 10 matter fields. In Appendix |C| we give a short 
discussion on its properties and justify the choice of the numerical coefficient v3/2 
in eq ( |4.24j ). By integrating out the colour triplets we obtain the following effective 
superpotential: 



w 



eff 



1 

ttit 



QAQ QCL + U C BE C U C DD C 



+ 



(4.39) 



where m^ has been defined in eq. 



B 



4y/3c 2 (l% 
&u+ c 4 (X) { ' 



50 



-2A 



G u 



2^02 (y) 

ca(X) 



G 



50 



(4.40) 



and dots stand for terms that do not violate baryon or lepton number. Minimal 

? d = 0. As discussed, in that case the matrices 



SU(5) is recovered by setting G-^j = F 

A, B, C and D are determined by G u (in which now also <pi and 

G d and therefore strictly related to the fermionic spectrum |78 

a distortion due to the terms proportional to G-^q and Ed. These distortions have 



2 play a role) and 
In our case we have 
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different physical origins. On the one hand the terms containing F^ are required to 
avoid the rigid mass relation of minimal SU(5). They are suppressed by (Y)/A and 



we expect a mild effect from them |^| . On the other hand the terms containing G-^ 
are a consequence of the missing doublet mechanism and they might be as important 
as those of the minimal model. 

To estimate the proton decay rates we should specify some important parameters. 
First of all the mass tut- From the discussion of the threshold corrections we know 
that a large value for m-r is preferred in our model. We should however check 
that this value can be obtained with reasonable choices of the parameters at our 
disposal. The unification conditions constrain (Y) in a small range around 10 16 GeV. 
More precisely, after the inclusion of threshold corrections and two loop effects, eq. 
( [4.1 5|) becomes a relation among the masses of the super-heavy particles and the 
leading order quantity Mq(j T . This relation limits the allowed range for the heavy 



gauge vector bosons and indirectly, through eq. flCJ.41 ), pushes the SU(5i) breaking 
vev (Y) close to 10 16 GeV. Then the cut-off scale and the X vev are fixed by the 
phenomenological requirements (Y) ~ 0.05 A and A = (X)/A = 0.25. Large values 
of rriT could be obtained either by taking m^m^ large or by choosing a small m^. 
This last possibility is however not practicable, since the gauge coupling a 5 blows 
up at approximately 20 mf it is not reasonable to push m^ below 5 • 10 15 GeV. The 
only remaining freedom to obtain the desired large value for rrir is represented by the 
coefficients c 2 and C3 that, however, cannot be taken arbitrarily large. We should also 
check that all the heavy spectrum remains below A and this requirement imposes a 
further constraint on our parameters. A choice that respects all these requirements 
is provided by []: 

Cl = 0.035 , c 2 = c 3 = 2.8, c 4 = 0.7 , (4.41) 

(X) = 3.0 • 10 16 GeV, (V) = 5.7 • 10 15 GeV (A = 1.2 ■ 10 17 GeV) , (4.42) 
which leads to: 

M GUT = 2.9 • 10 16 GeV , m = 2.0 • 10 16 GeV , (4.43) 

m Tl = 1.2 • 10 17 GeV , m Ta = 1.0 • 10 17 GeV , m T = 6 • 10 17 GeV . (4.44) 
The heavy sector of the particle spectrum is displayed in fig. (O. 



To evaluate the loop function / we also need the spectrum of the supersymmetric 
particles. As an example we take here the same spectrum considered in Table p.l 



4 The values of the coefficients C2 and C3 here adopted raise doubts on the validity of the per- 
turbative approach that has been exploited in several aspects of the present analysis. We adhere 
to this choice also to show the difficulty met to obtain acceptable phenomenological results within 
a not too complicated scheme. 
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Figure 4.3: Heavy sector of the spectrum. 
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with msusY = 250 GeV. This leads to a squark mass of about 800 GeV, a slepton 
mass of approximately 350 GeV a wino mass of 250 GeV and a charged higgsino 
mass of 125 GeV. 

The matrix G-^ is not directly related to any accessible observable quantity and 
the only constraint we have on it comes from the U(l) flavour symmetry that requires 
the following general pattern: 



G 



.-)0 



A 7 


A 6 


A 4 ] 


A 6 


A 5 


A 3 


A 1 


A' 3 


A 



(4.45) 



With the above numerical values we also obtain: 

2 v^ c 2 (Y) 



ca(X) 



2.7 



(4.46) 



The texture for Ggjj has an overall suppression factor A compared to G u . Therefore 
the contribution of G-^ to the matrices A and B is comparable or even slightly larger 
than the minimal contribution provided by G u . The dominant operators are then 
the ones already described and the calculation can be carried out following Table 

o. 

At variance of the minimal case, here the amplitude with G u and the one with 
G50 interfere. This interference can be either constructive or destructive, depending 
on the relative phases between the two terms. We have scanned several examples 
for GgQ, obtained by generating random coefficients for the order one variables in 



eq. ( 1.45J ). By keeping fixed all the remaining parameters we obtain a proton decay 
rate in the range 2 • 10 32 -r- 4 ■ 10 34 yr for the channel p —>■ K + V and a rate between 
3 • 10 32 yr and 5 • 10 34 yr for the channel p — ► ti + V. For comparison, considering 
the same choice of parameters but setting G§q = Fa = 0, we obtain 9 • 10 32 yr 
and 2 • 10 33 yr respectively, for the above channels. The present 90% CL bound on 
r/BR(p — > K + v) is 1.9 10 33 yr [fnj. These estimates have been obtained by setting 
to zero the two physical phases <fri, 02 contained in the matrix G u . These additional 
parameters may increase the uncertainty on the proton lifetime. For instance, in 
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minimal 577(5), the proton decay rates for the channels considered above, change 
by about one order of magnitude when (pi and 02 axe freely varied between and 
2-7T. Even when the inverse decay rates for the channels K + v and 7r + z/ are as large 
as 10 34 yr, they remain the dominant contribution to the proton lifetime. Indeed, 
since the heavy vector boson mass, Mqut, is equal to 2.9 10 16 GeV in our model, 
the dimension 6 operators provide an inverse decay rate for the channel e + 7r° larger 
than 10 36 yr. 

The effective theory considered here breaks down at the cut-off scale A. We 
expect additional non-renormalizable operators contributing to proton decay am- 
plitudes from the physics above the cut-off. We have already stressed that, by 
assuming dimensionless coupling constants of order one, in unified models without 
flavour symmetries the proton lifetime induced by these operators is unacceptably 
short, even when A = M Pi |S3|, [5TJ. In our case these contributions are adequately 
suppressed by the U(l) symmetry. If we compare the amplitude A nr induced by the 
new non-renormalizable operators with the amplitude A% coming from the triplet 
exchange, we obtain, for the generic decay channel, 

A nr ^ A8 c 2 c 3 f (Y)\ \ i (44?) 



A 3 c 4 \ A 

This supports the conclusion that a proton lifetime range considerably larger than 
the one estimated in minimal models is expected in our case. 



Conclusions 



We have considered the problem of obtaining the large atmospheric neutrino mixing 
in the framework of GUT models. 

First, we have analyzed the most interesting mechanisms which provide this large 
mixing. 

Usually the large atmospheric mixing arises because of the presence of large 
mixings in the Dirac and/or Majorana right mass matrices. In GUT models this sit- 
uation is naturally realized by exploiting the possible presence of large right mixings 
for down quarks, which correspond to large left mixing for charged leptons. 

We have pointed out an alternative mechanism, which does not require the pres- 
ence of any large mixings. In fact, even when both the Dirac and Majorana right 
matrices possess very small mixings, the large atmospheric mixing entirely arises 
from the see-saw mechanism if two conditions apply, that is if the spectrum of right 
Majorana masses is sufficiently splitted and if the left handed v^ and v T couple with 
comparable strength to the right handed v c . As opposite to the first mechanism, 
this situation can be realized even if all left and right mixings of quarks and leptons 
are very small. This feature is particularly attractive from the theoretical point of 
view, because it implies that this mechanism is well compatible with those GUT 
models where right and left mixings are similar. 

Despite the fact that GUT has been studied since a long time, still a lot of its 
aspects require a better understanding. It would be really improbable for this idea 
to survive if inside it there where not something true. But it is clear that minimal 
versions of SUSY GUTs are already ruled out because unrealistic. We have addressed 
the question whether it is possible to construct a not too complicated, relatively 
complete and realistic model which can serve as benchmark to be compared with 
experiment. This model must also correctly reproduce neutrino masses and mixings. 

We have then constructed an example of SUSY SU(5) GUT model, with an ad- 
ditional U(l) flavour symmetry, which is not plagued by the need for large amounts 
of fine tunings, like those associated with doublet-triplet splitting in the minimal 
model, and leads to an acceptable phenomenology. This includes coupling unifica- 
tion with a value of 0:3 (m^) in much better agreement with the data than in the 
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minimal version, an acceptable pattern for fermion masses and mixing angles, also 
including neutrino masses and mixings, and the possibility of a slower proton decay 
than in the minimal version, compatible with the present limits (in particular the 
limit from Super-Kamiokande of about 2 • 10 33 yr for the channel p — » K + v). In the 
neutrino sector the preferred solution is one with nearly maximal mixing both for 
atmospheric and solar neutrinos. 

The U(l) flavour symmetry plays a crucial role by protecting the light doublet 
Higgs mass from receiving large mass contributions from higher dimension opera- 
tors, by determining the observed hierarchy of fermion masses and mixings and by 
providing an adequate suppression for non renormalizable operators inducing proton 
decay. Of course, the U(l) symmetry can only reproduce the order of magnitude of 
masses and mixings, while more quantitative relations among masses and mixings 
can only arise from a non abelian flavour symmetry. 

A remarkable feature of the model is that the presence of the representations 
50, 50 and 75, demanded by the Missing Doublet Mechanism for the solution of the 
doublet-triplet splitting problem, directly produces, through threshold corrections at 
Mqut from the 75, a decrease of the value of a 3 (m^) that corresponds to coupling 
unification and an increase of the effective mass that mediates proton decay by 
a factor of typically 20-30. As a consequence the value of the strong coupling 
is in better agreement with the experimental value and the proton decay rate is 
smaller by a factor 400-1000 than in the minimal model. The presence of these large 
representations also has the consequence that the asymptotic freedom of SU(5) is 
spoiled and the associated gauge coupling becomes non perturbative below Mpi. 
We argue that this property far from being unacceptable can actually be useful to 
obtain better results for fermion masses and proton decay. 

Clearly such a model is not unique: our version is the simplest realistic model 
that we could construct. We think it is interesting because it proves that SUSY 
SU(5) GU is not excluded and offers a benchmark for comparison with experiment. 
For example, even including all possible uncertainties, it is difficult in this class of 
models to avoid the conclusion that proton decay must occur with a rate which is 
only a factor 10-50 from the present bounds. Failure to observe such a signal would 
require some additional specific mechanism in order to further suppress the decay 
rate. 

Finally, it is a generic feature of realistic models that the region between Mqjjt 
and Mpi becomes populated by many states with different thresholds and also non 
perturbative phenomena occur. This suggests that the reality can be more compli- 
cated than the neat separation between the Grand Unification and the string regime 
which is postulated in the simplest toy models of GUTs. 

There is then the need for further ideas on the theoretical side. It would be 

100 



of crucial importance to actually understand the origin of the flavour symmetry. 
The simple U(l) used in the present model is satisfactory at the level of orders 
of magnitude, but is not able to really predict the observed values of masses and 
mixing. Maybe the flavour has a deeper origin, in a fundamental theory of particles 
that includes also gravity beyond the strong and electroweak interactions. 
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Appendix A 

Formulae for the Calculation of 
Proton Lifetime 

A.l d=5 B Violating Operators 

The starting point is the most general renormalizable superpotential which could 
come from S77(5). Yukawa couplings of doublets and triplets will be considered un- 
corrected as well as down quarks and charged lepton mass matrices. The discussion 
can thus be adapted for any 577(5) based model. Omitting colour and weak isospin 
indices, in the interaction basis the relevant terms read: 

w 3 Q T y u U c H 2u + Q T y d D c H 2d + L T y e E c H 2d 

+ Q T AQH 3u + U cT BE c H 3u + Q T CLH 3d + U cT DD c H 3d (A.l) 

+ m T H 3u H 3d , 

where Q, L, U c , D c and E c denote as usual the chiral multiplets associated to the 
three fermion generations, y u ,d,e, A, B, ... are matrices in generation space and % is 
the triplets mass. By integrating out the colour triplets one obtains the following 
effective superpotential 

1 \qAQ QCL + U C BE C U C DD C 



W eff = 



+ ... , (A.2) 



dots standing for terms that do not violate baryon or lepton number. 

At lower scales these dimension five operators give rise to the four-fermion op- 
erators relevant to proton decay, via a "dressing" process. When considering the 
operator QQQL, the most important dressing comes from wino exchange Q. Charged 

1 Gluino dressing contributions cancel among each other in case of degeneracy between first two 
generations of squarks [g3|, Q . 
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0^=^-(y l d i )(d j u k )[L^A + A T )L d ] [LlCL e ] kl (U(u,e) + U(u,d)) 
Of = ^-(u i d j )(d k ui)[Ll(A + A T )L d UL^CL e } kl (U(u,e) + U(u,d)) 



o L 



rriT 

-.2 



[u^MiLKA + A T )L d ] ti [LlCL e } kl (U(d, v) + U(u, d)) 



Of = ^-(u 3 d k )(u iei )[Ll(A + A T )L U ] [L T d CL e } kl {U{lu) + U(d, v)) 

Tflj' J 



)R 



^1 ~~ 


m T 


o? = 


1 


rriT 


of = - 


1 


rriT 


of-- 


1 





(ta) (t#$ [L + d ylBytLl} kl [RlDR+^f^, e c ) 

(aw) Ken k^+l [l+^^+l:] . ./^(« c , j* ) 



m T 



^■4)(ex)[^^] fc jK^w^' 



c u c ) 



Table A.l: Four-fermion operators arising from dimension 5 operators. For each 
operator it is understood to also add the h.c. 



higgsino exchange provides instead the most important dressing of U C E C U C D C 

The integration over the loop has the effect of producing the multiplicative factor 
g 2 fw(d,b) for wino dressing (g being the SU(2) gauge coupling) and f^ia, b) for 
higgsino dressing. / is a function coming from the loop integration: 



fc(a, b) 



d 4 k 



(27r) 4 k 2 — m\ k 2 — m 2 ft — m 5 



VTir 



16tt 2 (m 2 - 



m~ 



m- 



2 2 6 2 

mi — mi mi 



log- 



m- 



mi mi 

2 b 2 lQg^ 

mi — mi mi 



. (A.3) 
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The operators are explicitly given in Table [A.l| , where we made the not too bad 
assumption of degeneracy between the generations of each sparticle (the dependence 
on sparticle masses is logarithmic), so that the final expression is more compact. 
We also turned to the mass eigenstate basis by indicating with L u ^ e and R u ,d,e the 
unitary matrices that diagonalize the fermion mass matrices: 

L T u yuRi = yi d) ~L T d y d R\ = yf -L T e y e R\ = y^ , (A.4) 

with (yf)(d) (f — u,d,e) diagonal and positive. The quark mixing matrix is then 
Vckm = L\L d . 

A. 2 The Relevant Operators 

As discussed in chapter ||, the most important operators are Of and Of, which are 
explicitly given by: 

0[ = Cf jkl (uidi) (djU k ) + h.c. 

O? = C* jkl {d k vi) (vFifrj) +h.c. , (A.5) 

where i, j, k, I are generation indices and 

c ijki = —9 2 [U(u, d) + U(u, e)] P i:j Q M 

Cijki = f~h\ uC i eC )SiiTij (A. 6) 



P = 2L T d AL d Q = L T U C L e 

S = L\ y* B yt L* T = R* D R.I 



(A.7) 



Considering Of, for the channel p — > K + u we have to analyze the contributions 
from both [yid) (su), which corresponds to C^u, and (z/js) (du), which corresponds 
to Of li; . Since C^ kl is proportional to PijQki, it is possible to identify immediately 
the important P, Q elements. For instance C± 21l is linked to P\2 and Qu. Note that 
P is symmetric, hence C\ 2 u — Of li; . For minimal 5*?7(5), P, Q greatly simplify, as 
can be seen from Table 13.21. In this case 



P12 = [V T yi d) K + V} 12 ~ y top X 5 , (A.8) 

and Q u = [V*y d %i, so that 

Q12 ~ Qis ~ yboA 3 > Qu . (A.9) 
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Thus, since 

VtopVbot ~ A 2 (tan/3 + l/tan/3) , (A. 10) 

one obtains 

Pi2Qu(3) ~ A 10 (tan/3 + l/tan/3) , (A.ll) 



as shown in Table |3]3|. For the channel p — > -k + u, we have only to consider (yid)(du) 



that is PiiQu- Since P u ~ APi 2 , then 

P u Qi2(3) ~ A n (tan/3 + l/tan/3) . (A.12) 

In the case of Of , the channel with K + is obtained for (dfi) (u c s c ), corresponding 
to S\{Ti2, and (sfj) (u c d c ), corresponding to S^Tn. In the minimal case, from 



>2, 



S=V + y^Vyf , (A.13) 



it turns out that the dominant contribution corresponds to / = 3, that is emission 
of u T , and S 13 ~ yt op y bo A 3 , S 23 ~ y 2 op y bot \ 2 . From 

T = KV*yf ] (A. 14) 

it follows T12 ~ 2/botA 3 , Tn ~ yboA 4 - Thus, 

5i 3 Ti 2 ~ y t 2 op2 /LA 6 ~ A 10 (tan/3 + l/tan/3) 2 ~ 5 23 T n . (A.15) 

For the channel with 7r + , we have to consider {dv{) (u c d c ), that is 

S 13 T n ~ A n (tan/3 + l/tan/3) 2 . (A. 16) 



These results are given in Table |3T3 . 



Note that we have also introduced Kl,r, which are constants accounting for the 
renormalization of the operators from the grand unification scale down to 1 GeV 



78| , |81| , [104J . In our estimates we take Ki = K 2 = 10. 

The next step is to parameterize the partial rates according to the results of a 
chiral Lagrangian computation [ |105|| . The rates for the dominant channels are given 



in Table Kl 



channel 


rate 


K + u, L 


x K p 


jS^W(l + ^P + 3F))c^ 


2 


K + u T 


*k P [f S^ma + (l + SkiP + 3F)) Gf 113 


2 


k + ^ 


X ff |/3(l + Z> + F)2C£ 12 | 2 


tt + v t 


X w |/5 (1 + D + F) 2C& 13 + a (1 + D + F) 2Cf 311 | 2 



Table A. 2: Proton decay rates. We define X K<7r = (mi — rn^ ^) 2 / (327rrripf^)] m p , 
rriK, ttItt are the proton, K + and 7r + masses; rriB is an average baryon mass, f n 
is the pion decay constant; D and F are coupling constants between baryons and 
mesons in the relevant chiral Lagrangian; /3 and a parameterize the hadronic matrix 
element. In our estimates we take m p = 0.938 GeV, itlb = 1.150 GeV, tjik = 0.494 
GeV, m n 
GeV 3 I88H 



0.140 GeV, U = 0-139 GeV, D = 0.8, F = 0.45 P§, (3 = -a = 0.003 



Appendix B 

Representative SUSY Spectrum 



SU{2) ® U(l) breaking effects are neglected. Numerical coefficients appearing in the 
previous Table are due to renormalization. 

The additional freedom related to the parameters m , miw, fi, van is here fixed 
by choosing 0.8m = 0.8rrii/2 = 2/i = van and taking as a definition msusY = ttih, 
so that all particle masses can be expressed in term of msusY- This parameterization 
leads to k(SUSY) = -0.510. 
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sparticle 


mass 2 


gluinos 


(2.7m 1/2 ) 2 


winos 


(0.8m 1/2 ) 2 


higgsinos 


^ 


extra Higgses 


m\ 


squarks 


ml + 6m 2 , 2 


(sleptons)^ 


ml + 0.5m 2 /2 


(sleptons)# 


ml + 0.15m 2 /2 



Table B.l: Representative SUSY spectrum. 



Appendix C 

Group Theory for the MDM 



In this appendix we summarize the method adopted to deal with tensors in SU(5). It 
is convenient to write large representations like 75, 50 and 50 as products of smaller 
representations, like 10 and 10. In this way, it is not difficult to recognize how the 
multiplets of the SM are embedded in those large representations, thus staying in 
touch with the physical role of the various superpotential terms. In particular, when 
considering the interaction terms, it is compelling to accurately take into account 
the numerical coefficients which follow directly from the tensorial structure of the 
representations involved. Since these coefficients can be large, it is not possible to 
omit them when performing a quantitative analysis. 



C.l Basic Conventions for the Tensors in SU(5) 

In SU(5) the matter fields are contained in just the 10 and 5 representations. We 
write them according to 



*m = —^ 



1 

V2 






m 


-u c 2 


u 1 


-W 





U{ 


u 2 


u c 2 


-m 





u 3 


-u 1 


-u 2 


-u 3 





-D 1 


-D 2 


-D 3 


-E c 



D 1 
D 2 
D 3 

E c 




(C.l) 



#« 



E 



—v 



(C.2) 
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The 10 representation tyf^ is thus represented by a tensor whose two indices A,B = 
1, ..., 5 are antisymmetric. Note that the SU(3) indices correspond to A, B = 1, 2, 3 
while the SU(2) ones are selected for A, B — 4, 5. To remember about this in an 
explicit way, we will sometimes write a (a) instead of A if A = 1, 2, 3 (4,5). Due 

it turns out that it is convenient to deal just with 
The 5 representation is instead written as a tensor with one low index 

In our conventions, the Higgs fields H a, H a belonging respectively to the 5 and 
5 representations are written as 



to the antisymmetry of ^iq, 






\A<B 



H 



H?,d2 
Hzd3 
H 2d 
n 2d. 



H 



hi 

HI 



(C.3) 



where (H* d ) = v d /V2, (Hi) = v u /^2. 
Let us underline a few points. 

• The factor l/v2 is introduced in the definition of \l/io in order to have the 
kinetic terms in the canonical form. In fact, the Kahler potential contains 



K 3 ^ w 



AB* 



(C.4) 



where summation over repeated indices has to be understood. Due to the 
antisymmetry of \J/ 10 under A <-» B, one immediately realizes that: 



< B < B * 



9 \Jl AB \V AB * I 

2 W 10 W 10 \A<B 



U c n U c * + U a U a * + D a D a * + E C E C 



(C.5) 



where a is the colour index. Thus, the fields U c , U, D, E c have a canonical 
kinetic term. On the other hand the pentaplet ^/g^ has no problems of multi- 
plicity, so that there is no need of introducing any coefficient in the definition 
of (|C.2|) to get the canonical form for the kinetic terms of D c , E, v. 



At the contrary of what happens for all the other \l/f B a<b-> we have introduced 
a minus sign in \1/Jq. The reason is the following. The mass term for up quarks 
(considering for simplicity one generation only) arises from the superpotential 
term 



w 3 -K B ^w D H E e A BCDE = 

x t,AB x t,CD ttE^ I 

_ ^10 ^10 n e ABCDE\A<B,C<D 



(C.6) 
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with E = 5. When one picks the indices corresponding to U%U 2 , that is 
(A,B,C,D) = (1,3,2,4), a minus sign arises from €13245, at the contrary of 
what happens for U±U l and U^U 3 . The minus sign in \1/}q is put to eliminate 
the minus arising from the e. Clearly there is nothing restrictive in doing this, 
since we are only exploiting the freedom of choosing a convenient phase in the 
definition of our fields. 



C.2 The 75 and 50 representations 

How to write the 75 

The 75 can be conveniently written through the product 

10 x TO = 1 + 24 + 75 . (C.7) 

Let us clarify how to isolate the 75 from the 24 and the 1. We define 10^ x 10c d — 
T^p, so that Tqq is antisymmetric in the upper and lower pair of indices and thus 
corresponds to 100 independent fields. It can be split according to 

rpAB _ qAB i sr^AB , \sAB i p q\ 

1 CD — D CD~r ^CD^~ * CD \^-°) 

1 24 75 

where we identify the 1 with S, the 24 with £ and the 75 with Y. 

Consider the SU(5) singlet first. One can easily guess it to be given by the 
combination Q 

a = T»% . (C.9) 



Thus, S has to be 



s cd = 2q( 5 c s d ~ 5 d s c) t mn ■ ( c - 10 ) 

In this case, in fact, substituting S^d into eq. ( |C.8|) it turns out that S^§ = T^g = s 
and that X and Y must be such that S^ = Y^f = 0. This indeed means that the 
singlet has been completely projected out by S. 

The 24 can be written as a tensor Oq satisfying the constraint o\ = . One can 
guess 

_A — rpMA * frArpMN /p ii\ 

°C = 1 MC ~ T°C 1 MN yy*- 11 ) 



-'^Here and in the following summation over repeated indices has to be understood. 
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and 



J CD 



' s,A(rpMB ^x Br T MN \ sAtrpMB * <;Br r MN\ 

o °C\ 1 MD T°D 1 MN! °D\ 1 MC T°C 1 MN) 

1 1 

rBjrpMA 1 cArpMN\ , sBlrpMA L rArpMN 

°C\ 1 MD T°D 1 MN ) ' °D\ 1 MC T°C 1 MN 



(C.12) 



Notice that this form for E satisfies E^f 



0, that is E doesn't contain the combi- 



nation corresponding to the singlet s. Substituting this expression in eq. ( |C.8| ), we 



get 



ipMA 



cr, 



Ci 



cMA 



and the requirement Y } 



MA 
MC 



0. This means that the 24 is 



completely projected out by E. 

Finally, with the expressions of Sqq and E^£ at our disposal, we know how to 
write the 75: 



Y 



AB 
CD 



rpAB 

1 CD 



1 

5' 



cA/rpMB L cBrpMN\ rAtrpMB L s,BrpMN\ 
°C\ 1 MD T°D 1 MN) °D\ 1 MC T°C 1 MN) 



-lMJV^ 



rB(rpMA rArp]\ 

°C\ 1 MD T°D 1 MN. 



^°D\ 1 MC cr C J J 



20 



(d c d D - d D d c )T 1 



MN 



(C.13) 



This would give 
0, we are 



Note that Y is antisymmetric in both the upper and lower indices 
100 independent fields but, since it satisfies the 25 constraints 1 
correctly left with only 75 independent fields. 

It is important to realize that the fields in Yq^ don't have canonical kinetic 
terms. This can be recognized by the following reasoning. Let us suppose to put the 
100 elements of Yq^ \a<b,c<d into a column vector, say y^, with A = 1, 2, ..., 100. 
Clearly only 75 of these 100 elements correspond to independent fields. The Kahler 
potential is then 



AB*vAB 



21 =? I CD 1 CD 

= iy + y. 



a V AB* V AB I 

4 Y CD Y CD \A<B,C<D 



(C.14) 



We can always exploit the possibility of rotating the fields by means of a 100 by 100 
unitary matrix U. For instance, we can get rid of the unpleasant redundancy by 
going into a basis where y has only 75 non-zero elements. Then, it is also feasible 
to find a unitary transformation of the 75 fields which carry us into a SM-like basis, 
that is a basis where the fields are split into multiplets of the SM. In a similar basis it 
is evident that the kinetic terms are not canonical because of the factor 4 appearing 
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in eq. (lOTiD . The remedy is to define Ygg = Y /A , B /2 (and y = y'/2), so that the 
fields associated to Y' CD are those actually possessing the correct coefficient 1 in 
the kinetic term. 

How to write the 50, 50 

By similar considerations, the 50 can be written from the product of 10 x 10. In 
fact 

10 x 10 = 5 S + 45 a + 50 s . (C.15) 

We denote the 10 x 10 with T AB T^ D , where T\^ are antisymmetric under exchange 
of the two upper indices. Again we want to separate T AB T^ D as 

T AB T CD = C ABCD + qABCD + pABCD ^ ^ 

5 S ^ 50 s 

It is useful to split T AB T^ into its symmetric and antisymmetric parts under the 
exchange 1 <-> 2, so that 

rpABrpCD I rpABrpCD rpABrpCD _ rpABrpCD 

rpABrpCD _ J l 1 2 ~r- f 2 1 1 , 1 1 1 2 1 2 1 1 



2 2 

V ' * -s- 



5 S ,50 S 45a 

= s ABCD + A ABCD . (C.17) 

It follows that Q is antisymmetric while C, F are symmetric under the exchange 
(AB) *-> (CD) (or equivalently 1 «2). Thus, clearly the 45 a , q ABCD ; is given by 

q ABCD = A ABCD ^ (c>lg) 

and the simple identification Q ABCD = A ABCD holds because {Q ABCD -Q CDAB )/2 = 
q ABCD while ( C abcd _ c CDABy 2 = (F ABCD - F CDAB )/2 = 0. Then, we have only 
to separate the 5 S from the 50 s . One can guess that the 5 S is given by 

c E = e ABCDE S ABCD (C.19) 

and verify that the following expression for C ABCD 

qABCD = ^ABCDE^^gFGHI {Q2Q) 

satisfies €abcdeC abcd = c E - Since €abcdeQ ABCD = 0, C completely projects out 
the 5 if also €abcdeF abcd = holds. This property is automatically satisfied by 

pABCD rpABrpCD 



_ aABCD 
1 



-e 



ABCDE^ qFGHI 



24 



£fghieS . (C.21) 
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Note that F is antisymmetric in [AB] and [CD] but is symmetric for the exchange 
(AB) <-► (CD). 

By writing the Kahler potential as K 3 f abcd *F abcd , it turns out that F ABCD 
hasn't canonical kinetic terms. The proof is very similar to the one already discussed 
for the case of the 75. Also in the present case, defining F ABCD = F' /2, one 

realizes that the fields F' are those with the kinetic terms in the canonical 

form. 

The 50 is obtained following the same procedure but starting from the product 
of \Qx\Q = Y lAB Y 2CD . 



C.3 The Method for Catching SM Multiplets 

The construction of the 75, 50 and 50 in terms of products of 10 and/or 10 is 
particularly useful because one can apply a simple recipe to recognize how the SM 
multiplets are included in these large representations. The starting point is the 
identification of the SM multiplets in the 10 and 10 themselves. Keeping A < B: 



(C.22) 



Let us consider the case of the 75 (50) in order to be specific. One starts by writing 
the SM multiplets obtained by multiplying each SM multiplet of T with each SM 
multiplet of T (T). All the 100 fields obtained are thus separated in SM multiplets. 
Next we have to further correctly separate the various multiplets into the 1, 24 and 75 
(5 S , 45 a and 50 s ). This is easy to do for such multiplets which are obtained only once, 
because from Table |C1| we immediately understand in which SU(5) representation 
they lie. But for the multiplets which are obtained more then once, like e.g. the 
singlet (the color antitriplet), we have to understand which linear combination of 
them is the one actually embedded in the 75 (50). As an example, we show the 
procedure which permits to single out the SM singlet contained in the 75 and the 
SM antitriplet contained in the 50. 

The SM singlet of the 75 

There are three singlets (1, 1, 0) arising from the product T x T. 



rpAB 


= 10 a/3 + 10 afc + 10 afc 




(3,1,-2/3) (3,2,1/6) (1,1,1) 


Tab - 


= 10 a/3 + 10 Qfe + 10 ab 




(3,1,2/3) (3,2,-1/6) (1,1,-1) 
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SU(5) D SU(3) x SU(2) x 17(1) 

1 = (1,1,0) 

24 = (1,1,0) + (1,3,0) + (3, 2, -5/6) + (3, 2, 5/6) + (8, 1,0) 

50 = (1,1, 2) + (3, 1,1/3) + (3, 2, 7/6) + (6, 3, 1/3) + (6, 1,-4/3) + 
+ (8,2,-1/2) 

50 = (1,1, -2) + (3, 1,-1/3) + (3, 2 -7/6) + (6, 3, -1/3) + (6, 1,4/3) + 
+ (8,2,1/2) 

75 = (1,1,0) + (3, 1,5/3) + (3, 1,-5/3) + (3, 2, -5/6) + (3, 2, 5/6) + 
+ (6, 2, -5/6) + (6, 2, 5/6) + (8, 1, 0) + (8, 3, 0) 

Table C.l: The SM multiplets contained in some of the SU(5) representations. 
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We call Si the singlet which is obtained by multiplying 

(3.1.-2/3) (3^3) V ^ ' 

+ T^-hs^-8^)T^ . (C.23) 

N v ' 

(8,1,0) 

It is then evident that si oc Tj^ = 2T^\ IX<U . Thus, choosing to normalize to 1 each 
combination of independent fields: 

*i = -^(Tll + Tl! + T%) = ^fHX<» . (C.24) 

The second singlet s 2 comes from the product 

rpab __ rpab ~rp X a xbrpij,n 

-yd — v^^ L id — „ "7 °d 1 ixn 

(3,2,1/6) (o o _i /fil v 

V v ' 

(1,3,0) 

+ ^X:~Wd T P (C25) 



2 a 7 " 6 

(8,1,0) 



1 1 1 

i T-iafc sctrpixb rbrpan , A"l' , 7 1 /"' 









(8,3,0) 




Thus, s 2 oc T^™. Precisely: 










s -7! (r " + T 


24 j_ rpM _, rpl5 

24 "T J 34 i- -1 15 


I rp25 

"T" J 25 


+ ^11) = 


rpfm 


Finally the last singlet 


S3 is obtained from 






rpab 


-^T^ T^ = 

(1.1,1) (1,1,-1) 


^(«-% 6 )' 


-pmn 
mn 






(1,1,0) 




and is simply 


S3 - J 45 - 


mn 1 


m<n 





(C.26) 



(C.27) 



(C.2* 
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The next step is to understand which linear combination of si.,2,3 is contained in 
the 75. This can be accomplished by using eq. ( |C.8| ) and the expressions for S, S, Y, 
that is eqs. flC.10| , |C.12| , |C.13|) . Indeed, by constructing on the l.h.s. of eq. ( |C.8| ) 
the combination relative to si.2,3, we obtain respectively: 



Sl 



s 2 = 



rpa/3 
1 af3 

rpab 
1 ab _ 




10 
6 



Sl 




T 



S3 



ab 

Y 




Si 



Si 



+ 



+ 




S> 



21 



+ 



s- 



75 



Y 




S- 



21 




15 



S> 



21 



+ 



+ 



s> 



75 



Y 



-s> 



75 



where 



^1 = 



5*24 = 



5*75 = 



a/3si + \ZQs 2 + s 3 



10 



a/8si - s 2 - VGs 3 



15 



si - V2s 2 + V3s 3 



(C.29) 



(C.30) 



It is clear that by taking in the l.h.s. of eq. ( |C.8| ) the combination S75 , one obtains 
5*75 — [0]s + [0]s + [5 , 75]y- The SM siglet lying in the 75 is then S 75 . One can check 
that the combinations Si, S'24 are instead respectively the singlets of the 1 and 24 
representations. Note that S 75 can be written as 



S- 



75 



3V2 



rpCXp , cyrpab rpAB 

1 ' + Zi ab ~ 9 AB 



-a/3 



(C.31) 



The antitriplet of the 50 

There are three multiplets in the product 10 x 10 having the quantum numbers of 
the antitriplets, that is (3, 1, 1/3). Two arise from the products 



T[ 



ab 



T, 



7<5 



t 



ab-yS 



(1,1,1) (3,1,-2/3) 



and 



na/3 



-•cd 



(3,1,-2/3) (1,1,1) 



rpca __ j. 

J-2 — l 2 



aftcd 



(C.32) 



(C.33) 
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The other one is contained in 



rpab n-r/d 
(3,2,1/6) (3,2,1/6) 



(C.34) 



and can be written as 



/rpabrpjd , rrrfdrpab rrnadrrryb rp-ybrpcid\ 

.1. — L -_ _L — J. — 



,ab~fd 



(C.35) 



Following the same procedure already described in some detail for the extraction 
of the SM singlet of the 75, we take on the l.h.s. of eq. QU.lfip the combination 
corresponding to each antitriplet: 



,abfS 



.■ySab 



,ja8b 



6 ( § ) 



- rpab'yS 

6 < 5 > 



c 



c 



_ rpab"/S 

3 © 



rpab-yS 
2 (^ 



+ 



c 



+ 



Q 



-WrT! 



-iafe7<5 



2 ( 4 5) 



-*- rpab^S 
3 i (50) 



rpab-yS 
3 i (50) 



+ 



ab-yS 



3 T (50) 



(C.36) 



where 



rpab-yS 
J (5) 



rpab-yS 
(45) 



rpab-yS 
i (50) 



f 



0675 . ,-ySab 



+u 



Vq 



t 



ab-yS i-ySab 



t 



V2 



t 



ab-yS . ,-ySab 



It 



-yaSb 



+ t-T" + £ 



-yaSb 



Vs 



(C.37) 



Thus, one realizes that the linear combination of the above antitriplets which is 
contained in the 50 is TZ2 . Note that, as required, it is symmetric under the 



1/4 e a7 se a bT^2 S . Note that, 



exchange T\ <-> T^. 

To get rid of the indices it is useful to define £(5o) a — -v -* ■-cryd'-ao-t (5,, 
as already discussed, £( 50 ) hasn't canonical kinetic terms. It is H' 3da (introduced 
in section |4.5|) , defined as £(50) = H 3da /2, that possesses the correct coefficient 1: 



K 3 F ABCD* pABCD 3 ^ - 



(50) a «-(50) c 



I1 3da rl 3da- 
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C.4 Basic Formulae for the MDM 

SU(5) breaking 

The part of the superpotential responsible for the breaking of 577(5) reads 

WHiggs = * Y$E y eF y aE + My y£> Y™ , (C.38) 

where, as usual, A,B, ... = 1,2, ...,5. This is actually the most general renormaliz- 
able superpotential that can be written by means of a 75, the cubic term picking the 
combination (75 x 75)75 x 75. The other cubic term that one could write, namely 
y ed y fc y ab, turns out to be proportional to Ygg Ygf Yfg. 

The SUSY vacuum of the 75 breaks SU(5) uniquely to SU{3) <g> SU{2) U(l): 



(Y C A B) = -00 



(s*Sl + 2S a J b d - is&i) ~{A~B) 



(C.39) 



where, as usual, Greek indices assume values from 1 to 3, Latin indices from 4 to 5 
and 

(Y)=M Y /{2 Cl ) . (C.40) 

The above superpotential gives a mass to all physical components of Y, i.e. those 
that are not absorbed by the Higgs mechanism. The calculation is not difficult and 



leads to Table 4.2. 



Gauge vector bosons instead acquire a mass m v : 

m 2 v = 24g 2 (Y) 2 . (C.41) 

The triplet mass matrix 

The doublet-triplet splitting is realized by the following superpotential terms, 

tAv-DEtt BCFG^ , „ T7 \/BCtj ADEFG 



WMDM — C 2 H Y bc H 50 tADEFG + C3H aY de H-^BCFG € 

^■Wabcd 



+ M 50 H 50 H 50ABCD ■ (C42) 



The multiplet in Y which gets a vev is the SM singlet. This means that in the 
term 

c 2 H Y BC H 5Q cadefg (C.43) 

only the interaction between the triplet contained in H A and the antitriplet con- 
tained in 7/50 is allowed. Since H^ Q BCD can be identified with F ABCD ; using eqs. 
( |C.39|) and ( |G.36| ), one obtains 



c 2 H A (Y B D c E )H 50 BCFG e ADEFG = -4V3c 2 (r> H' 3d H 3u . (C.44) 
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Finally, 



TT ABCDtt_ 

-"50 n bQABCD ^ x± 3d ±x Zu 



9 H 3d H k 
so that the coefficients appearing in eq. ( ^.9| ) are obtained. 



(C.45) 



Interaction of 10 10 50 

We are interested in the interactions of the colour triplet of the 50 with the an- 
titriplet combination of 10j x 10j, 10 being matter fields and i, j family indices. It 
is convenient to recognize that 

1 — ^10 i ^Wij ^10 j HMABCD ^ 4 G~50ij t(50)aij ^50) , (C.46) 



where by t^ ai j we denote the antitriplet combination of 10j x 10j, that is t(so)aij 

pab-yS 
(50)ij : 



1 I A rpalryS 

1/4 Ea^ygdaf)! n 



rpab'yS 
(50)ij 



,r r ab mt6 , , t ,7* \T,ab , W 10i W 10j + ^lOi^lOj ^lOi^lOj *10i*10j 

w i0i w i0j "•" w i0i w i0j "i 2 

(C.47) 



Substituting the above expression in eq. (|C.46|) and exploiting the symmetry of GgQ, 
one obtains 



jjc 1 2( ^50 ij 



E i--oQ 



\Q?(-^)Qle 



ap-y 



m 



3« 



This is the cause of the introduction of the coefficient Vo/2 in eq. Q4.24Q . 



(C.48) 



